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XXII. BAND ZWEITES HEFT 1954 


Der selbsterregte unsymmetrische Kreisel. 


Von R. Grammel. 


1. Einleitung. Ein starrer Kérper, der sich um seinen Schwerpunkt frei drehen kann, soll 
als selbsterregter Kreisel bezeichnet werden, wenn seine Antriebsdrehkraft nicht 
von auSen kommt, sondern in ihm selbst erzeugt wird (etwa durch RiickstoBvorrichtungen, 
die aber seine Masse und seine Massenverteilung nicht merklich andern migen). Je nach- 
dem der Momentvektor Jt einer solchen Antriebskraft eine feste Lage im Kreisel hat oder sich 
im Kreisel auf vorgeschriebene Art bewegt, spricht man von achsenfester oder 
achsenbeweglicher Selbsterregung; und bei achsenfester Selbsterregung wird man 
dann noch die beiden Falle unterscheiden miissen, daB® der Betrag von St zeitlich unverdnder- 
lich oder veranderlich ist: zeitunabhaingige und zeitabhangige achsenfeste 
Selbsterregung. 

Wahrend nun die Bewegung des achsenfest selbsterregten symmetrischen Krei- 
sels sich durch Quadraturen ausdriicken ]4Bt 1, weiB man iiber den selbsterregten uns ym - 
metrischen Kreisel bisher nur wenig. Zwar liefern die Eulerschen Gleichungen zu einem 
Drehvektor 0, welcher im Kreisel festliegt oder in ihm auf vorgeschriebene Weise wandert, das 
(i. a. dann weder zeitunabhangige noch achsenfeste) Moment St der erforderlichen Drehkraft ; 
aber die umgekehrte Fragestellung, zu einem im Kreisel vorgegebenen Erregungsmoment Jt 
den Drehvektor 0 zu ermitteln (und daraus vollends in allerdings dann bekannter Weise etwa 
die Eulerschen Lagewinkel des Kreisels im Raume), bereitet erhebliche Schwierigkeiten, weil 
die Eulerschen Gleichungen den Drehvektor 0 nicht linear enthalten. Schon der klassische 
Euler- Poinsotsche Fall ))¢ = 0 erfordert bekanntlich elliptische Funktionen; die allgemeine 
Lésung fiir Jt =-0 ist meines Wissens bislang nicht einmal fiir zeitunabhangige achsenfeste 
Selbsterregung, geschweige denn fiir zeitabhangige achsenfeste oder fiir achsenbewegliche 
Selbsterregung gefunden worden, — ausgenommen den trivialen Fall des Antriebs um eine 
raumfeste Hauptachse des Tragheitsellipsoids des Kreisels beziiglich seines Schwerpunkts, und 
weiter ausgenommen die Gesamtheit aller permanenten Drehungen? des unsymme- 
trischen Kreisels, die durch ein zeitunabhangiges achsenfestes Moment zwar nicht erzeugt, aber 
wenigstens unterhalten werden (man kennt auch das Stabilitatsverhalten solcher permanenten 
Drehungen und weiB, wie ihre Stabilitat notfalls durch geeignete Steuerung der Selbst- 
erregung erzwungen werden kann). 

Im folgenden werden wir die allgemeine Lisung fiir den zeitunabhangig und zeitabhangig 
achsenfest und den achsenbeweglich selbsterregten unsymmetrischen Kreisel zwar noch nicht 
vollstandig erhalten; wir werden ihr jedoch sukzessiv sehr nahe kommen, sie innerhalb ge- 
wisser Konvergenzbereiche erreichen und in einigen Sonderfallen sogar exakt darstellen. 
Zuerst behandeln wir mit einer iterativen Methode den Fall, daB der Momentvektor Jt 
der erregenden Drehkraft in der kleinsten oder in der gré8ten Hauptachse des Tragheits- 
ellipsoids beziiglich des Schwerpunkts liegt und der anfangliche Drehvektor 0, nur eine kleine 
Komponente senkrecht zu jener Achse besitzt (Ziff.2). Dann kommt der Fall, daB der 

~Momentvektor Jt nur schwach von einer jenex beiden Achsen abweicht und der anfangliche 
Drehvektor entweder in jener Achse liegt oder Null ist (Ziff. 3). Eine geeignete Kombination 
dieser beiden Falle, die aber wegen des nichtlinearen Geprages der Eulerschen Gleichungen 
keineswegs eine additive Uberlagerung sein kann, liefert innerhalb gewisser Konvergenzbereiche 
die allgemeine Lésung fiir zeitunabhangig und -abhangig achsenfeste sowie fiir achsenbeweg- 
liche Selbsterregung (Ziff. 4). Alsdann nehmen wir das Problem von Ziff. 2 erneut vor und 


1 Fiir den zeitunabhangig achsenfest selbsterregten symmetrischen Kreisel hat U. T. Bédewadt, 
Math. Z. 55 (1952), S. 310, die Lésung angegeben; fiir den zeitabhangig achsenfest selbsterregten sym- 
metrischen Kreisel kann die Lésung in gleicher Weise gefunden werden, nur liegen dann die dabei auf- 
tretenden Integralfunktionen i. a. noch nicht tabuliert vor. 


2 R. Grammel, Ing.-Arch, 21 (1953), S. 149. 
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lésen es in exakter Form (Ziff. 5), und zwar auch fiir einen Momenivektor S)¢ in der mittleren 
Hauptachse. Wir werden dabei drei Bewegungstypen kennenlernen: einen periodischen und 
zwei asymptotische Typen. Da weder diese Lésung noch der Konvergenzbereich der Lésung von 
Ziff. 4 bis zu den beiden Ebenen der sogenannten trennenden Polbahnen des (zum Tragheits- 
ellipsoid koaxialen) Poinsotellipsoids heranreichen, so behandeln wir weiter den Fall, daf der 
Drehvektor 0 in diesen Ebenen infolge eines zeitunabhangig achsenfesten Momentes S{ wan- 
dert, welches entweder in der mittleren Hauptachse des Tragheitsellipsoids liegt (Ziff. 6) oder 
in einer bestimmten Weise beliebig stark von ihr abweicht (Ziff. 7). SchlieBlich beweisen wir 
noch drei spezielle Satze: namlich erstens, daB es (abgesehen von den trivialen Fallen der 
drei Hauptachsen) keine zeitunabhangige achsenfeste Selbsterregung gibt, welche eine achsen- 
feste Drehung des Kreisels erzeugen oder vernichten oder beliebig beeinflussen kénnte (Ziff. 8) ; 
zweitens, dai es (wieder abgesehen von jenen trivialen Fallen) nur eine in ganz bestimmter 
Weise von der Zeit abhangige achsenfeste Selbsterregung gibt, welche eine zeitabhangige 
achsenfeste Drehung in ihrer Geschwindigkeit beeinflussen kann (Ziff. 9); und drittens, daB 
zur Erzeugung einer achsenfesten Drehung (wieder abgesehen von jenen trivialen Fallen) 
eine achsenbewegliche Selbsterregung nétig ist, deren Momentvektor Yt dabei jedoch stets 
in einer kérperfesten Ebene bleibt (Ziff. 10). 


Die fiir den selbsterregten unsymmetrischen Kreisel zustandigen Eulerschen Gleichungen 
lauten 
Aw, — (B— C)o,o, = M,, | 
Bo, +(A—C)o,0, = M,, ! (1) 
Cw, —(A— B)o,wy = M., | 


worin A, B, C die Hauptdrehmassen (Haupttragheitsmomente), w,, wy, @, die Komponenten 
des Drehvektors 0, ferner M,, My, M, diejenigen des Momentvektors Jt der erregenden Dreh- 
kraft, je in dem korperfesten, rechtshandigen Hauptachsenkreuz (x, y, 2) durch den Schwer- 
punkt, bedeuten. Wir schreiben die Gleichungen (1) mit der Rangordnung 


Ay> BEG (2) 
und den somit positiven Tragheitszahlen 
Prec Meare Vee 
und den drei MomentgréBen 
Tae a, Oe a (4) 
zweckmafig in der Form 
Dx — 1 Wy OW, —= My , (5) 
Dea p On, == ity (6) 
WO, — Y Wx Wy = M,. (7) 


Zeitunabhangige achsenfeste Selbsterregung liegt vor, wenn M,, My, M, und also auch m,, 
my, m, Festwerte sind, zeitabhangige achsenfeste Selbsterregung, wenn mit einer gemeinsamen 
Zeitfunktion F(t) und drei Festwerten M,, M,, M; 

M, = M, Fit), M, = M, F(), Me aM EG) (8) 
ist, endlich achsenbewegliche Selbsterregung, wenn M,(t), M,(t), M,(t) und also auch m,(t), — 
m,(t), m,(t) drei nicht zu einander proportionale Zeitfunktionen sind. | 


Die (beschleunigten oder verzégerten) Drehungen des Kreisels um eine seiner drei Haupt- | 
achsen infolge eines Momentvektors St in der gleichen Hauptachse wollen wir weiterhin ein | 
fiir alle Mal als triviale Falle ausschlieBen. 


2. Momentyektor in einer Hauptachse; iterative Lésung. Der erste Fall, den wir unter- | 
suchen wollen, ist gekennzeichnet durch die Voraussetzungen 


ity KOUSL ee Om my, =0, ie, =D) . (9) | 
das heiBt: der achsenfeste Momentvektor JM soll in der Hauptachse der gréBten 
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Drehmasse A (x-Achse) liegen und zunachst zeitunabhangig sein. Die Voraussetzung, 
da8 m, positiv sein soll, ist keine Einschrankung der Allgemeinheit, da man die positive 
Richtung der x-Achse stets entsprechend wahlen kann. Von den Anfangswerten xo, Dy, On 
der Komponenten des Drehvektors zur Zeit t =0 darf dann q@,, positiv, Null oder epee 
sein. Damit nicht der triviale Fall vorliegt, miissen wir aber voraussetzen, daB von den beiden 
anderen Anfangskomponenten @,. und @,) mindestens die eine von Null verschieden sei. Und 
nun ist die wesentliche Voraussetzung, die wir machen miissen, um die Konvergenz des fol- 
genden Verfahrens der sukzessiven Approximation zu erméglichen, 


Wyy <M, Ge, aM, . (10) 
(Man stellt leicht fest, daB w} und w} dimensionsgleich mit m, sind.) 

Unsere Voraussetzungen besagen, da ein um seine erste Hauptachse (x-Achse) durch 
.Selbsterregung dauernd angetriebener Kreisel zur Zeit t = 0 einen schwachen seitlichen Dreh- 
stof erhalt, der seine Drehachse ein wenig aus der ersten Hauptachse herauswirft, genauer 
gesagt (wenn man auf die Bedeutung (4) von m, achtet, und die drei DrehmassenA, B, C 
als von gleicher GréSenordnung ansieht), da die zusitzliche Bewegungsenergie, die der seit- 
liche DrehstoS dem Kreisel erteilt, klein gegen das (damit dimensionsgleiche) Selbsterregungs- 
moment bleibt. Wir fragen nach dem weiteren Verhalten des Drehvektors 0. 

Wir werden bestatigen, dai infolge der Voraussetzungen (10) das Produkt |wy o,| 
dauernd klein gegen m, bleibt. Unterdriicken wir daher dieses Produkt fiir die erste 
Naherung in (5), so folgt durch Integration 


oD = Wy + m,t, (11) 

und damit gehen die Gleichungen (6) und (7) wegen (9) iiber in 
oy + B (@x9 + ms t) aa 0, | 
() . 


Os I (Og 2m, t) OP = 0°, J 


wobei der obere Zeiger bei w,, wy, w; hier und weiterhin jeweils den Grad der Naherung aus- 
driicken soll. 

Man bestatigt leicht, daB die allgemeinen Integrale dieses gekoppelten linearen Systems (12) 
mit zeitabhangigen Koeffizienten die Gestalt haben 


ao = Wy VB cos [VB y (Wag t + 4 m,t?) + o| : 


(12) 


= tone (13 
wo — a, Vy sin [VB y (xt + 4 m, t*) +o] , | ) 
wobei die Konstanten @, und o durch 
w? ow? 
Ff Ls ae oF 14 
O; B ee y (14) 
Ke Wzo 
te oO = — — LS 
g ar oe (15) 


definiert sind, und wobei voraussetzungsgem4B w,-++ 0 ist und die Quadratwurzeln in (13) 
und (15) als positiv festgesetzt werden diirfen. 

Man bestatigt, daB zufolge der Voraussetzung (10) das Produkt jot) o)| in der Tat klein 
gegen m, ist, so daB seine Vernachlassigung in (5) fiir die erste Naherung berechtigt war. 

Die Deutung dieser ersten Naherung ist einfach. Wahrend die x-Komponente ow) (11) des 
Drehvektors 0 gleichférmig mit der Zeit anwachst, so beschreibt die Pfeilspitze des zusatz- 
lichen (kleinen) Drehvektors 0* mit den Komponenten ies cw) (13) in der (y,z)-Ebene immer 
schneller eine Ellipse mit den festen Halbachsen w, JB und a, Vy. Mithin beschreibt der ganze 
Drehvektor 0 immer schneller einen elliptisch-spiraligen Kegel um die x-Achse mit anwachsen- 
den Mantellinien, derart, da® die Pfeilspitze von 0 auf einem elliptischen Zylinder, dessen 
Achse die x-Achse ist, und dessen Halbachsen die angegebenen Langen haben, eine raumliche 
Spirale mit monoton anwachsenden x-Werten aufzeichnet. Demgem48 nahert sich der Winkel 
zwischen Drehvektor 0 und x-Achse mit immer rascher werdenden Schwankungen asymptotisch 
dem Wert Null. 

Damit ist zugleich gezeigt, daB die Hauptachse der gréSten Drehmasse A, deren Stabilitat 
meines Wissens bisher nur fiir die permanenten Drehungen des kraftefreien unsym- 

6* 
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metrischen Kreisels bewiesen wordenist, auchhbei beschleunigtem selbsterreg- 
tem Antrieb stabil gegen kleine seitliche Drehsté8e bleibt. 

Wenn man sich an das Poinsotsche Bild der Bewegung des kraftefreien unsymmetrischen 
Kreisels erinnert 1, so erkennt man leicht, daB die Bewegung des Kreisels, von auBen her 
betrachtet, so verlauft: der Kreisel dreht sich beschleunigt um seine x-Achse, und diese um- 
tanzt immer schneller und immer enger ihre urspriingliche Raumlage vor der Stérung und 
nihert sich dieser nach der Stérung wieder asymptotisch. 

Man kann dieses Bild der Bewegung noch verfeinern, wenn man von der ersten zu einer 
zweiten Naherung ibergeht. Hierzu setzt man die gefundenen Werte (13) von On 
und win die Gleichung (5) ein und erh4lt so statt (11) 


eo! = wa» +mt +l @raypy S® © ae) 


mit der Integralfunktion 
t 


QM) (t) = | sin 2 [VB y (mut + 4m, #) to] dt, (17) 


deren Auswertung wir spater nachholen werden. Hiernach gehen die Gleichungen (6) und (7) 
wegen (9) iiber in 
op Ba. o, — Ol i 
2? — y w@?) wo?) = f (18) 
mit den al] gemeinen Integralen 


cof) =o, Boos [VB 7 K @) +0]. | 
o =o Vy sin [JB 7 K® () +0]. J 


welche die weitere Integralfunktion 


(19) 


t t | 
K® (t) = f w® dt = ont +4 m2 +1 08a VB y [ 3% @ at (20) | 
0 0 


enthalten. 
Die dritte Naherung, in gleicher Weise aus der zweiten hergeleitet, lautet 


Of) =O» + met + FopayBy S(t), 


col) = 1 VB cos [VB y K(9) (¢) + o| A (21) | 
af =o, Jy sin [YBy KO) +o] J 
mit 
3® (t) = [sin 2 [PB y K® (t) +o] dt, (22) 
0 
K® (t) = [ of dt =a ot +4im,2 +102 a/B y [ 3° (ode, (23) 
0 0 


und man sieht nun ohne weiteres, wie die Lésung durch gleichartig fortschreitende Naherungen 
mehr und mehr verfeinert werden kann. Wesentlich ist dabei, daf® in allen diesen Naherungen 
die Konstanten @, und o ihre etwa vorgeschriebenen Werte (14) und (15) nicht andern. 

Auch an dem Verlauf der Bewegung andert sich gegeniiber der ersten Naherung nicht mehr 
viel. Der elliptische Zylinder, der die spiralige Bewegung des anwachsenden Drehvektors 9 
begrenzt, bleibt véllig erhalten, und nur das, was wir vorhin als zeitliches Schwanken bezeich- | 
net haben, wird durch die weiteren Naherungen verfeinert wiedergegeben: dieses Schwanken 
wird nun schon in der Hauptkomponente w, des Drehvektors 0 beriicksichtigt, und zwar immer | 
genauer mit jeder weiteren Naherung, und iibertragt sich dann natiirlich auch auf den zeit- | 
lichen Ablauf der Nebenkomponenten w, und @, von 0. Der elliptische Zylinder hat iibrigens | 
das gleiche Halbachsenverhaltnis VBly = yc (A — C)/B (A — B) wie schon beim kraftefreien | 


Kreisel, bei welchem er die Polbahn, die die x-Achse umschlieBt, aus dem Poinsotellipsoid | 
ausschneidet 2, 


1 Siehe etwa R. Grammel, Der Kreisel Bd. 1, S. 125, 2. Aufl. Berlin-Gétti -Hei 
2 Siehe R.Grammel, Der Kreisel Bd. 1, S. 131. anh ean erence aS 
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Jetzt ist noch nachzutragen, daB wenigstens die in der zweiten Naherung auftretende Inte- 
gralfunktion 3) (t) (17) sich auf bekannte Funktionen zuriickfithren ]aBt, némlich auf die 
Fresnelschen Integrale, welche in der Form 


3 
aL cos & 
V2a J Vé 
0 


tabuliert sind’. Man transformiert $((t) auf eine neue (dimensionslose) Zeitvariable € und 
ihren Nullwert &: 


g . 
C(é) dé, Sf) = — | ue ae (24) 


E are. Mxg 2 5 a wr 
=m, VBy(6+28),  & = B72 (5) 
und hat dann zunichst 
é 
Y(1) (7) — sin (¢ Sot 29) 4 26 
ee (26) 


mit 
: = tS (27) 
2 mx\/B y 
und also schlieBlich 


S(t) = RALS (6) — S (£)] cos (E&) — 2 0) — [€ (£) — €(&)] sin (& — 2 0)} . (28) 
Auf die mégliche Auswertung auch der Funktion K®) (20) wollen wir hier nicht weiter eingehen?. 

In genau der gleichen Weise kann man den Fall durchrechnen, daB der Momentvektor Jt 
in der Hauptachse der kleinsten Drehmasse C (z-Achse) liegt, und daB der Drehvektor 
anfanglich eine kleine Komponente senkrecht zu dieser Achse besitzt. Offenbar braucht man 
in den bisherigen Ergebnissen lediglich die Buchstaben entsprechend zu vertauschen. AuBer- 
dem ist so dargetan, daB die Hauptachse der kleinsten Drehmasse C auch bei beschleunigtem 
selbsterregtem Antrieb ihre vom kraftefreien Kreisel her bekannte Stabilitat beibehalt. 

Dagegen laBt sich der Fall, dafS der Momentvektor Jt in der Achse der mittleren 
Hauptdrehmasse B (y-Achse) liegt und der Drehvektor 0 anfanglich eine kleine Komponente 
senkrecht zu dieser Achse besitzt, nicht auf diese Art behandeln, weil dann, wie man leicht 
erkennt, die (beschrankten) Kreisfunktionen cos und sin durch die unbeschrankt anwachsen- 
den Hyperbelfunktionen Wpj und Gin zu ersetzen waren, im Gegensatz zu der Annahme (auf 
der die Rechnung doch beruhen wiirde), daB die anfanglich kleinen Komponenten @, und w, 
in der Gleichung (6) dauernd klein bleiben sollen. Hiermit ist zugleich die Erkenntnis gewonnen, 
daB die mittlere Hauptachse, die schon beim kraftefreien unsymmetrischen Kreisel] eine labile 
Drehachse vorstellt, auch bei beschleunigtem selbsterregtem Antrieb labil bleibt. 

Wie weit sich der Konvergenzbereich unserer Rechnungen fiir die beiden Falle erstreckt, 
dai der Momentvektor $f in der x- bzw. z-Achse liegt, d.h. fiir welche Werte @yo, Wzy D2W. Wx, Wyo 
unsere sukzessiven Naherungen gegen die exakte Lésung konvergieren, lasse ich hier uner6rtert. 
Der Gang der Rechnung lat vermuten, daB dieser Konvergenzbereich ziemlich groB sein wird, 
weil jeder weitere Naherungsschritt immer nur Verbesserungen von der GréBenordnung wf 
aufzeigt. 

Formal 1aBt sich die ganze Rechnung sofort auf den viel allgemeineren Fall der zeitab - 
hangigen achsenfesten Selbsterregung erweitern, wobei also m, nicht mehr ein 
Festwert, sondern eine vorgeschriebene Funktion der Zeit ist. Offenbar hat man dann iiberall 
die Ausdriicke m, t und } m, t” zu ersetzen durch die Integrale 


| m,(t) dt , ie m,(t) dt? , 


und ganz Entsprechendes gilt fiir ein zeitabhangiges m,(t). Aber natiirlich muf man nun 
beziiglich der Konvergenz iiberhaupt und insbesondere beziiglich einer guten Konvergenz in 
der Auswahl der Funktionen m,(t) baw. m,(t) vorsichtig sein. 


1 Jahnke-Emde, Tafeln héherer Funktionen S. 36, 4. Aufl. Leipzig 1948. 
2 Vgl. W. Magnus u. F. Oberhettinger, Formeln und Satze fiir die speziellen Funktionen der mathe- 


matischen Physik, S. 127, 2. Aufl. Berlin-Géttingen-Heidelberg 1948. 
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Wir werden auf das hier durch Iteration geléste Problem wenigstens fiir zeitunabhangige 
Selbsterregung in Ziff. 5 noch einmal zuriickkommen und dann auch noch eine geschlossene 
exakte Lésung aufsuchen; und diese wird zeigen, da® die Mannigfaltigkeit der Bewegungs- 
formen gréBer ist, als die [an die Voraussetzungen (10) gebundene] iterative Lésung angeben 
konunte. 


3. Momentyektor nahezu in einer Hauptachse. Wahrend bisher der Momentvektor {t 
genau in einer Hauptachse lag (und zwar entweder in der x- oder 2-Achse, aber nicht in der 
y-Achse), der anfangliche Drehvektor 0, aber von ihr abweichen durfte, so wenden wir uns nun 
dem Gegenstiick zu, nimlich dem Falle, daB der anfangliche Drehvektor 0) genau in einer 
Hauptachse liegt (einschlieBlich des Sonderfalles, da er Null ist), zunachst in der x-Achse, 
also mit den Komponenten @xo, 0, 0 (einschlieBlich des Sonderwertes Wy = 0), daB dagegen 
der achsenfeste Momentvektor, zunachst zeitunabhangig, von der x-Achse 
schwach abweicht. Wir bezeichnen demgem48 die drei MomentgréBen (4) mit m, und 


aD Is —=1 ll, We (29) 


wobei fy und yz hinreichend kleine Zahlen sein sollen. Wir wollen und diirfen m, wieder als | 


positiv voraussetzen, wogegen fu, und jl, positiv oder negativ sein kénnen (aber nicht beide 
Null sein mégen, damit nicht der triviale Fall vorliegt). 


Auch hier wird sich die Vermutung, daB® alsdann |wy w:| dauernd klein gegen m, bleibt, | 
durch die weitere Rechnung bestatigen, und somit unterdriicken wir fiir die erste Nahe- | 


rung wiederum das Produkt w,@, in (5) und erhalten fiir w) das frithere Integral (11), 
so daB die Gleichungen (6) und (7) mit (29) iibergehen in 


io +B (my + me t) of = py my, | | 
a) ee a oh ( (30) | 
Oz” Y (@xo My t) Wy == gn) 
Thre an die Anfangsbedingungen 
; wy =0, 0, =0 fir 1=0 (31) 
angepaBten Integrale lauten, wie man durch Einsetzen bestatigt, | 
a , . 
Oy) =m, IM Se (i) + Hs yz Oe ( cos Ve y (oro ono Les ) is 
++ m, Ms SP) (t) — Ms ye Se? (t)| sin \VBy (or t i My | > 
ios LS : (32) 
of) = — mele) % 30) — 43 920] 005 17 (oat +4m6#)] + 
+ mela |/% 920 +9 990] sin [VF 7 (Ot +m) 
mit den Integralfunktionen 
t 
SO (2) =| cos VB y (@xg t +4 m, #)] die | | 
= ( (33) | 
32) = | sin [YB 7 (Oxo +4 me P)] dt, J 


welche man in gleicher Weise wie (17) auf die tabulierten Fresnelschen Integrale (24) zuriick- 
fiihren kann. Mit 


/ 1 sae Dx, : / 1 oe 
g =F mVBr( 4-< ) > & = VP y— (34) 


Mx Mx 
und 


4) (t) ye {Lc (é’) — C (&)] cos > + [S (é’) — S (€)] sin Eo} : 
Os” (1) = F’ {[S (€) — S (&)] cos & — [C (é’) — C(&)] sin &} . 


findet man 


| 
( (36) 
J 


Weil die Fresnelschen Integrale fitr alle Argumente nicht nur beschrankt sind, sondern sogar | 
kleiner als 1 bleiben und sich mit wachsendem Argument in immer kleiner werdenden Schwan- | 
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kungen dem Wert 1/2 nihern, so bleiben auch || und |o| gemaf (32) und (36) beschrankt 
und wegen der kleinen Faktoren w. und yz dauernd klein, und zwar bleibt ihr Produkt in der 


Tat klein gegen m,, wie man erkennt, wenn man auf den Faktor k’ ~ 1/y/m, in 3 und 3 (36) 
achtet. 


Wenn man (32) in der einfachen Form 


oy) = mg [e, cos y (t) +e, sin g (t)], 


ow) = m, 2 [c, sin ~ (t) — c, cos p (t) | 


schreibt, mit den Abkiirzungen 


P()=VBY (rot +3 meP), (37) 
a(t) = 320) +o /2 9, 
- (38) 
a= 390 —mV2 300 | 
und dann noch 
c, (t) =A(t) cos o'(t), cy (t) = h(t) sin o’(t) 
setzt, so daB 
mah + =(o8 +8 Z) (92 +90), (39) 
io See Nee os (40) 


wird, so nehmen die beiden zusatzlichen Drehkomponenten wy) und w, die zu wY hinzu- 
kommen, die Gestalt an 


ol) = m,h cos (py —o’), | 


ao” = m, nye sin (y — 0’) | 


wobei k(t) gemaB (39) eine Funktion von klein bleibendem Betrag ist. 

Demgem48 iiberstreicht der zusatzliche Drehvektor p* mit den Komponenten Oy) wo) 
in der (y, z)-Ebene eine Ellipse mit den zeitlich pulsierenden Halbachsen m, h und m, h)y/B, 
und zwar mit einer Umlaufgeschwindigkeit, die nach (37) im ganzen mehr und mehr zunimmt, 
jedoch nach (40) von zeitlichen Schwankungen iiberlagert ist. Man wird sich so nun auch 
wieder den zeitlichen Ablauf des wachsenden Gesamtdrehvektors 0 (w , w\"), wl?) wenig- 
stens qualitativ gut vorstellen kénnen. Die Halbachsen des pulsierenden elliptischen Zylinders, 
innerhalb dessen der Drehvektor 0 im Kreisel sich vergréBert, naihern sich schwankend asymp- 
totisch den festen Werten 


(41) 


a. = M, KY (a3 + 134) R, 


| 
oat: ( (42) 
by = m, k’ V (8% + 1) Ie | 
wobei zur Abkiirzung 
R = [} — € (£0)? + [4 —S (20)? (43) 
gesetzt ist. 
Auch hier kann man diese erste Naherung durch beliebig weit vorgetriebene Naherungen 


héheren Grades verbessern und verfeinern. Die zweite Naiherung_ beniitzt die 
Werte (32) von Oy und @) in der Gleichung (5), so daB man statt (11) 


t 
wo) = Weg +mzt +a f wl) wf) dt (44) 
0 


erhalt und damit statt (30) 


wy) - p wo) wo?) = Ug Mx » | 


(2 (2 


(45) 
ao?) — y wf) of) = ps mz f 
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samt den an die Anfangsbedingungen (31) angepaBten Integralen 


oP = m|n 3 O +a VL 3% 0] 08 BF KO] + 


$ melt Os" (t) — Bs ae Oe (0] sin [VBy K® @], 
~ (46) 
Ge — m|p40|/ Hee 3) — Hs 3 (a eos [¥By K® (@] + 
sm f/E 9000 +m a0 oa 87 x00). 
aa = 2 By K@)(t)! dt, 
[VB y Ko] es 
9° @ = Jsin [By KOO] ae 
und 
y(t) = if wo) dt (48) 
ist. 


Fir die weiteren N&aherungen hat man in diesen Formeln lediglich die oberen 
Zeiger bei wx, Wy, Mz, Yor Ys und K, beginnend mit (44), jeweils um eine Einheit zu erhéhen. 

Man darf in allen Formeln von (30) bis (48) ohne weiteres auch Wy = 0 setzen, wobei dann 
insbesondere in (43) R = 4 wird, und hat so den Fall, daB der Kreisel von der R uh e aus 
durch eine Drehkraft angetrieben wird, deren Momentvektor 9J¢{ mit der x-Achse einen hin- 
reichend kleinen Winkel bildet. 

Ferner macht es keine Schwierigkeiten, die ganze Rechnung auf einen zeitabhan- 
gigen achsenfesten Momentvektor St zu verallgemeinern, also m, als gegebene 
Zeitfunktion anzusetzen, womit dann zufolge (29) auch my und m, sich in gleicher Weise 
zeitlich andern. Man muf dann nimlich nur die bisherige Konstante m, als Zeitfunktion 
m,(t) bei w, behandeln und sie bei w, und w, in die Integrale S® und J mit hereinnehmen 
und hat also als erste Naiherung 


t 
Of On Ve (t) dt 


Uy 
On = 


PB (0) Ey i tus (] cos |/By 7 (oot + f met an)| + 
+ [3 — nef & 30%] sin [YP 7 (o.0t + / moa (49) 
of = — fin 90%) — 15 9*C0] 005 1F 7 (ot + | mete a] + 
+ |o/% 32 *0 +5 900] sin [VF 7 (ont + | / m.(0 a) 
mitiden! Intesealfuuleionen 


gins = [ml ) cos By (0.98 a Jim, (t) i) dt 


| 
U 
= | m, (t) sin Ver (0.08 aE Fi my (t) ay) a | (50) 


Die zweite Naiherung geht aus den Formeln (44) bis (48) hervor, wenn man in 


(44) m,t ersetzt durch fmt dt und in (46) den Faktor m, unterdriickt und ihn dafiir in die 


Integranden der Fan krioaen (47) hereinnimmt, wie das schon in (50) geschah. Entsprechendes 
gilt fiir die weiteren Naherungen. 
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Man erkennt, daB dieses Verfahren sogar auf achsenbewe gliche Selbsterregung 
angewandt werden kann. Man mu8 dann eben auch die nunmehr zeitabhangigen wu, und fs 
in die Integranden der Funktionen (47), (50) usw. hereinnehmen [wo man sie natiirlich mit 
M, ZU Ms M, = my (t) und us m, = m, (t) vereinigt ]. 

Auch hier kann man alle Ergebnisse sofort auf den Fall iibertragen, da® der Moment- 
vektor St der Antriebsdrehkraft in der Nahe der Achse der kleinsten Hauptdrehmasse 
(z-Achse) liegt. Dagegen versagt die Rechnung bei Momentvektoren in der Nahe der mitt - 
leren MHauptachse aus dem gleichen Grunde wie in Ziff. 2. 


4. Die allgemeine iterative Lésung. Zicht man den kraftefreien unsymmetrischen Kreisel 
zum Vergleich heran, so umfassen die iterativen Lésungen von Ziff. 2 und 3 fiir den selbst- 
erregten Kreisel diejenigen Bewegungen, deren Drehvektoren 09 in der Umgebung der 
kleinsten und der gréBten Hauptachse des Tragheitsellipsoids verlaufen, was also den Pol- 
bahnen erster und zweiter Art in der Nachbarschaft jener beiden Achsen beim kraftefreien 
Kreisel entspricht. Wie groB diese beiden Bereiche sind, hangt lediglich von der Konver- 
genz unserer Rechnungen ab, die nicht leicht abzugrenzen ist, sicherlich in der Nahe jener bei- 
den Achsen gut sein wird, aber ebenso gewiB nicht bis zur mittleren Hauptachse reicht. 


Allerdings sind die Lésungen von Ziff. 2 und 3 insofern noch von spezieller Art, als in Ziff. 2 
der Momentvektor 9){ genau in eine jener beiden Hauptachsen fallt, und als in Ziff. 3 der 
anfangliche Drehvektor 0, entweder genau in die zugehérige Hauptachse fallt oder verschwin- 
det. Man kann nun aber — vermutlich innerhalb des gemeinsamen Konvergenzbereiches 
unserer Rechnungen von Ziff. 2 und 3 — den allgemeinen Fall, daB der Momentvektor St und 
der anfangliche Drehvektor 0) beide nicht genau in einer jener beiden Hauptachsen, sondern 
wenigstens innerhalb eines Zylinders liegen, der die betreffende Achse hinreichend eng umgibt, 
offenbar durch eine geeignete Kombination der Lisungen von Ziff. 2 und 3 erhalten. 


Weil die urspriinglichen, exakten Gleichungen (6) und (7) in (30) linearisiert sind, so kann 
man ihrer Liésung (32) die Lésung (13) der verkiirzten Gleichungen (12) additiv hinzufiigen und 
hat dann in ihrer Summe, die wir mit ay) und @) bezeichnen wollen, zusammen mit oJ) (11) 
die erste Naherung der allgemeinen Lésung (in der Umgebung jener beiden Haupt- 
achsen). Fiir die zweite Naherung _ hat man dann in (44) die Werte wy) und 
statt oe und «*) einzufiithren und wiederum zu der mit diesem Wert @®) gebildeten Lisung (46) 
die Lésung (19), ebenfalls mit diesem Wert w) gebildet [also mit @% statt w@ in K()(t) (20)], 
additiv hinzuzufiigen. Ganz entsprechend verfahrt man bei den weiteren Naherungen; wir 
brauchen die Formel!n nicht explizit anzuschreiben, da sie sich nach der jetzt deutlich erklarten 
Vorschrift aus unseren friiheren Formeln von selbst ergeben. ; 

Weil zu vermuten ist, daf der Konvergenzbereich dieser Rechnungen ,,in der Nahe“ der 
kleinsten und der gréBten Hauptachse keineswegs sehr eng ist, so haben wir mit dieser 
Loin binierton Uosune wahrscheinlich einen -erheblichen Teil 
der allgemeinen Lisung des Problems des achsenfest zeit- 
unabhaingig oder zeitabhangig und sogar des achsenbeweglich 
selbsterregten unsymmetrischen Kreisels gefunden. Sie besitzt 
auch dieselbe oo*-fache Mannigfaltigkeit wie die allgemeine Lésung selbst (beziiglich der 
Anfangslage und des Anfangswertes des Drehvektors 0 und beziiglich der Lage und GroéBe des 
Momentvektors Jt der erregenden Drehkraft). 

Wir wollen die Bewegungen, die mit dieser iterativen Liésung darstellbar sind, als as ym p - 
totische Bewegungen erster Art des selbsterregten unsymmetrischen Kreisels 
bezeichnen, weil bei ihnen der Drehvektor 0 mit t—> 0 iiber alle Grenzen wachsend sich 
asymptotisch der x- bzw. z-Achse nahert. Daneben ist, wie die Ergebnisse von Ziff. 5 zeigen 
werden, noch ein zweiter Typ von asymptotischen Bewegungen zu erwarten, bei denen sich 
der Drehvektor asymptotisch einem endlichen Grenzvektor nahert. Diese as ympto- 
tischen Bewegungen zweiter Art kénnen mit der geschilderten iterativen 
Methode nicht dargestellt werden. 


5. Momentvektor in einer Hauptachse; exakte Lésung. Wie schon am Schlufi von Ziff. 2 
angekiindigt, wollen wir nun fiir die drei Falle, daB der Momentvektor St in einer der drei 
Hauptachsen liegt, der Drehvektor 0 aber anfanglich davon abweicht — und zwar jetzt 
beliebig stark davon abweichen darf —, auch noch eine exakte Lésung in 
geschlossener Form angeben, wobei wir uns auf zeitunabhangige Momentvektoren beschranken. 
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Wir beginnen mit dem Fall m, == 0, my = m, = 0 eines Momentvektors in der x-Achse und 
bezeichnen die beliebigen Anfangskomponenten des Drehvektors 0 wieder mit @xo, Wyo, Wzo 
Der Schliissel zur exakten Lésung besteht darin, daf man die Variable 


t 
On Wee dt (51) 
0 


einfiihrt, welche allerdings keine kinematische Bedeutung (etwa als Drehwinkel) hat, weil w, 
nur eine Komponente des Drehvektors 0 ist, und weil dieser zudem im Kérper wandern kann. 
Mit (51) und my = m, = 0 haben die beiden Gleichungen (6) und (7) die allgemeinen Integrale 


ee EL | 
oz = 0, Pysin (YFP G« +9). J 


wobei w, und o die fritheren, durch (14) und (15) definierten Abkiirzungen sind und voraus- 
setzungsgemaB auch hier @,=+ 0 ist; die Quadratwurzeln in (52) darf man wieder als positiv 


festsetzen. 
Um ¢, als Funktion der Zeit und damit nach (51) auch vollends w, zu finden, fiihren wir die 
Integrale (52) in die Gleichung (5) ein und erhalten 


Wy =m, +ho%a VBP y sin 2 (VB y ex +24) (53) 
und hieraus durch Multiplikation mit w, = dy,/dt und Integration 


oy = VU) | (54) 
mit der Funktion 


Uo] 0% 208; Gy Oe eos 2 (VB yo2 +o)—cos2o]. (55) 
Weiter folgt aus (54) mit (51) durch Integration 


(52) 


Vx 
dfx 
7 ee (56) 
| /-U (gx) 
Die Umkehrfunktion sei 
Px = Px (t) » (57) 


und hiermit ist die Lésung (52) und (54) als Funktion der Zeit gefunden. Uber das Vorzeichen 


von JU in (54) und (56) wird nachher noch eine Festsetzung zu treffen sein. 

Das Integral (56) ist offensichtlich die Erweiterung des elliptischen Integrales erster Gat- 
tung, auf welches man bei der Behandlung des kraftefreien unsymmetrischen Kreisels stéBt, 
und in das es mit m, = 0 iitbergeht; die Umkehrung (57) entspricht der elliptischen Amplituden- 
funktion, und die Liésungen (52) und (54) entsprechen den Jacobischen Funktionen cn, sn und 
dn. Weil im Gegensatz zu diesen Funktionen die Funktion ¢, (t) bisher nicht tabuliert ist, so 
wird die zahlenmaBige Auswertung unserer jetzigen exakten Lisung vielleicht miihsamer sein 
als die frithere iterative Lésung von Ziff. 2, wenigstens soweit jene gut konvergiert (also fiir 
kleine Werte von w,). Dafiir wird uns aber die exakte Lésung weitere Bewegungstypen auf- 
zeigen, die in dem beschrankten Geltungsbereich (Konvergenzbereich) der iterativen Lésung 
noch nicht enthalten waren. 

Man mu8 im wesentlichen zwei Bewegungstypen sowie einen Grenzfall unterscheiden und 
erhalt diese, wenn man die beiden Bestandteile der Funktion U (55) iither einer Abszisse gp, 
als Kurven aufzeichnet, namlich die Gerade wi, +2m,qy, und die Cosinuslinie c (px) 


=y wi a [cos 2 (VB Y Px + o) == (oS 2 ol. Hierbei kénnen wir ohne Kinschréankung der All- 
gemeinheit Wx als positiv voraussetzen, miissen dann aber fiir m, positive und negative Werte 
zulassen (diese Vorzeichenfestsetzung im Gegensatz zu Ziff. 2 hat hier gewisse Vorteile fiir die 
Darstellung). Je nachdem, zunachst fiir m, > 0, die Gerade die Cosinuslinie bei positiven 
Abszissen y, nicht schneidet oder doch schneidet, erhalt man die beiden Typen von Abb. 1 
oder von Abb. 2. Dort ist jeweils der Bereich der positiven U-Werte schraffiert, welche wegen 
(54) und (56) fiir die Bewegung allein in Betracht kommen, und die Schraffurstriche stellen 
die Werte qm? dar. 

Im ersten Fall (Abb. 1) nimmt die Drehkomponente w, im ganzen, wenn auch von mehr 
oder weniger starken Schwankungen tiberlagert, mit der Zeit iiber alle Grenzen zu, wie es im 
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unteren Teile von Abb. 1 skizziert ist, und dies ist offenbar der schon durch die iterative Lé- 
sung gefundene Typ derasymptotischen Bewegungen erster Art (Ziff. 2). 
Die Quadratwurzel yu ist dabei in (54) und (56) positiv zu nehmen, da @,) voraussetzungs- 
gema8 positiv ist und dann ¢y, in der Tat mit t wachst. Aber man kann natiirlich jU und 
also m, und damit nach (56) auch dy, negativ wahlen, und dies bedeutet, da® nun die Gerade 
Wx + 2 m, py von beliebig groBen Werten ~, aus nach links hin durchlaufen wird: die positive 
Selbsterregung m, zieht einen Drehvektor, der anfanglich in der Nihe der negativen x-Achse 
lag und einen beliebig groBen Betrag hatte, in die Lage Wxy > Wyg, Wz Zuriick. 

Als Beispiel fiir den zweiten Fall nehmen wir an (Abb. 2), daB8 wi, kleiner als der Héchst- 
wert von c(px), und da8B m, schwach positiv sei. Dann nimmt @,, mit Wx beginnend, nach 
Uberschreitung eines Hichstwertes bis zu Null ab (fiir y, =, und t =t,), und bis dahin war 
yU positiv zu nehmen. Da beim Weiterschreiten auf der Geraden U negativ, also JU und 


damit auch wm, imaginiar wiirden, so ist jetzt das Vorzeichen von jU umzudrehen: p, nimmt 
nach (56) ab, w, wird negativ, geht iiber einen negativen Héchstwert zu Null (bei 9, =Qxo 
und t =t,), kehrt dann nach abermaligem Vorzeichenwechsel von }'U zum Anfangswert @zo 
(fiir ¢ = T) zuriick, und nun beginnt das Spiel von neuem. Da nach (52) auch die Dreh- 
komponenten @, und @, entsprechend hin und her schwanken (allerdings nur mit Qy2 < x 
<x, also nicht unbedingt zwischen positiven und negativen Werten), so ist die Bewegung 
jetzt hinsichtlich der Wanderung des Drehvektors 0 im Kérper periodisch, und zwar mit 


! 4 ih 
Abb. 1. Asymptotischer Bewegungstyp fiir 
Ps Oat, on — 0. Vr, 
ees a Wi, 
ie Pat 
einer Periode T, die durch das iiber die ge- 
schlossene @,-Kurve von Abb.2 unten er- 
streckte Integral 
Abb. 2. Periodischer Bewegungstyp fir 
Appx (58) m= = 0, my =m,= 0. 


VU(ge) 
dargestellt ist. Diese Periodizitat des Drehvektors 0 im Kérper bedeutet aber keineswegs auch 
eine villige Periodizitat der Bewegung des Kreisels im Raume (wie man ja schon von der 
Poinsothewegung des kraftefreien Kreisels wei). Trotzdem wollen wir diesen Bewegungstyp 
periodisch nennen. 


Die Quadratwurzel J Uist in (58) positiv oder negativ zu wahlen, je nachdem jeweils dy, 
positiv oder negativ ist. 

Man erkennt, da®B Abb.1u. 2 und mit ihnen der Bewegungsverlauf mannigfacher Abwand- 
lung fahig ist. Wenn etwa in Abb.2 die MomentgréBe m, stark negativ gewahlt wird, so bleibt 
zWar (x1 reell, wogegen die Gerade wyo + 2 m, y, fiir negative p, keinen reellen Schnittpunkt 
xo mit der Kurve c(y,) besitzt. DemgemaS nimmt in diesem Fall die urspriinglich positive 
Drehkomponente w, nach ihrem Durchgang durch Null negative Werte an, deren Betrag dann 
fiir fortschreitende negative Werte von ¢,) wie im Falle von Abb. 1 mit mehr oder weniger 
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starken Schwankungen iiber alle Grenzen wiachst. Ist wi) groBer als der Héchstwert von c(qz), 
so kann es fiir nicht allzu stark negative m, vorkommen, daB die Gerade Wry |+ 2m, x die 
Kurve c(gy) erst rechts von ihrem ersten Héchstwert schneidet (in Abb. 2 strichpunktiert). 
Man wird sich leicht klarmachen, wie alsdann der Wert von w, ttber eine (oder auch mehrere) 
Schwankungen zu Null abnimmt und darnach schwankend zu immer starker negativen Werten 
aufsteigt. 

Offenbar sind insbesondere die Bewegungstypen von Abb. 1 und 2 die Verallgemeinerungen 
der Poinsotbewegung des kraftefreien unsymmetrischen Kreisels. Wenn Wx. groB genug ist — 
dies entspricht den Polbahnen erster Art der Poinsotbewegung —, so hebt jede von Null ver- 
schiedene (positive oder negative) MomentgréBe m, die Periodizitat der Bewegung auf und ver- 
wandelt diese in eine asymptotische mit |w,| > co. Wenn dagegen der (positive) Wert @yp hin- 
reichend klein ist — dies entspricht den Polbahnen zweiter Art —, so wird durch eine Moment- 
gréBe m, von hinreichend kleinem Betrag die Periodizitat der Bewegung nicht aufgehoben. 
Man kann mit (52) und (53) beweisen, daB der erste oder zweite Fall vorliegt, je nachdem 


O> &  C(B—C) 
oy 4(A—B8) Wee) 


ist, je nachdem also der anfangliche Drehvektor 0) in demjenigen Keil zwischen den beiden 
Ebenen der sogenannten trennenden Polbahn! (des kraftefreien Kreisels) liegt, welchem 
die x- oder die z-Achse angehért. In den Darstellungen gemaf Abb. 1 und 2 unterscheiden 
sich die beiden Falle in der Form w%o 2 ¢(Px)max3 Vgl. die ausgezogene und die strichpunktierte 
Gerade von Abb. 2.) Kurz gesagt: umtanzt der Drehvektor 0 des zunachst kraftefreien 
Kreisels die x-Achse, so 1aBt je des Selbsterregungsmoment M, den Drehvektor mit der Zeit 
zu unbeschrankt hohen Betragen anwachsen, wobei aber seine Komponenten m, und w, be- 
schrankt bleiben, so daB der Drehvektor nach (52) wenigstens innerhalb eines elliptischen 
Zylinders mit der x-Achse als Achse und den Halbachsen @, VB und wo, Vy bleibt, wie dies in 
Ziff. 2 geschildert war; umtanzt der Drehvektor dagegen urspriinglich die 2-Achse, so bleibt 
bei hinreichend schwacher Selbsterregung M, die Periodizitat des Drehvektors im 
Kreisel erhalten. Es ist daher verstandlich, da die iterative Methode von Ziff. 2 nur den 
asymptotischen Bewegungstyp liefern konnte: nach den dortigen Voraussetzungen (10) durfte 
die Selbsterregung m, nicht ,,schwach“ sein. 

Falls o nicht, wie in Abb. 1 und 2, positiv, sondern 
negativ ist, wird man sich die entsprechenden Figuren 
mit ihren beiden Bewegungstypen leicht klar machen 
koénnen, ohne da wir sie anzeichnen. 

Nun bleibt noch der Grenzfall, daB die Gerade 
wx? +2 m, pz die Cosinuslinie im Punkt ¢,, berihrt; 
er ist in Abb. 3 dargestellt fiir den Fall, daB die Be- 
rithrung ,,von oben“ erfolgt. Die Zeit t wird jetzt 
fiir ~, = Px, logarithmisch unendlich, und die Be- 
wegung nahert sich nun asymptotisch dem Dreh- 


vektor 0, mit der Komponente w,, = 0, also einer 

@. 7 7 
a Drehung um eine feste Achse in der (y, z)-Ebene, 
und dies ist der Typ der schon in Ziff. 4 erwahnten 
Oxo . asymptotischen Bewegungen zweiter 

Pe Art. Weil dort dp,/dt = 0, also 

Gat ° dwx divx 
Abb. 3. Grenzfall der stationiren Bewegung. Ee dyx dp =0 


wird, so ist diese Drehung dort pe r manent geworden, und somit ist nach (5) 
Mz, = — & Wy O, 
oder wegen (3) und (4) 
M, = — (B— C)o,o,, (60) 
und dies bedeutet, daB M, gerade das Moment ist, welches diese permanente Drehung aufrecht- 
erhalten kann?. Weil die Drehung (60) aber nicht nur durch das Moment M, asymptotisch er- 


* R. Grammel, Der Kreisel, Bd. 1, S. 131, Formel (16). 
2 R. Grammel, Ing.-Arch. 21 (1953), S, 150. 
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reichbar ist, sondern auch (vgl. die weitere Schraffur in Abb. 3) bei der geringsten Stérung 
wieder verlassen werden kann, so ist diese stationdre Bewegung alsinst abil zu bezeichnen, 
und da der Drehvektor dann (schwankend) iiber alle Grenzen wachst, so kann also eine asymp- 
totische Bewegung zweiter Art in eine solche erster Art iibergehen, und umgekehrt. 

Wenn hingegen die Berithrung der Geraden Ox, + 2m, y, an der Cosinuslinie ,,von unten“ 
erfolgt, wie dies in Abb. 3 durch die gestrichelte Linie angedeutet ist, so hat man es offenbar 
dort mit einer st a bile n permanenten Drehung zu tun. 

Man kann zeigen, daB dies die gleichen Bereiche der Instabilitat und Stabilitat sind, die 
von den permanenten Drehungen um Achsen in der (y, z)-Ebene bekannt sind!. Sie werden 
getrennt durch die Wendepunkte der Cosinuslinie; und von denjenigen permanenten Dre- 
hungen, die durch diese Wendepunkte dargestellt werden, geht nur eine asymptotische Be- 
wegung erster Art aus, aber keine solche zweiter Art. 

Es versteht sich von selbst, da man in der gleichen Weise die exakte Lésung fiir den Fall 
m, = my = 0, m, + 0 eines Momentvektors St in der z-Achse findet ; man braucht nur wieder 
die maBgebenden GréBen entsprechend zu vertauschen. Auch hier gibt es die soeben diskutier- 
ten drei Bewegungstypen samt ihrem stationaren Grenzfall. 

Aber auch die exakte Liésung fiir den Fall eines Momentes St in der mit tl eren Haupt- 
achse (y-Achse) kann man auf diesem Wege finden. Mit m, = m, = 0, my =- 0 folgt aus (5) und 
(7) entweder 


Wy = 0 Va Sof (Vaygy +o), | 


a 1y (61a) 
| @eia= (os Vy Sin (Vo y Py ae °) j 
mit ae w3, 2 way Fao == o (62a) 
aaa. Wes bau tee 
oder = Fae 
Qs =O), Vo Sin (Jo Y Py + 0) 2 | (61b) 
| O55 Vy Sof (Vou y Py 7 0) J 
oe We Ox oe oe 
Oe = G0 =) — ae 
» i a Wz 
und dann aus (6) 
Oy = VV @y) oe 
mit t 
Py al @y dt me 


und der Funktion 


V(py) = yy + 2 my Py — § OF BLGoj 2 (Vay gy + 0) — Gof 2a). (65) 
Von den beiden Lésungen (61a) und (61b) hat man jeweils diejenige zu nehmen, fiir welche @, 
und o nach (62a) und (62b) reell werden, und das heift: die erste oder die zweite, je nachdem 
Ox & _ C(B—C) (66) 

of, ~ 7% A(A—B) 
ist, je nachdem also der anfangliche Drehvektor 0, in demjenigen Keil zwischen den beiden 
Ebenen der sogenannten trennenden Polbahn liegt, welchem die x- oder die z-Achse angehért. 
Der Fall w2o/w2) =a/y, daB der anfangliche Drehvektor 0, in einer der beiden Ebenen der 


trennenden Polbahn liegt, ist hier auszuschlieBen, da fiir ihn mit @, —0 und o =+ © die 
Lésungen (61a) und (61b) unbrauchbar werden; wir werden ihn in Ziff. 6 behandeln. 
Weiter ist 
ie 
nee | peae is (67) 
/ Vey) 


und umgekehrt 
Py = Pylt) - (68) 
1 R. Grammel, Ing.-Arch,. 21 (1953), S. 152. 
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Die Quadratwurzeln in (61) und (62) darf man wieder positiv festsetzen, diejenigen in (63) und 
(67) richten sich nach dem Vorzeichen von dy. 

Um den Bewegungstyp zu erkennen, trigt man die beiden Bestandteile der Funktion V (65), 
von denen wir den zweiten kurz mit C(py) bezeichnen wollen, iiber der Abszisse yy auf und er- 
halt so Abb. 4. Man erkennt, daf jetzt die Bewegung durchweg dem Typ von Abb. 2 angehort 


und periodisch ist, und zwar mit der Periodendauer 
ie f Be TEES (69) 
/Vigy) 


Die Pfeilspitze des Drehvektors p lauft, wie aus (61a) oder (61b) folgt, auf einem hyperbolischen 
Zylinder, der die y-Achse zur Achse hat und entweder der Gleichung 


ae Ore (70a) 
ox wiry 
oder der Gleichung 
2 2 
=e ss (70b) 


GIRO OG 


gehorcht und also das gleiche Achsenverhaltnis hat wie die 
hyperbolischen Zylinder gleicher Achse, welche beim krafte- 


T freien Kreisel die Polbahnen erster oder zweiter Art aus dem 
ete Poinsotellipsoid ausschneiden. Im Gegensatz zu jenen Pol- 


bahnen haben aber die jetzigen Polbahnen infolge des Mo- 
mentes Jt ihre Symmetrie beziiglich der (x, y)- und (y, 2)- 
Ebene verloren. 

Wiahrend ein Momentvektor )J¢ in der x-Achse oder in 
der z-Achse, wenn er nur gro genug war, die vom krdafte- 
Py freien unsymmetrischen Kreisel her bekannte Periodizitat 
des Drehvektors iitberwinden und den Drehvektor p iiber alle 
WV Grenzen anwachsen lassen konnte, so ist dies bei einem Mo- 
fy mentvektor in der y-Achse, auch wenn er noch so gro wird, 


fe ie Pyr nicht mehr der Fall (mit einer singuladren, in Ziff. 6 fest- 
az zustellenden Ausnahme), — ein unerwartetes und sehr para- 
dox erscheinendes Ergebnis. Allerdings erkennt man aus 


Abb. 4. periodischer Bewegungstyp fir Abb. 4, da ein Momentvektor { von ungemein hohem Be- 

my= 0 Mg = m, = 0. trag (also sehr steilem Anstieg der dortigen Geraden) den 

Schnittpunkt ~,, oder gy, beliebig weit vom Nullpunkt weg- 

riicken kann, so da auch hier beliebig hohe Betrage von w, erreichbar sind; aber schlieBlich 

muffs, wenn auch nach langer Zeit, w, doch wieder abnehmen. [Wir werden in Ziff. 6 sehen, 

da8 erst, wenn die Zylinder (70) in die beiden Ebenen der trennenden Polbahn ausarten, die 
Periode T (69) unendlich gro8& werden kann. ] 

Es sei der Vollstandigkeit halber noch erwahnt, da auch hier der Grenzfall einer per ma - 
nenten Drehung vorkommt: wenn mit wy)= 0 die Gerade w},-+ 2 my py in Abb.4 die Coj- 
Kurve im Nullpunkt beriihrt, so schrumpft namlich das schraffierte Gebiet in diesen Punkt zu- 
sammen. Diese stationare Drehung mit w, = 0 um eine Achse in der (x, z)-Ebene ist sta bil }, 
da kein schraffierbarer Nachbarbereich existiert, und allgemein stellt nicht nur der Nullpunkt, 
sondern jeder Punkt der Gpj-Kurve eine der bekannten stabilen permanenten Drehungen um 
alle Nullpunktsachsen der (x, z)-Ebene dar. 

Wie schon festgestellt, kann unsere Lisung nicht bis zur trennenden Polbahn des krafte- 
freien unsymmetrischen Kreisels vordringen. Auch der Konvergenzbereich der Lésungen von 
Ziff. 2 bis 4 kann dies nicht. Es ist deshalb erwiinscht, jetzt auch noch diejenigen Bewegungen 
des achsenfest selbsterregten unsymmetrischen Kreisels zu finden, bei denen der Drehvektor p in 
einer der beiden kérperfesten Ebenen lauft, welche aus dem Poinsotellipsoid die trennende Pol- 
bahn ausschneiden. Da dieser Fall beim kraftefreien Kreisel mit elementaren Funktionen ge- 
lést werden kann?, so ist zu vermuten, daB er auch beim selbsterregten Kreisel, wenigstens 
bei zeitunabhangig achsenfester Selbsterregung, eine verhaltnismaBig einfache Lésung zulaBt. 


* Genauer gesagt: parastabil, vgl. R. Grammel, Ing.-Arch. 21 (1953), S. 152. 
2 Siehe R. Grammel, Der Kreisel Bd. 1, S. 147. 
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6. Der Fall der trennenden Polbahn. Man erhiilt die Lésung fiir diesen Fall, wenn man sich 
die Frage vorlegt, ob ein kiérperfester zeitunabhaingiger Momentvektor QM in 
der mittleren Hauptachse den Drehvektor 9 in einer kérperfesten Ebene bewegen 
kann, welche durch die mittlere Hauptachse hindurchgeht, das heiBt: ob die Gleichungen (5) 
bis (7) mit 

m, = 0, my = konst. =0, m, = 0 (71) 
durch einen Ansatz von der Form 


Ox = Wy f(t), wy =Wg(t), @, =A, f(t) (72) 
befriedigt werden kénnen, wo q@y ein noch offener positiver Faktor von der Dimension einer 
Drehgeschwindigkeit ist und f(t) und g(t) noch zu bestimmende Funktionen der dimensions- 
losen Zeitvariabeln 

enya (73) 
sind und 4 ein ebenfalls noch unbekannter Festwert sein soll, der offenbar den Tangens des 
Keilwinkels y zwischen jener kérperfesten Ebene und der (x, y)-Ebene angibt. 

Die Gleichungen (5) bis (7) gehen mit (71) bis (73) iiber in 


f’—adfg=0, (74) 
s+ pif? ==, (75) 
Ta ee (76) 


wobei Striche fortan Ableitungen nach 1 bedeuten sollen. Damit (74) und (76) miteinander 
vertraglich sind, muB «A = y//, also 


pet ey) ete 8) 
A +2 =+)eG=G 2) 
sein, und das besagt, daB in jedem unsymmetrischen Kreisel nur zwei solche Ebenen existieren, 
symmetrisch zur (x, y)-Ebene, und in der Tat gerade die beiden Ebenen der sogenannten tren- 


nenden Polbahn des kraftefreien Kreisels sind. 
Dann folgt aus (76), unter Ausschaltung des trivialen Falles f = 0, 


g() 2 F, (78) 


und mit diesem Wert von g geht (75) iiber in 


ieee get year 5 (79) 


wobei zur Abkiirzung die dimensionslosen GréBen 


f= poe 0’, r=aA—% (80) 


eingefiihrt sind und auf (77) geachtet ist. 

Weil die nichtlineare Differentialgleichung zweiter Ordnung (79) die Zeitvariable 7 nicht 
explizit enthalt, kann man sie in eine Differentialgleichung erster Ordnung verwandeln, und 
zwar durch die Substitution 


FD a) uy) (81) 


in eine solche mit separablen Veranderlichen, namlich 
du 2 
ae — = 82 
He i +2 qf 7s (82) 


und diese hat das allgemeine Integral 
u=D—gqf?+2rlnf. (83) 


Die Integrationskonstante D bestimmen wir, indem wir, offenbar ohne wesentliche Kin- 
schrankung der Allgemeinheit, verabreden, daB die Bewegung fiir t = + = 0 mit einem Dreh- 
vektor beginnen soll, der in der (x, z)-Ebene liegt, also senkrecht zum Momentvektor M steht, 
und dessen Komponente w, die GréBe wy haben soll. Wir setzen also 


Cais) fa 0). stir r= 0, (84) 
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wobei die letzte dieser drei Bedingungen wegen der vorletzten aus (74) oder (76) folgt. Dies gibt 
dann zufolge (81) auch u = 0 fiir 7 = 0 und daher nach (83) D = q. Mit diesem Wert folgt aus 
(81) und (83) durch eine weitere Integration 

Ap 


Ean if teary (85) 
é lA" 


1 


wobei tiber das Vorzeichen der Quadratwurzel nachher noch eine Vorschrift folgen wird. Wir 
denken uns diese Funktion umgekehrt; dann ist mit 


f=f) (86) 
die Zeitfunktion f gefunden. Aus (78) folgt mit (81) und (83) sofort vollends 


g(c) =T qd —f) F2rinf, (87) 


wiederum vorbehaltlich einer noch zu treffenden Vorschrift tiber das Vorzeichen der Quadrat- 
wurzel. 

Schon aus der Differentialgleichung (79) ist, wegen der verschiedenen Vorzeichen von r (80) 
in den beiden Ebenen (77), ersichtlich, daB man auch bei gleichem Wert my und bei gleichen 
Werten a, B, Y, Wp, also gleichen Werten q und |r| verschiedene Verlaufe des Drehvektors 0 in 
jenen beiden Ebenen zu erwarten hat. Dies sieht zunachst paradox aus, da doch die beiden 
Ebenen vollig gleichwertig zu sein scheinen, wird aber verstandlich, wenn man beachtet, daB 
im einen Falle (A > 0) die positive x-Achse, der (fiir my > 0) in die positive y-Achse fallende 
Momentvektor 3Jt und der anfangliche Drehvektor (in der (x, z)-Ebene) fiir @) > 0 ein rechts- 
handiges (teilweise schiefwinkliges) Achsenkreuz bilden, im anderen Falle (A < 0) dagegen ein 

linkshandiges. Wir werden jetzt durch eine genauere 
Untersuchung der Liésung bestatigen, daB in der Tat 
die Bewegung in beiden Fallen (nicht nur quantitativ, 
he sondern auch qualitativ) vollstandig verschieden ver- 
lauft, und wir werden sehen, daB man im Falle r > 0 
Inf auch noch die drei Unterfalle r= q mit verschiedenem 
Geprage der Bewegung unterscheiden muf. 


u 
i 3 
J LLL | if 4 ae 
2 3 1 é ; 
Abb. 5. Die Funktionen 1 — f? und Inf. Abb. 6. Uherlagerung der Funktionen 1 —f? und Inf. 


Wir gehen aus von der Funktion u(f) (83) mit D = q: 
u =q(1—f?) +2 rInf (88) 


und betrachten zunachst den Fall r > 0, also auch 2 > 0, da man ohne Einschrénkung der 
Allgemeinheit weiterhin my als positiv voraussetzen darf. Tragt man die beiden Funktionen 
1 — f? und In f itber der Abszisse f als Kurven auf (Abb. 5), so erkennt man, daB ihre Uber- 
lagerung mittels der positiven Faktoren q und r je nach der GréBe von q und r drei verschiedene 
Kurventypen fiir u(f) ergeben kann, namlich fiir r< q die Type 1 von Abb. 6, fiir r = q die 
Type 2 und fiir r > q die Type 3. Man erhalt namlich aus (88) 


du . 
(#),_ = 20-90 fir rSq 


und 
(a). —? (qQtrn<0. 


Die Kurventype 1 schneidet die f-Achse auBer in f = 1 noch in einem zweiten Punkt mit der 
Abszisse f,, wobei 0< f,< list. Ebenso hat die Kurventype 3 auBer f = 1 noch einen Ab- 
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szissenschnittpunkt f, > 1. Beide Typen besitzen positive u- Werte in dem Bereich zwischen 1 
und f;, und nur in diesem Bereich. Weil nun nach (85) und (87) 


ot (emus e(r) =-yulf) 89 
= | Tach a(t) = Yulf) > 


ist, so kommen fiir eine reelle Bewegung nur diese (in Abb. 6 wieder schraffierten) Bereiche in 
Frage. 

Beginnt die Bewegung, wie wir es verabredet hatten, mit f = 1 zur Zeit t = 0, so wandert 
u auf den Kurven 1 (fiir r<_ q) bzw. 3 (fiir r > q) vom Punkt f = 1 gegen den Punkt f, hin und, 
sobald dieser erreicht ist, wieder zu f = 1 zuriick, und dann beginnt das Spiel von neuem. Dem- 
gema hat man fiir 7 (89) und die erste Halbperiode der Bewegung folgende Vorschrift iiber das 
Vorzeichen von yu zu treffen: 


negativ| .... 

eee fare rg und) PSs fp 
Die gleiche Vorzeichenregel gilt natiirlich auch fiir 
fi=fyu und also wegen (78) auch fiir g (r)(89). Ge- 


rade die umgekehrte Vorzeichenvorschrift gilt in der 
zweiten Halbperiode: 


Ss 


—{positiv | |. 
yu “) negativ | fur Sq und f—>1. 


Die Periodendauer ist im ganzen 


fi 


ay 
T=2 90 
| fu(f) ie) 


Abb. 7. Die beiden Polebenen I (A > 0) und II (4 < 0). 


Wir wollen nun den Verlauf des Drehvektors 0p fiir 1 > 0, also r > 0 in den beiden Fallen 
r S q, also die zugehérige Polbahn figiirlich darstellen. Abb. 7 zeigt zunadchst schematisch das 
Hauptachsenkreuz samt den beiden Ebenen (77) und je einer f-Achse in diesen Ebenen: I fiir 
A > 0Ound II fiir, < 0. In Abb. 8 sind diese Ebenen, genauer gesagt, die ,,oberen‘* Halbebenen 


-3 


Ss 
S, 


1 My 


Abb. 8. Die Polbahnen in den beiden Polebenen J und J7, fiir my 20,m,=m,=0. 
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in die (y, 2)-Ebene hineingeklappt und dann zunachst in der rechten Halbebene I die (f, g)- 
Kurven aufgezeichnet, und zwar durch Auswertung der Gleichung (87) je fiir die Parameter- 
werte r=} q, 44,2 4,4 q, und versehen mit den Pfeilen der Richtung, in der sie durchlaufen 
werden, wenn my > 0 ist. Dabei sind folgende Werte zugrunde gelegt: A: B: C=4 3332, 
also A =Y2,q=4,A/y =2 /2. In Wirklichkeit miBte man die Abszissen f noch mit dem 


Faktor w) J1 +/?, die Ordinaten g noch mit dem Faktor a multiplizieren; da dies aber nur 
eine Sache des Mafstabs ist, so kann man die gezeichneten Kurven schon als die Polbahnen 
(Kurven der Pfeilspitze des Vektors 0) ansehen. (In Abb. 8 ist der Ma®stabsfaktor V1 +2? bei 
f beriicksichtigt; die Kurven sind also die Polbahnen fiir @) = 1.) 

Man stellt zundchst fest, daB® der selbsterregte Kreise] mit einem Erregungsmoment in der 
y-Achse bei einem anfanglichen Drehvektor 0) in der Ebene I nur Bewegungen vom pe - 
riodischen Typ besitzt. 

Der Vollstandigkeit halber ist in Abb. 8 auch noch fiir r = 0, also my = 0 die bekannte 
trennende Polbahn des kraftefreien Kreisels als Viertelsellipse 


f? ste pao enal (91) 


eingetragen, wobei nach (85) mit f = 0 (und r = 0) $T = oo wird, und man stellt anhand 
des bekannten Poinsotschen Bildes der Bewegung leicht fest, daB der eingetragene Pfeil dieser 
Kurve richtig ist. [Wir wollen hier nur noch erwahnen, daB die Lésung (86), (87) von (72) fiir 
r = 0 in die bekannten Formeln der Bewegung im Falle der trennenden Polbahn des krafte- 
freien Kreisels ibergehen muB8, und in der Tat bestatigt man dies durch eine einfache Rechnung. | 

Aus Abb. 8 liest man ab, daB fiir r< q (schwaches Moment Jt) der Drehvektor 0, der zu 
Beginn der Bewegung senkrecht zum Momentvektor (y-Achse) steht (09), sich zunachst 
entgegengesetzt zur Richtung des Momentvektors Jt weiterbewegt, dann in eine 
Richtung senkrecht zu Yt zuriickkehrt (0;) und dabei kleiner als Dy geworden ist, und daB er 
nunmehr im Sinne des Moments St weiterwandert und nach der Periode T (90) wieder 
in seine Anfangslage (0,) zuriickgekommen ist. Fiir r > q (starkes Moment) bewegt er sich von 
Anfang an zunadachst im Sinne des Momentvektors Jit, spaterim entgegengesetz- 
ten Sinne iiber eine Lage (oj) senkrecht zu 9J{ und dann wieder zuriick zu seiner Anfangs- 
lage (0p). 

Man bemerkt aber, daB ein grundsatzlicher Unterschied zwischen diesen beiden Bewegungs- 
typen eigentlich nicht besteht. Wiirde im Falle r< q die Bewegung mit dem Anfangsvektor 
0; beginnen, so ware sie vom zweiten Typ, und das ist einfach so zu erklaren, daB dieser An- 
fangslage ein kleinerer Wert @,) entsprache (wenn dort f = | sein soll), um so viel kleiner, daB 
dann r (80) gréBer als q ware; und Entsprechendes gilt fir den zweiten Typ, der mit pj als 
Anfangsvektor sich in den ersten Typ verwandelt. Es ist eben nicht die Gré8e von my, allein 
fiir die Fallunterscheidung rS q mafgebend, sondern der Quotient m,/wé und also der Quotient 
|Me|/o§, und man kann das Verhalten des Kreisels etwa dahin kennzeichnen, daf die fiir w) > 0 
schon beim kraftefreien Kreisel vorherrschende Bewegungstendenz des Pols auf der Polbahn 
entgegengerichtet zur positiven y-Achse durch einen Momentvektor St in der positiven y-Achse 
nur mit einem hinreichend grofen Verhaltnis Me| /we itberwunden werden kann, und darin liegt 
auch die sachliche Erklarung fiir dieses eigenartige und zundchst wohl paradox erscheinende 
Verhalten der Polbahn in der Ebene J. Das Grenzmoment My, das beide Falle trennt, hat 
zufolge (80) mit (3), (4) und (77) den Wert 

A(A—B ‘ 
Mea Vee? ag. (92) 
Ist M, kleiner als dieser Wert, so hat man den Kurventyp r < q von Abb. 8, ist es gréBer, 
so hat man den Kurventyp r > q. 

Jetzt bleibt (fitr A > 0) noch der Grenzfall r = q. Hier ist nach Abb. 6 nur der Wert f = 1 
méglich. Dazu gehort wegen f’ = 0 nach (74) oder (76) g = 0 und also nach (75) my, = Blo 
oder zufolge (4), (3) und (72) 

M, = (A— C)o, a, (@, —=10) 7, (93) 


und dies ist wieder gerade das Moment, welches eine solche permanente Drehung unterhalten 
kann +. [Natiirlich ist (93) nur ein anderer Ausdruck fiir das Moment (92).] 


1 R. Grammel, Ing.-Arch, 21 (1953), S. 150. 
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Der Fall r< 0, alsod < 0 der Ebene IT erledigt sich vollends rasch. Hier ist nach (88) die 
Funktion u(f) positiv fiir 0 < f <1 und negativ fiir f > 1; in 7 (89) und g (r) (89) ist nun der 
Quadratwurzel Yu das negative Vorzeichen zu geben, und man erhilt eine nichtperiodische 
Polbahn, wie sie in Abb. 8 links in der ,,vorderen‘* HalbebeneII fiir die friiheren Verhialtnisse 
A:B:C und fiir r = — q aufgezeichnet ist: Der Drehvektor 0 wird von dem Momentvektor 
Wt asymptotisch in die positive y-Achse hineingezogen. Dies riihrt davon her, da® schon beim 
kraftefreien Kreisel (r = 0; vg]. die Viertelsellipse in Abb. 8 links) seine Bewegungstendenz fiir 
@y > 0 zur positiven y-Achse hin gerichtet ist. Somit gibt es in der Ebene II nur Bewegungen 
vomasymptotischen Typ erster Art und dies ist zugleich der in Ziff. 5 an- 
gekiindigte einzige Ausnahmefall, daB ein Selbsterregungsmoment in der y-Achse den Dreh- 
vektor iiber alle Grenzen anwachsen lassen kann. 

Wenn man auch den zeitlichen Lauf des Poles auf den Polbahnen in Abb. 8 kennen, die 
Polbahnen also mit Zeitmarken versehen will, so mu8 man natiirlich auch noch das Integral 
(85) auswerten und umkehren, was wir aber hier nicht auszufiihren beabsichtigen. 

Wohl aber sei noch erwahnt, daB die Polbahnen in den ,,unteren‘‘ Halbebenen I und II 
von Abb.7 einfach dadurch erhalten werden, da®B man die Kurven in Abb. 8 samt ihren Pfeilen, 
je getrennt fiir die beiden Falle 1 > 0 und/< 0, an der g-Achse spiegelt. Dann wenn man das 
Gebilde von Abb. 7 um die y-Achse so dreht, daB® ,,unten‘* und ,,oben‘‘ vertauscht werden, 
so andert sich sachlich nichts. 

Da die einzigen stabilen permanenten Drehungen durch Selbsterregungsmomente 
unterhalten werden kénnen, deren Vektoren 9){ in einer der drei Hauptachsen liegen?, und 
da die Scharen der periodischen Bewegungen wohl stets als Nachbarbewegungen von stabilen 
permanenten Drehungen angesehen werden diirfen, so ist mit groBer Sicherheit zu vermuten, 
daB die in Ziff. 5 und 6 gefundenen periodischen Bewegungen die ein zi gen periodischen 
Bewegungen sind, deren der selbsterregte unsymmetrische Kreisel fahig ist. 


7. Erweiterung des Falles der trennenden Polbahn. Man kann die vorangehende Lésung 
noch erheblich erweitern, wenn man zulaBt, daB der Momentvektor Jt nicht mehr in die mitt- 
lere Hauptachse fallt, so daB auBer my auch m, und m, von Null verschieden sein kénnen, jedoch 
immer noch verlangt, da®B der Drehvektor 0 in einer der beiden Ebenen verbleibt, die durch 
(77) gekennzeichnet sind. Jetzt hat man mit (72) statt (74) bis (76) 


f—aafe=z- (94) 
Scat ae (95) 
Be as (96) 


Damit (94) und (96) miteinander vertraglich sind, mu8 aufer (77) nun auch noch 
es ae (97) 


sein, und das heiBt zufolge (4) und (77), daB die beiden Komponenten M, und M, des Momentes 
Nt der erregenden Drehkraft in dem festen Verhaltnis 


Me Cy Cla Bl 
Mz a) =+/$4=3 (98) 
stehen miissen. Zwar fallen die Projektionen des in der Ebene (77) wandernden Drehvektors 0 
und des kérperfesten Momentvektors Jt auf die (x, z)-Ebene nicht in die gleiche Gerade, 
aber sie liegen, wie man leicht erkennt, stets im gleichen Quadranten der (x, z)-Ebene, und 
zwar wegen C < A die Projektion des Momentvektors stets naher bei der x-Achse als die- 
jenige des Drehvektors. 
Man fiihrt auBer den Abkiirzungen (80) noch die weitere dimensionslose Grobe 
De (99) 


4 


1 R.Grammel, Ing.-Arch, 21 (1953) S. 152. 
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“ 


ein und erhalt dann aus (96) mit (97) und (99) (da hier wegen m, = 0 oder m, + 0 sicherlich 
fa Vist) 


y= (100) 
und hiermit aus (95) 
SS LEP gin) cei atoms) (101) 


als Erweiterungen von (78) und (79). ‘ 
Wendet man auf die nichtlineare Differentialgleichung zweiter Ordnung (101) die Sub- 


stitution 
f(t) =f (t) off) +P (102) 


an, so geht sie iiber in eine ebensolche erster Ordnung 
d 
Get ry ames oe (103) 


Sei v = v(f) ihr Integral, etwa fiir die Anfangsbedingungen f(0) =f) und g(0) = go, so folgt 
aus (102) 


f 

eee ek 104. 

; eae ae) 
B 


hieraus durch Funktionsumkehrung f = f(r) und dann vollends aus (100) wegen (102) 
A 
arr ll (105) 


[ Obwohl die Funktionen g und v sich nur um den festen Faktor A/y unterscheiden, ist es zweck- 
maBig, in der Substitution (102) lieber v statt (y/A) g zu verwenden, weil dann die Differential - 
gleichung (103), deren — mdglicherweise nur numerisch durchfiihrbare — Integration das 
Keinstiick der ganzen Liésung bildet, nur drei statt vier Parameter enthalt. | 

Wahrend man sich in Ziff. 6 auf die Anfangswerte f(0) = 1, g(0) =0 beschranken durfte, 
weil dort sowohl fiir 2 > 0 wie fir A < 0 die beiden Halbebenen f > 0 und f< 0 und ebenso 
die beiden Halbebenen g >0 und g< 0 unter sich véllig gleichberechtigt waren und also 
sich gegenseitig an der g- bzw. f-Achse spiegelten, so muB man hier die al] gemeineren Anfangs- 
werte f(0) =f, v(0) =v» und entsprechend dann nach (105) g(0) = gy =(A/g)vp zulassen, 
und somit stellt jetzt wp) im allgemeinen nicht mehr den Anfangswert von @, dar. 

Die Differential gleichung (103) liefert in Verbindung mit (105) (und wieder mit dem MabBb- 
stabsfaktor yl + A? auf der f-Achse) die gesuchten Polhahnen. Diese haben in den Punkten 


mit den Abszissen 
fate (r > 0) 
waagerechte Tangenten und schneiden die Hyperbel 
UE Ser =p = 0 


mit lotrechter Tangente. Fiir r< 0 haben die Polbahnen keine waagerechten Tangenten, 
fiir r = 0 schneiden sie die g-Achse mit waagerechter Wendetangente. Fiir 


fy=+ 2. g =F 4p (r > 0) (106) 


schrumpfen die Polbahnen auf je einen Punkt zusammen, und dieser stellt die zugehdérige 
permanente Drehung dar, und zwar stets eine 1a bile fiir p+ 0. 

Weil wir wieder my = 0 voraussetzen diirfen, so liegen gemaB (80) die Polbahnen r > 0 
in der Ebene J, diejenigen fiir r << 0 in der Polebene IJ und sind in den beiden Ebenen wieder 
von ganz verschiedenem Aussehen. Lediglich fiir r —0 haben sie in den beiden Ebenen 
gleiche Gestalt. 

Da die ,,]abilen‘* Punkte (106) alle in der Ebene I liegen, so sind in der Ebene I fir r > 0 
lauter asymptotische Bewegungen zweiter Art zu erwarten, welche 
in den Punkten (106) asymptotisch miinden oder aus ihnen asymptotisch hervorgehen. Fiir 
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r = 0 und in der der Ebene II fiir r < 0 sind dagegen lauter asymptotische Be- 
wegungen erster Art zu erwarten. Ich lasse die Frage (als Vermutung) offen, ob 
es in der Ebene I analog zu Ziff. 5 auch zwei Typen von asymptotischen Bewegungen zweiter 
Art gibt, namlich solche, die von einem Punkt (106) asymptotisch ausgehend im Gegenpunkt 
(106) asymptotisch miinden, oder solche, die dann in eine asymptotische Bewegung erster Art 
tibergehen. Die Mannigfaltigkeit aller dieser Polbahnen ist, da die Differential gleichung (103) 
fiir einen gegebenen Kreisel («, 6, y) zwei wesentliche Parameter p und r enthalt, sehr groB 
und soll daher in einer spateren Mitteilung untersucht werden. 


8. Achsenfeste zeitabhingige Drehungen bei achsenfester zeitunabhangiger Selbsterregung. 
_ Wenn man von den trivialen Fallen der Drehungen um eine der drei Hauptachsen absieht, 
so war in den vorangehenden Lisungen (Ziff. 2 bis 7) keine einzige enthalten, bei welcher eine 
Selbsterregung mit achsenfestem zeitunabhaingigem  Momentvektor I 
eine achsenfeste Drehung hatte erzeugen oder vernichten oder auch nur in ihrer Dreh- 
geschwindigkeit 0 beeinflussen kénnen. Wir vermuten, daB das iiberhaupt nicht méglich ist, 
und wollen diese Vermutung nun beweisen. 

Jede zeitabhangige Drehung p um eine kérperfeste und also auch raumfeste Achse ]aBt sich 
in der Form 

Ox =aY(t), wy, =Oog(t), wo, =cg(t), p(t) = konst. (107) 

ansetzen, wobei a, 6, c Festwerte sind. Wir fiihren zwei kérperfeste Vektoren D und D ein mit 
den Komponenten 


D=(a,b,c), D=(Aa, Bb, Co). (108) 


Diese beiden Vektoren fallen beim unsymmetrischen Kreisel wegen (2) dann und nur dann in 
die gleiche Richtung, wenn zwei der drei Festwerte a, 6, c verschwinden, wenn also die Drehung 
um eine Hauptachse erfolgt. SchlieBen wir dies weiterhin aus, so definieren die beiden Vektoren 
Dund D eine kérperfeste Ebene E, und somit auch einen nichtverschwindenden Produktvektor 


E== [dD], (109) 


welcher auf den Vektoren D und D senkrecht steht und also insbesondere zusammen mit dem 
Vektor D eine weitere kérperfeste Ebene E, definiert, welche den Vektor D nicht enthalt. 

Nunmehr kann man fiir einen kérperfesten Drehvektor 0, der die Komponenten (107) hat, 
die Eulerschen Gleichungen (1) in der vektoriellen Gestalt 


DP) +EPM=M ue, 


schreiben. Hieraus geht hervor, daB der Momentvektor St, der zu einer Drehung (107) gehért, 
sich aus zwei nicht dauernd verschwindenden zeitabhangigen Teilen zusammensetzt, also selbst 
zeitabhangig sein mu. Damit haben wir folgenden Satz bewiesen: 

Abgesehen von den Momentvektoren in einer der drei Hauptachsen gibt es keinen kérper- 
festen zeitunabhangigen Momentvektor, der eine Drehung des Kreisels um irgendeine kérper- 
feste und also auch raumfeste Achse erzeugen oder vernichten oder auch nur in ihrer Dreh- 
geschwindigkeit 0 beeinflussen kann. 

Wir wollen dem vektoriellen Beweis dieses Satzes noch einen analytischen Beweis hinzu- 
fiigen, welcher zeigen wird, daB es zu kérperfesten zeitunabhingigen Momentvektoren Jt 
wohl Lésungen (107) gibt, daB diese aber ebenso wie Net selbst nicht reell sein kénnen. Zu diesem 
Zweck gehen wir mit (107) in die Gleichungen (5) bis (7) ein und erhalten 


ag—abcy*=™m,, 
bo +Bpecag?=m,,, (111) 
cp—yabg=m, i 

Damit diese drei Gleichungen miteinander vertraglich sind, d. h. dieselbe Differentialgleichung 


fiir die Funktion 9g (¢) darstellen, mu8 mit zwei Konstanten yw und x 


a b Cc 


und 


Se Be pea Ss (113) 
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sein, wofiir man auch schreiben kann 
CDG = FEC) « 
pea=—xb, (114) 


01D terete 
Durch Multiplikation dieser drei Gleichungen entsteht die Beziehung 
apy a bc? =— abe. (115) 


Die drei Gleichungen (114) sind fiir den ausgeschlossenen Fall der Drehung um eine der drei 
Hauptachsen infolge eines Momentvektors in der betreffenden Hauptachse erfiillt [z. B. b = 0, 
c =0und dann nach (111) auch m, =m,= 0 und nach (113) x =0]. Der Fall, daB nur eine der 
drei GréBen a, b, c verschwindet, fiihrt auf einen Widerspruch mit der Voraussetzung 
p == konst. im Ansatz (107) [z. B. gibt a =0, 6 +0, ¢ +0 nach (111) sofort py = konst. |. 
Somit bleibt nur noch der Fall a + 0, 6 + 0, ¢ + 0 tibrig. Fir ihn folgt aus (115) 


aByabe=— # 
und daraus durch Division mit der ersten oder zweiten oder dritten Gleichung (114) 
ipo Joy palate a ose 
Ge = eee b ss c oe (116) 


Wenn ~ reell ist, so miissen also a und c imaginar und 0 reell sein und folglich nach (112) | 
entweder m, und m, imaginar und my, reell oder m, und m, reell und m, imaginar, je nachdem yu 
reell oder imaginar ist. Wenn dagegen x imaginar ist, so verhalten sich a, b, c, m,, my, m, 
gerade umgekehrt. Den Fall, daB x und 4 komplex sind, brauchen wir nicht weiter zu erértern, 
auch nicht die Lésung selbst, welche z. B. fiir reelles x und reelles w aus einer beliebigen der 
Gleichungen (111) @ (t) als reelle arc tg- oder Wy Tg-Funktion der Zeit ergibt, aber eben keine 
reelle Bewegung bedeutet. 


9. Achsenfeste zeitabhangige Drehungen bei achsenfester zeitabhingiger Selbsterregung. 
Jetzt wenden wir uns dem allgemeineren Falle zu, daB die Selbsterregung zwar immer noch 
achsenfest, aber zeitabhangig ist, daB also M,, My und M, von der Form (8) mit F (¢) == konst. 
sind. Dann gibt es fiir a + 0, b + 0, ¢ + 0 eine einzige Méglichkeit, die Vektorgleichung (110) 
fiir eine achsenfeste zeitabhangige Drehung (107) und einen achsenfesten zeitabhangigen 
Momentvektor St zu erfiillen: Man muf namlich nur verlangen, daf in (110) und (8) mit einer 
reellen Konstanten 0 


P=0Ff" (117) 


FQ) = g(t) (118) 
sein soll. Die Konstante @ darf man als von Null verschieden voraussetzen, weil 0 = 0 zufol ge (117) 
auf den Fall  (t) = konst. fiihren wiirde, welcher die schon anderweitig erledigten stationaren 
Bewegungen des Kreisels bedeutet. Und in (118) kann man eine an sich mégliche multiplikative 
Konstante unbeschadet der Allgemeinheit weglassen, weil diese schon in den Konstanten M,, 
M,, M; von (8) als gemeinsamer Faktor angesehen werden diirfte und fiir das Endergebnis, wie 
man leicht erkennt, ohne Wirkung wire. 

Die Gleichheiten (117) und (118) wandeln (110) oder lieber gleich wieder die Komponenten- 
gleichungen (1) mit (8) und (107) um in 


0 Aa—(B—C)bc=M,, 


0 Bb+(A—C)ca=M,, (119) 
o Chey aed 


und 


wenn wir den ganzlich trivialen Fall g (t) = 0 der Ruhe [o0 = 0 und also nach (110) We 07 
ausschlieBen. 

Diese drei Gleichungen sind zwar bei vorgegebenen Werten M,, M,, M, algebraisch nicht 
ohne weiteres nach a, b,c aufzulésen (jedenfalls wollen wir dies hier nicht tun); wohl aber 


folgen aus ihnen bei vorgeschriebenen Werten von a, b, c und o ohne weiteres die Moment- 
konstanten M,, M,, M3. 
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Aus (117) ergibt sich schlieBlich, wenn wir der Zeit t = 0 den Funktionswert ~ = oo zu- 
ordnen, was, wie wir sehen werden, ohne wesentliche Kinschrankung der Allgemeinheit erlaubt 
ist, 

p(t) =——= (120 
Gute ) 

Man kann das Ergebnis, das in den Forme!n (120) samt (107) und in (119) samt (8) und (118) 
enthalten ist, in dem folgenden Satz zusammenfassen: 

Eine achsenfeste zeitabhangige Drehung des Kreisels mit einem kérper- und also auch 
raumfesten Drehvektor 0, der die Komponenten 


b c 
Wx ane Wy Sacre Oz eaep ir (= 0) (121) 


besitzt, kann durch eine achsenfeste zeitabhangige Selbsterregung mit einem (seiner Lage nach 
korperfesten, aber im allgemeinen nicht raumfesten) Momentvektor Mt, der die Komponenten 


1 

OE Re Ba Obie as 
1 

ERR ark Cece OK Norra (122) 
WV 

M, =[9 Co—(A— B)ab] 5, 


besitzt, unterhalten werden, und die Momente mit den Komponenten (122) sind — abgesehen 
von einer unwesentlichen Zeittransformation t = t’ — t, — die einzigen zeitabhangigen, ihrer 
Lage nach kérperfesten Momente, die bei einem unsymmetrischen Kreisel eine zeitabhangige 
Drehung um eine kérper- und raumfeste Achse, die keine Hauptachse ist, zu unterhalten 
vermégen. 

Man bestatigt auch schon durch Kinsetzen von (121) und (122) in die Eulerschen Gleichun- 
gen (1) die Richtigkeit des ersten Teiles dieses Satzes; der vorgetragene Beweis zeigt aber ins- 
besondere die Giiltigkeit des zweiten Teiles des Satzes (Einzigkeit jener Momente). Im iibrigen 
mu betont werden, daB die Momente (122) niemals den Kreisel aus vélliger Ruhe in Bewegung 
zu setzen oder eine schon vorhandene Bewegung vdéllig zur Ruhe zu bringen vermégen, weil 
diese Ruhezustande zu t = + oo gehéren: es gibt zufolge des zweiten Teiles unseres Satzes 
iiberhaupt kein kérperfestes Moment 9)¢, das fiir eine Drehung um eine kérperfeste Achse, die 
keine Hauptachse ist, hierzu in der Lage ware. 

Das Ergebnis (122) 1aBt eine einfache Deutung zu, die wir fiir die Annahme formulieren 
wollen, da8 der Anfangspunkt der Zeit bei einem hinreichend stark negativen Wert t, liegt, und 
daB die Zeit dann bis in die Nahe vont = 0 zunimmt. In diesem Falle stellt der erste Bestand- 
teil It, = Dep des Momentvektors (110) mit seinen Komponenten 
Aa eee) aCe 
ot ° i 0 ’ Ds ow ° 


Nes (123) 


die die ersten Glieder rechts in (122) sind, das Antriebsmoment vor, welches von einer kleinen, 
schon vorhandenen Anfangsdrehung 


a b c 
- set ie 124 
os pn eee ae oar oe 
aus die Drehgeschwindigkeit 9 des Kreisels immer mehr vergréSert und fiir t> 0 itiber alle 
Grenzen wachsen lat. Der zweite Bestandteil Jt, = & y? des Momentvektors (110) mit 
seinen Komponenten 


(B—C)be 
Mas = a 5 ee 

A—C 
My, = +o ee (125) 
M.. (ae 5) obs 


Oi 0? 
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die die zweiten Glieder rechts in (122) sind, stellt das Moment vor, welches notig ist, um das bei 
dieser beschleunigten Drehung geweckte Flichkraftmoment zu kompensieren und so die Dreh- 
achse in ihrer Lage festzuhalten. 

Bei positiven Werten von a, 6, c und @ liegt fiir t< 0 der Drehvektor 9 in dem Oktanten 
(x >0, y >0, > 0), ebenso der Momentvektor St, (123), wogegen der Momentvektor mM, 
(125) in dem Oktanten (x< 0, y > 0, 2< 0) liegt und sowohl auf dem Drehvektor D0 wie aut 
dem Momentvektor Jt, senkrecht steht, weil 0 richtungsgleich mit D, ferner Mt, richtungsgleich 
mit D und Jt, richtungsgleich mit & ist, und weil zufolge (109) € auf D und D senkrecht steht. 
In welchem Oktanten der resultierende Momentvektor It = Nt, + Mt, liegt, das hangt von den 
Werten a, b,c, A, B, C und @ ab. 

Wenn der Drehvektor in einem anderen Oktanten liegt, so liegt auch der erste Momentvek- 
tor Nt, im gleichen Oktanten wie 0, wogegen die oktantenmaBige Zuordnung zwischen 9 und 
dem zweiten Momentvektor Jt, genau die gleiche ist, die ich schon fiir die stationaren Drehun- 
gen und ihr unterhaltendes Moment gefunden habe }. 

Fiir positive t bedeutet die Lésung (121), (122) das Abbremsen einer kérper- und raumfesten 
Drehung durch ein kérperfestes zeitabhangiges Moment, wobei aber, wie schon erwahnt, die 
Drehgeschwindigkeit 0 erst mit t—> oo gegen Null geht. 

Will man negative Zeiten t beim Anlaufen des Kreisels vermeiden, so hat man lediglich die 
Zeit t in allen Formeln durch eine Zeit t/ vermittels der Transformation t = t’ — t, mit hin- 
reichend groBem positivem Festwert t,; zu ersetzen, und man erkennt dabei, daf in der Tat 
die bei der Integration von (117) beniitzte Zuordnung von t = 0 und m = oo ganz unwesentlich 
ist und nur zur Vereinfachung der Formeln diente. 


10. Achsenfeste zeitabhingige Drehungen bei achsenbeweglicher Selbsterregung. Jetzt 
bleibt noch die Frage ttbrig, wie eine achsenfeste zeitabhang’ge Drehung — da sie, abgesehen 
von den trivialen Fallen der Drehungen um eine der drei Hauptachsen und abgesehen von dem 
Fall von Ziff. 9 durch kein achsenfestes Moment 9)¢ beeinfluBbar ist und, wieder abgesehen von 
jenen trivialen Fallen, insbesondere nicht durch ein achsenfestes Moment St erzeugt oder ver- 
nichtet werden kann — wenigstens durch ein achsenbewegliches Moment beein- 
fluBt und insbesondere erzeugt oder vernichtet zu werden vermag. Die Antwort auf diese Frage 
ist schon in der Forme] (110) vollstandig enthalten und kann in dem folgenden Satz aus- 
gesprochen werden (der fiir praktische Anwendungen offenbar wichtig ist): 


Abgesehen von den trivialen Fallen der Drehungen um eine der drei Hauptachsen kann eine 
Drehung um irgendeine kérperfeste und also auch raumfeste Achse nur durch eine achsen- 
bewegliche Selbsterregung erzeugt oder vernichtet oder allgemein zeitlich anders regu- 
liert werden als die Drehungen von Ziff. 9, und zwar wandert dabei der Momentvektor Jt nicht 
beliebig im Kreisel, sondern lediglich in einer kérperfesten Ebene E,, welche durch die Vektoren 
®D (108) und € (109) bestimmt ist. 

Der erste Bestandteil DY von Yt (110) muB dabei, wie in Ziff. 9, die Drehgeschwindigkeit 
regulieren; der zweite Bestandteil © g* muB das Flichkraftmoment kompensieren, und seine 
oktantenmaBige Zuordnung zum Drehvektor 0 ist auch in diesem allgemeinen Falle wieder die 
gleiche wie bei den stationdéren Drehungen. 


Wenn insbesondere die Drehung 0 von der Ruhe aus gleichférmig beschleunigt anlaufen 
soll, wenn also  (t) = 1 ist, so wird der erste Bestandteil D@ von Yt der konstante (kérper- 
feste) Vektor D, wogegen der zweite Bestandteil © g? = € # mit konstanter Beschleunigung 
zunimmt. Die Pfeilspitze des Momentvektors Jt lauft dann also mit konstanter Beschleunigung 
2 € auf einer kérperfesten Geraden, die parallel zum Vektor © durch die Pfeilspitze des Vektors 
® geht, und der Vektor 9)t nahert sich mit unbegrenzt anwachsendem Betrag asymptotisch der 
Richtung des Vektors ©, der senkrecht zur Drehachse steht. 


11. Zusammenfassung. Wir kénnen alle bisherigen Erkenntnisse iiber die Selbsterregung 
des unsymmetrischen Kreisels dahin zusammenfassen, da — abgesehen von den trivialen Fallen 
der Drehungen um jede der drei Hauptachsen und abgesehen von den permanenten Drehungen 
um beliebige Achsen— zu jeder achsenfesten zeitunabhangigen Sel bsterregung achsenbewegliche 
Drehungen gehéren, wie wir solche in Ziff. 2 bis 7 untersucht haben; ferner daB es, wieder 
abgesehen von jenen trivialen Fallen, nur den einzigen Typ der achsenfesten zeitabhangigen 


1 R. Grammel, Ing.-Arch, 21 (1953), S. 151. 
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Selbsterregung von Ziff.9 gibt, welche eine achsenfeste zeitabhangige Drehung zu_beein- 
flussen, wenn auch weder zu erzeugen noch zu vernichten vermag; endlich daB der Moment- 
vektor Jt der Selbsterregung, welcher eine beliebige achsenfeste Drehung erzeugen oder ver- 
nichten oder in beliebiger Weise beeinflussen soll, im allgemeinen (also abgesehen von den 
trivialen Fallen und von dem Fall in Ziff. 9) zeitabhangig in einer kérperfesten Ebene wandern 
mu. Die achsenbeweglichen Drehungen, welche von achsenfester zeitunabhangiger Selbst- 
erregung erzeugt oder unterhalten werden, sind je nach der Lage (und Gréfe) des Momentes Nt 
der Selbsterregung periodisch hinsichtlich der Bewegung des Drehvektors 0 im Kreisel oder 
asymptotisch mit unbeschrankt wachsendem Drehvektor 0 oder asymptotisch mit einem 
Drehvektor 0, der sich einem permanenten Drehvektor (von endlichem Betrag) annahert 
(oder auch sich von ihm entfernt). 


(Hingegangen am 23, Januar 1953.) 
Anschrift des Verfassers: Prof. Dr. Richard Grammel, Stuttgart-N., Robert-Bosch-Strafe 101. 
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Theorie der elastischen Ringe starker Kriimmung. 
Von F. K. G, Odqvist. 


1. Einleitung. Die technische Theorie stark gekriimmter Stabe geht auf die ersten Jahr- 
zehnte des vorigen Jahrhunderts zuriick und zahit als Urheber Navier, Bresse und Resal. Eine 
vollstandige Behandlung findet sich bei EF. Winkler und wurde spater in den meisten Lehr- 
biichern der Festigkeitslehre abgedruckt +. Hier handelt es sich ausschlieBlich um Stabe mit in 
bezug auf die Schwerpunktsebene symmetrischen Querschnitt. Die duBere Belastung findet sich 
auch in dieser Ebene. 

Neuere Untersuchungen der Stabe starker Kriimmung behandeln meistens nur spezielle 
Belastungsfalle. Eine Reihe von Untersuchungen kniipfen an eine Arbeit von O. Géhner iiber 
Spiralfedern ? an. Wichtige Arbeiten insbesondere tiber geschlossene Ringe und deren Um- 
still pung rithren von R. Grammel her *. 

Im Folgenden wird das Problem des stark gekriimmten, offenen oder geschlossenen Kreis- 
ringes vom belie bigen Querschnitt neu aufgegriffen, teils um allgemeinere Lastsysteme 
zuzulassen, teils um eines der Ergebnisse von R. Grammel noch ein wenig zu erweitern. 


2. Die Gleichgewichtshedingungen. Um die Geometrie des gekriimmten Stabes in der Form 
eines Kreisringstiickes vom Halbmesser r zu beschreiben, wird ein bewegliches Koordinaten- 
system laut Abb. 1 eingefitthrt mit x =r@p als Bogenlange der Schwerpunktslinie, y in radialer 

Z Richtung positiv nach aufen und z in achsialer Richtung nach 
iF oben. Am Bogenelement greifen die SchnittgréBen My, M, (Biege- 
momente), M, (Torsionsmoment), Q, T (Schubkrafte in axialer 
bzw. radialer Richtung) und N (Normalkraft) an (Abb. 2 und 3). 

Die streckenweise stetig verteilten dauBeren Krafte werden 
je Langeneinheit der Schwerpunktslinie aufgelést in p, q, t (in 
radialer, bzw. axialer und tangentialer Richtung). Dazu wird 
noch ein stetig verteiltes Drillingsmoment m je Langeneinheit 


ks 


Abb.1. Allgemeines Koordinaten- Abb, 2. Krafte am Abb. 3. Momentenzeiger am 
system fiir ein Kreisringstiick. Bogenelement. Bogenelement. 


angenommen. Sodann Jauten die Gleichgewichtsbedingungen am Bogenelement 


dN+ Tdp+trdy =—0, (1) 
dT— Ndy+prdy=0, (2) 
dQ+qrdp=0, (3) 
dM,+ M,dp+mrdp=0, (4) 
dM, — M,dp— Qrdp =0, (5) 
dM,— Trdp = 0. (6) 


1 Siehe I. Todhunter u. K. Pearson, A History of the Theory of Elasticity ..., Bd. II, 1, 5. 422, Cam- 
bridge 1893. j 

2 Siehe C. B. Biezeno u. R. Grammel, Technische Dynamik, S, 305—404, Berlin 1939, in der nach Ab- 
schlu8 des Manuskripts erschienenen 2. Auflage: Bd. I, 5. 339—438. 
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Hier ist r der Kriimmungshalbmesser der Schwerpunktslinie. Es wird wie allgemein in der Statik 
die Verzerrung des Stabes so klein angenommen, daB die SchnittgréBen ebenso gut auf das ver- 
zerrte wie auf das unverzerrte Element bezogen werden kénnen. 


Im statisch bestimmten Falle lassen sich die Gleichungen (1) bis (6) unmittelbar durch 
Quadraturen integrieren. Insofern unterscheidet sich der von uns behandelte Fall des stark 
gekriimmten Ringes keineswegs von dem von Biezeno und Grammel behandelten Fall des 
schwach gekriimmten Ringes. Wir wollen jedoch die Integration etwas verschieden durch- 
fiihren. Wenn die AnfangsgréBen fiir y = 0 der betreffenden SchnittgréBen mit iy Peal BST dP 
Mio, My. und M,, bezeichnet werden, bekommen wir mit elementaren Methoden 


Q 
N = N, cos » — T, sin p + [ [prsin (p — A)—trcos(p— A)] dad, (7) 
0 
¢ 
T = T, cos g + Ny sin y — f [pros (py — A) + trsin (p—A)] dA, (8) 
3 
~ 
Q=—fard, (9) 
0 


Pp Pp 
M, = (Mi) + Qy 7) cos p— Myosin p — Qor + [ gr? dA— | (mr + qr’) cos (p—A) di, (10) 
0 0 


7 
M, = (Mio + Qo 7) sin gp + M,. cos — [{ (m r-+ qr’) sin (p — A) dd, (11) 


0 


M, = M,,+ Tyrsing + Nr (1— cos g)— | {presin (p—A) + tr? [1—cos (p—A)]}} dd. (12) 
6 


In diesen Integralen sind p, q, t und m als Funktionen der Integrationsverdnderlichen / auf- 
zufassen. Die Liésungen (7) bis (12) lassen sich in zwei Gruppen zerlegen (7), (8), (12) bzw. (9), 
(10), (11), wobei die Anfangswerte, die in je einer Gruppe auftreten, unabhangig von denjenigen 
der anderen angenommen werden kénnen. Belastungsfalle der einen Gruppe lassen sich daher 
unabhangig von denjenigen der anderen liésen, jedenfalls sofern es sich um statisch bestimmte 
Lastfalle handelt. 


3. Deformationshypothese, Spannungsverteilung. Um die innere Spannungsverteilung des 
Stabes zu bestimmen, mu man etwa wie in der technischen Biegetheorie gerader Balken auf 
die Deformation naher eingehen. Fiir die Ringdehnung ¢ in der x-Richtung geniigt hierbei die 
Bernoullische Annahme des Ebenbleibens der Querschnitte. Der Verschiebungsvektor der 


Abb. 4. Verschiebungsvektor. Abb. 5. Winkelainderungen. am Bogenelement. 


Schwerpunktslinie sei mit u, v, w bezeichnet (Abb. 4). Wenn Differentiationen in bezug auf 
mit Strichen angedeutet werden, bekommen wir fiir die Winkelabweichung y der Tangente der 
Schwerpunktslinie aus ihrer Ebene 


wu’ 


: Sire 
Die Ringdehnung der Schwerpunktslinie selbst wird sein 


du v uw +ov 


Fie r 


Wenn der Stab dazu noch eine Drillung um die Schwerpunktslinie um einen Winkel # erfahrt 
(positiv im Uhrzeigersinn, wenn in der Richtung positiver p betrachtet; Abb. 5), dann wird sich 
eine beliebige Ringfaser von der Lange (r + 7) dp beim Flachenelement dF des Querschnittes 
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um den Betrag 


Sanh oe t 8 dp —Cdy—rdon 


dehnen als Folge der Verschiebung und der Winkel anderungen @ und dy. 
Die Bernoullische Hypothese driickt sich somit durch die folgende Gleichung aus: 
1 v 


c= |W too + 9) n—|- (13) 


Unter Heranziehung des Hookeschen Gesetzes bekommt man sodann fiir die Spannungsver- 


teilung 
o=Ese (14) 


unter der Voraussetzung, daf u, v, w und & bekannte Funktionen von ¢ sind. 


4. Die Grundgleichungen. Die Aquivalenz auferer und innerer Krafte beim Querschnitt 


ergibt 


N=/fodF, (15) 
Mpa (ot dk; (16) 
M, =fondF. (17) 
Wir fiihren die Bezeichnungen 
[ake be (18) 
CndF 
co . (19) 
ny? dF rs 
| 1 + y/r 3 A, 
acid La F 
| eae (21) 
ein und bekommen alsdann 
dF I dj 
[n=s(F+ 4): (22) 
C1 DD: 
r+n r2 * (23) 
nade a] 
as (24) 


Wird schlieBlich 


5 lI D)=J,, | 
= (IJ — D*) = Jy, (25) 
= (IJ — D*) = J; | 


gesetzt, so bekommen wir nach Einfithrung von (13) und (14) in (15) bis (17 t 
von (18) bis (24) (13) (14) in (15) bis (17) unter Benutzung 


Pawo) FPP) +r Oty) D+ws (26) 
Mo? ’ 

oT = —(w +0) D—r(O+y)I—o"D, (27) 
Ms? 


E Sa (Ue) Pr (ni ap) Beet 2 3 (28) 
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Wird nach den DeformationsgréBen aufgelést so bekommen wir 


Nr+M; 
wteao (29) 
ee Mor M:r 
je. Mer M; r? Nr-+ M; 
Coes eee Pp Pree. (31) 


Wenn (29) bis (31) in (13) und (14) eingetragen werden, so ersieht man sofort, daB die Verteilung 
der Normalspannung durch die Werte von N, M, und M, allein bedingt ist [vg]. spater Glei- 
chung (33’)]. 

Aus der Abb. 5 ist ersichtlich, da der Zuwachs d? — y dq des Drillungswinkels am Element 
r dy der Schwerpunktslinie durch das Torsionsmoment M, bedingt sein mu. Wir setzen dem- 
nach an 


(32) 


wo in erster Annaherung GK als die Verdrehungssteifigkeit des geraden Stabes vom selben 
Querschnitt angenommen werden kann. 

Die in (32) steckende Hypothese ist zweifellos die fragwiirdigste der ganzen Theorie. Sie 
wird im Folgenden nur deswegen eingefiihrt weil wir eine Kupplungsbedingung zwischen # und 
wy notig haben. Noch schlechter begriindet ware eine nachtragliche Bestimmung der Schub- 
spannungen iiber den Querschnitt aus M, und der St-Venantschen Theorie der Verdrehung 
prismatischer Stabe. Nachdem die rechterhand in den Gleichungen (29) bis (32) auftretenden 
GréBen laut (7) bis (12) bekannt sind, ersieht man sofort, daB wir ein vollstandiges System von 
Grundgleichungen zur Bestimmung der Gréfen u, v, w und @ als Funktionen von y bekommen 


haben. 


Der Kiirze halber fithren wir im Folgenden die Bezeichnungen 


1 r r ip r 
ep a i ee ey oY (33) 


ein und wollen uns auf den Fall beschranken, da «, B,, By, Bz und y konstant sind, obgleich dies 
an sich nicht notwendig ware. Mit diesen Bezeichnungen ergibt sich nach (14) und (29) bis (31) 
fiir die Spannungsverteilung im Querschnitt 


o =| al + M,) —5(— Ba My + Bx Me) — | Ba My — Ba Ma — 5 (Nr + M,)]| (384 


ary 


5. Integration der Grundgleichungen. Als Anfangsbedingungen fiir y = 0 fiir die Defor- 
mationsgréBen werden angenommen 


u(0) =u, (0) =%, v(0) =m, w(0) =m, w(0) = (0) =Y, H0) =H. (34) 


Wir fangen mit (31) an und bekommen mit Riicksicht auf (33) nach ein- und zweimaliger Inte- 
gration 


g. WA P 
v = + 8,7 | M,dA— («+ Bsr) f M.dA—arf Ndi, (35) 
0 0 


0 
gQ Ay P Ay P Ay 
v=y+tup+pirf [M,dAda,— («+ Bsr) f [ M.dAdjy—arf f[Ndadi. (36) 
0 0 0 0 0 0 
Die Integrale in (35) und (36) kénnen mit Hilfe von (7), (11) und (12) berechnet werden und 
durch Benutzung der Integralformel von Lejewne-Dirichlet vereinfacht werden. Fir die An- 
wendung hier und im Folgenden erhalten wir so die Formelsammlung 


2 
[Na = N, sin yg —T,(1 — cos 9) 
6 


+ f{pr[1 —cos(p—A)] —trsin (p—A} di, (37) 
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pg A 


i i, N da dd, = N, (1— cos y) — Ty(yv — sin 9) 
0 


+ [ip r[p —A—sin (p —A)] —t r[1 — cos (p — A} a, (38) 
0 


aah, 
ff JN diddy dry = No(p — sing) — T($ g? — 1 + €08 9) 


000 


Jip Tithe A) Le cos (p—A)]—tr[p —A—sin ie he 


JM, dA = (M,. + Q)r) (1 — cos y) + Myosin 


— fiumr + 479) [1 — cos (p—A)] aA, (40) 
0 


g, Aa 
Re dd dA, = (Min + Q 7) (vy — sin y) + My, (1 — cos ?) 


0 0 


Mm grip en ee (41) 


A 
| M, dA dd, dd, = (Mi. + Qo r) (4 ? — 14 cos ~) + Moo (Y — Sin ~) 
0 


— fiimr + ar) 1b (9 — A? — 1 + 008 (p — AN) aA, (42) 
0 
J M.di = M oP + Tyr (1 — cos y) + Nor (vy — sin @) 
—— [fp 2 [1 — cos (p — A) Ft? [p —A— sin (P—A} AU, (43) 
0 


eA 
J [ -M. a dd, =4 May gp +1 (p — sin p) + Nor (yy? — 1 + cos 9) 
00 


P 
— f{pr[y—A—sin (yp —A)] +t [s (p—AP—1 
, : + eos (p— Ay dd. am 
"f° M, dd da, dig =} Mig? + Tor (4g? — 1+ 008 9) + Nor (kp? —@ + sing) 
0 


j 


Ss 


c—} 


= Jern (p — A — 1+ cos (p— A] + tr? [k (p —A)$ 
—(p—A)+sin(p—A)}} dA. (45) 
Tm nichsten Schritt bestimmt man u aus (29): 


YP QP 
u=u—foditarfNdi+af Md, (46) 
0 0 0 


worin die letzten beiden Integrale durch (37) und (43) gegeben sind und das erste Integral 
folgendermafen zu bilden ist: 


foda=not4 v9 + Byr ff fm, dj dA, di, 
0 0 0 
avis 
~ tot bof fia dna aya ff faa dy (a7 
0 0 0: 20)50 


Die hier auftretenden neuen Integrale sind unter den oben berechneten enthalten, namlich (42), 
(45) und (39). 
Aus (30) und (32) erhalten wir 


yp’ +y=—yM,——M+8,M,, (48) 
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eine Differential gleichung, deren Integral nach Befriedigung der Anfangsbedingung 
Y = Yo cos p — (% + fo M,, — B, M 20) sing — 9 | Mz sin (y — 2) dj 
B. | My sin (p — A) da + B, i M, sin (p — A) dd (49) 
0 0 

lautet. Gleichzeitig bekommen wir 
8 = (Io + Ba Myo — fr Mi } eos.» + yosin yp — 7 ['M,c0s (p — 2) dp 

+ a M, cos (p — A) dA — psf M; cos (py — A) dA — B, M, + B, M,. (50) 
Durch partielle Integration ergibt sich 


9 g 
x, sin (p — A) d4 = — vosing + { M, cos (p — A) dA (51) 
o 0 


und 


9 
fey cos (yp — A) dA = M, — Myo cos 9 — | M, sin (p — A) dA (52) 


und zwei ahnliche Formeln fiir M, anstatt M,. Durch Benutzung von (51) und (52) lassen sich 
(49) und (50) vereinfachen, so daB kommt 


py = yo cos p — Dysin p—y [ M, sin (p — A) dd 
0 


es Micon (ee) di fi) M con = Ayal, (53) 


Q 
B= 3, cosm + yy sing + y { M, cos (gy — A) dA 
0 
P 
— fi. f M, sin (p A) da + By f M, sin (p — A) dA. (54) 
0 0 


Wenn (10) bis (12) in (53) und (54) eingetragen werden, bekommt man fiir die dort auftretenden 
Integrale 


fim, sin (p — A) dA = 4 (M,, + Qo r) (pv sin y) — 4 My, (sin » — ¢ cos @) 
é 
— Q,r (1 — cos y) + ic r? [1 — cos (y — A)] dA 
0 


(mr + q7°) [(p — A) sin (p — A)] dd (55) 


oc 


ne (Mio + Qo 7) (sin g + y cos vy) — 4 M,, (vy sin vy) — Qyrsin 
' y 


+ [artsin(p—2) )dA—} $f (mr + qr’) [sin (g — A) 
"E(w A)eos(p— Aa, (56) 


[M, sin (p — A) dA = 4 (Mi) + Qo r) (sin y — pcos py) + $ My (v sin g) 
0 


— 4 f (mr + 14) [sin (p — a) — (py — A) 008 (p — ANI AA, (57) 


fom, ©05 (p — A) #2 = 4 (Mig + Q 7) (p sin g) +4 Myo (sin y + 9 008 9) 
0 


— 4 [mr + qr) (p—A) sin (py — A) a2, (58) 
0 
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°M, sim (p — A) dd = Mg (1 — cos @) +} Tyr (sin p — 9 €08 ) 
+ Nyr[1— cos — } (psin 9) 
— fig pe lsin (9 —2) —(P— A) c08 (P— AI] 
4 t19(1 — cos (p — A) — 3 [lp — 4) sin (9 — AD} (59) 


a . 
J M, cos (py — A) dA = M,ysing + > T, 1 (psin g) + 4 Nor (sing — 7 05 ) 
0 


3 Ip r? [(p — A) sin (p — )] 
0 
+ tr® [sin (p— A) — (p — A) cos (y — Aj} dA. (60) 


SchlieBlich bekommen wir 


? 
w=utrfyd (61) 
0 
oder mit Hilfe von (53) und partieller Integration 


Pp Ai 
B= 2+ yy sin p — Bp (1 — cos p) — adele sin (A, — A) dd 
0 


vp 


Ay 2 Ay 
— By f da, [My cos (Ay — 2) dA + By | @, | M, cos (2, — 2) @ 
0 0 0 0 
Q 
= 0+ yy sin p — (1 — cos y) — 7 | M, [1 — cos (p — A)] dA 
0 
P 
— fy My sin (p — A) dA + By f Mz sin (p — A) aA. (62) 
0 0 
Die in (62) auftretenden Integrale lassen sich durch (56), (57), (59) und die Formel 
@. 
| M, dA = (M,, + QT) sin gy — My, (1 — cos g) —Q TP 
: g, ?. 
+ [qr (p—A) d— J (mr + qr?) sin (p— A) dd (63) 
0 0 


ermitteln. Die Gleichungen (35) bis (62) bilden ein vollstandiges System von Liésungen der 
Grundgleichungen (29) bis (32). 

Bis jetzt waren die Anfangswerte Nj». Les Qes Mies Meo ize BOS die GréBen (34) (insgesamt 
zwolf GroéBen) als von einander unabhiangig gegebene GréBen angenommen. Wie wir sehen wer- 
den, sind die zwélf GréBen bei statisch unbestimmten Fallen durch weitere Bedingungen ver- 


knipft. 


Die obige Lésung schlieBt konzentrierte Lasten aus. Im Falle konzentrierter Lasten oder 
Momente wird man den Ring in Intervalle zerlegen, wobei innerhalb jedes Intervalls die obige 
Methode verwendbar sein wird mit im allgemeinen verschiedenen Werten der zwolf Integra- 
tionskonstanten N, bis M,o, ug bis wo, J) in jedem Intervall. 


6. Spezielle Lastfalle. Wir betrachten eine Reihe von sechs Lastfallen eines offenen Ring- 
stiickes vom Mittelpunktswinkel 2g, (Abb. 6), wo auBere Krafte ausschlieBlich an den End- 
flachen des Ringstiickes angreifen. Die sechs Falle finden sich in der Tabelle 1 verzeichnet, und 
die entsprechenden Werte von N,... Mz, fiir den Fall angegeben, da8 sich der Winkel p von 
der Mitte des Ringstiickes gerechnet wird. Alsdann wird die entsprechende Lésung aus (7) bis 
(12) fir p =q =t =m = 0 gewonnen. Ferner lassen sich auch die Deformationsgréfen u, v, 
w, # nach (34) bis (63) integrallos darstellen. Der Fall IV entspricht dem Falle einer zylindri- 
schen Spiralfeder kleiner Steigung mit zentraler Last. Wir werden in den nachsten Abschnitten 
einige der Lastfalle benutzen, um bestimmte Probleme von praktischer Bedeutung zu lésen. 


7. Der geschlossene Ring mit konzentrierten radialen Kraften. Wir beschranken uns auf 
den Fall gleichmaBig verteilter, gleich groBer Krafte nach Abb. 7. Durch Kombination der 
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Lastfalle III und VI der Tabelle gewinnt man sofort 


= P re 1—cos 
ore: Mz» M, - Pr Santor T, = Q = M,, = M,, =0, (64) 


worin M, den Wert von M, fiir YP = 7, eine statisch unbestimmte GréBe bedeutet. Aus (64) 
und (7) bis (12), (13), (14) und (29) bis (33’) wird die Spannungsverteilung bis auf die GréBe M, 
vollstandig bestimmt. Aus (35) bekommen wir 


, M, 2 ysing, +Pr (cosy, —sin ) sin 
ef == (x r 1 1 1 Pp 2 
ata ( 13 Bs 7) 2 sin g, se aae 2sin y,° 1) 
Tabelle 1. 


Lasf- 
tale 


-F, r(1-cos gy) 


—— 


Aus Symmetriegriinden gilt v’ = 0 fiir gy = 0 und P =, somit v9 = 0 und 


ae Pie Wey 
M, = 2 Pipe er! ° (66) 


Fiir in bezug auf die Ringebene symmetrischen Querschnitt (D = 0) reduziert sich (66) auf eine 
wohlbekannte Formel. 


Abb.7. Geschlossener Kreisring mit 
konzentrierten radialen Kriften. 
Ausgeschnittenes Teilstiick. 


Abb. 6. Koordinatensystem fiir 
ein offenes Ringstiick vom 
Mittelpunktswinkel 2 ¢,. 


Integriert man (65) zweimal, so bekommt man 


(gy? cos y,) — 1 + cos p 


1 ; 1 — cos 
v= — (+6, [2 Mg? +rPS See gD a (67) 


2sing, ” 


4 vy’ cos yp, —p +sing Y —siny 
w= — mp +t («+ Ba) /b Mp feat mos 2 sin @, arP Pang, (68) 


Die Konstanten u, und v, bestimmen sich aus der Symmetriebedingung u = 0 fiir gm = 0 und 
Y =. Somit bekommen wir u, = 0 und einen komplizierten Ausdruck fiir vy), den wir hier 


unterdriicken. Weiter ergibt sich aus (53) und (54) mit Hilfe von (55) bis (60) 


; J] — : in ~— 
yp = yw cosy — % sing + py (14 Ppa St =) sing + P po ‘ ae "| . (69) 


2 sin 9, 
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1—cos mp —} sing 


: 1— cos 9 
9 = Bycosp + yosing + fa Mi Pe. aaa ‘J (1 — 08 9) holt, 2sin gy 


(70) 


Die Konstanten y, und #, bestimmen sich aus den Bedingungen p = 0 fix p = 0 und yg = 9, 
d. h. wir bekommen y, = 0 und einen komplizierten Ausdruck fiir 9. SchlieBlich ergibt sich 


aus (62) 
(1 — cos y) —psin y 


| 1 ? P a) 
bi ¥o Iy(1 cos ?) + By (4 Pr WY, ae ‘ a tee ) f 4sin 9 


und hier wird w, aus der Bedingung berechnet, daB w = 0 fiir p = 9, ist. Somit ist der Defor- 
mationszustand vollstandig bestimmt. Aus der Lésung (64) bis (71) ist zu entnehmen, daB sich 
die Schwerpunktslinie des Ringes im allgemeinen aus der Ringebene heraus verbiegen wird, 
auch wenn die auBeren Krafte selbst auf die Schwerpunktsebene der Ringquerschnitte be- 
schrankt sind. 


8. Der geschlossene Ring mit konzentrierten umstiilpenden Momenten. Wir beschranken 
uns auf den Fall gleich groBer konzentrierter Momente M, gleichmaBig verteilt mit dem Tei- 
lungswinkel 2 y,. Wir bemerken, daB dies genau dem Grammelschen Falle entspricht und zwar 
fiir sogenannte schwache Umstiilpung. Hier 1a6t sich eine Liésung durch Kombination der 
Falle IL und III der Tabelle herstellen. Aus Gleichgewichtsgriinden gilt (man vg]. Abb. 8) 


Meee eee (72) 


aM 2sin gy, ’” 


LGTY 


Abb. 8. Geschlossener Kreisring mit konzentrierten Abb. 9. Geschlossener Kreisring mit gleichmaBig verteilter Last quer 
umstilpenden Momenten. zur Ringebene und dquidistanten Stiitzen. 


wogegen M,, statisch unbestimmt bleibt, und im tibrigen N, T, =0) = My —=0> Ausdee 
Bedingung v’ = 0 fiir py = 0 und mp = q, gewinnen wir mit Hilfe von (35) 
_ Br sing, 
smears ay De | 


Aus (53) ergibt sich ferner 


y =— Ising + $y Myo (sin p — pcos p) — } Ba Myo (sing + pcos Y) + By Mio sin y. (74) 
Die Konstante #, wird aus der Bedingung yp = 0 fiir y = —¢, bestimmt, so daB 
M 2 p2 : 
0 ~ 4sing, 6. (1 + oy, ctg y,) — y (1 — 9, ctg g1) — oe oa (75) 
wird. SchlieBlich ergibt sich der Wert 0, von @ fiir yp = y, zu 
et Y ~ 2p2r sing 
ve 4 sin Q P Ga Py eee 01) + Here Pa 1) PEG +B, Tas i ; (76) 


Der volistandige Deformationszustand wird wie im vorigen Abschnitte bestimmt aus den Be- 
dingungen u = 0 firg = 0und@ = y,, w = 0 fiir y = q,. Die Formel (76) ist fiir den Sonder- 
fall B, = 0 schon von Grammel angegeben worden }. 


1 Technische Dynamik, 2. Auflage, Bd. 1, S. 437, Formel (7). Der Fall 8, = 0 wird dort als symmetrie- 
artiger Ring bezeichnet. 
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9. Der geschlossene Ring mit gleichmifig verteilter Last quer zur Ringebene und di quidistan- 
ten Stiitzen. Dieser Fall wird durch Abb. 9 illustriert. Er ist im Falle schwach gekriimmter 
Ringe von E. Volterra! behandelt worden, wobei ein mit der Durchbiegung —w proportionaler 
Bettungsdruck q auch beriicksichtigt wurde. Wir behandeln hier den Grenzfall stark ge- 
kriimmter Ringe unter Beschrinkung auf den Fall q =konst. Wir rechnen ~ wie frither von 
der Intervallmitte positiv nach rechts in Abb. 9. Es tritt ein statisch unbestimmtes Biegemo- 
ment M, uber den Stiitzpunkten auf, hingegen verschwindet dort das Torsionsmoment aus 
Symmetriegriinden. Eine Gleichgewichtsbetrachtung ergibt sofort 


My) = 7? (sin py, — P1 cos 9) — M, sin g, , 


Myo = qr’ (g, sin y, — 1+ cos 91) — M, cos y, , (77) 
0, == 0% Neely SE = 05 


wonach sich alle SchnittgréSen wie vorher aus den Formeln (7) bis (12) ergeben. Unbekannt 
sind noch die Groen Up, LO, Wy, Yo: Jo und M,. Zuihrer Bestimmung dienen die Bedingungen 


= 3) ssa (0) 5 1) = 0, 1) 1). 
P=Piu v ~ | (78) 
yA = OF. v=0, y=), J 
also genau die geniigende Anzahl. Die Bedingung v’ = 0 fiir p = 0 ergibt v) = 0 und v’ =0 
fiir » = q, also 


J M, di = 0, 
0 
woraus sich M, berechnen 1]aBt: . 
M, = qr’? (1 — 9, etg ¢). (79) 


Wir unterdriicken die weitere Ausfiihrung der Deformationsberechnung mit der Bemerkung, 
da8 auch in diesem Falle alle drei VerschiebungsgréBen u, v, w, sowie der Verdrehungswinkel 9 
von Null verschieden sind. Die Spannungsverteilung ergibt sich wie frither durch Einfithrung 
der SchnittgréBen in (33). 


(EHingegangen am 29. Juni 1953). 
Anschrift des Verfassers: Professor Dr. F. K. G. Odgqvist, Djursholm 2, Schweden. 


1 FE, Volterra, Journ. Appl. Mech., Bd. 19, S. 1—4, 1952. 
8* 
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Uber spezielle Systeme von Tragfliigelgittern 1. 
(Theorie der tragenden Linie.) 


Von K. Nickel. 


1. Einleitung. Wahrend in einer friiheren Veréffentlichung } Systeme von endlich vielen 
Tragfliigeln untersucht wurden, soll in der vorliegenden und einer noch folgenden Arbeit der 
Grenziibergang zu unendlich vielen Tragfliigeln vollzogen werden. Eine jalousieartige Zu- 
sammenstellung von unendlich vielen gleichartigen Tragfliigeln mit gleichen Abstanden wird 
seit langem als Tragfliigelgitter bezeichnet; die Betrachtung wird im folgenden auf endlich 
viele solche Tragfliigelgitter beschrankt, denen iiberdies spiater noch weitere Bedingungen auf- 
erlegt werden sollen. Die Behandlungsmethode ist dieselbe wie dort, es wird also von reibungs- 
freier, inkompressibler und stationarer Potentialstrémung ausgegangen und der Strémungs- 
verlauf aus einer Wirbelbelegung der Tragfliigel berechnet, wie das in der Prandtlschen Theorie 
der tragenden Linie und der Birnbaumschen Theorie der diinnen Profile geschieht. Um den 
Umfang nicht zu sehr anwachsen zu lassen, werden diese beiden Teile getrennt behandelt; die 
vorliegende Note bringt die Idealisierung der Fliigelgitter als tragende Linien, waihrend iiber 
ihre Auffassung als diinne Profile demnachst an derselben Stelle berichtet werden soll. In 
beiden Fallen werde die Anzahl der betrachteten Fliigelgitter mit n bezeichnet. 

Im Gegensatz zu Gittern aus Fliigelprofilen sind die hier betrachteten Gitter aus Trag- 
fliigeln mit endlicher Spannweite und verschwindender Fliigeltiefe theoretisch noch kaum 
behandelt worden, da sie nur bei wenigen Anwendungsgebieten auftreten. Die wichtigste 
Verwirklichung liegt wohl in den Sonderfallen der waagerechten und senkrechten Gitter vor, 
die man sich aus Fliigeln zwischen parallelen Wanden mit Spalt durch Spiegelung entstanden 
denken kann. Trotz dieser geringen praktischen Bedeutung mégen die nachfolgenden Unter- 
suchungen der Vollstandigkeit halber erlaubt sein, da die verwendeten Untersuchungsmethoden 
ohne weiteres auch auf den Fall der Profilgitter iibertragen werden kénnen, wie demnachst in 
Teil II gezeigt werden soll. 

Wie bei endlich vielen Tragfliigeln! besteht auch bei Systemen von Fliigelgittern die 
Schwierigkeit, da die am Ort des Fliigels berechneten Abwindgeschwindigkeiten nur fiir 
spezielle Gitteranordnungen reell sind, so daB sich gewisse Probleme, wie etwa die ,,dritte 
Grundaufgabe“: der Tragfliigeltheorie, nur fiir diese besonderen Formen in einfacher Weise 
lésen lassen. Auf diesen Sachverhalt wird an den entsprechenden Stellen noch besonders hin- 
gewiesen werden. 

Der Inhalt der nachfolgenden Ausfithrungen gliedert sich folgendermaBen: Zunichst wird 
in Ziff. 2 ein einziges Fliigelgitter betrachtet (n = 1). Daim Gegensatz zu dem vielbearbeiteten 
Fall eines Einzelfliigels hieritber bis heute noch wenig vorliegt, wird diese einfachste Anordnung 
etwas ausfiihrlicher behandelt. Es werden zunachst Formeln angegeben, nach denen sich 
Zirkulations- und Abwindverteilung entwickeln lassen und die denen fitr den Hinzelfliigel voll- 
stindig entsprechen. Weiter ergibt sich eine Umkehrformel fiir die Zirkulationsverteilung und 
ein geschlossener Ausdruck fiir den Gesamtauftrieb eines aus dem Gitterverband heraus- 
segriffenen Fliigels. Als Beispiel wird schlieBlich in zwei Sonderfallen das Gitter mit geringstem 
induziertem Widerstand durchgerechnet. Systeme von endlich vielen solchen Fligelgittern 
folgen in Ziff. 3 (n> 1); in Ziff. 4 wird aus ihnen der Sonderfall ,,alternierender‘’ Tragfligel- 
gitter herausgegriffen. Der einfachste Fall eines einzigen solchen Gitters wird in Analogie zu 
Ziff. 2 ausfiihrlich dargestellt, auch hier 148t sich (bei noch einer zusatzlichen Annahme) eine 
Umkehrformel angeben. Die etwas umstindliche Berechnung eines vorkommenden Integrals 
ist schlieBlich im Anhang Ziff. 5 untergebracht. 

Lat man den Abstand zweier Nachbarfliigel im Gitterverband immer gréfer werden, so 
streben die gefundenen GesetzmaBigkeiten gegen die entsprechenden im Falle eines Einzel- 
fliigels bzw. endlich vieler solcher. Durch diesen Grenziibergang 1]4Bt sich somit an jeder Stelle 
der Rechnung der Zusammenhang mit bereits Bekanntem herstellen. Auf Fragen rein mathe- 


1K. Nickel, Ing.-Arch. 20 (1952), S. 363. 
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matischer Natur, wie z. B. die nach der Existenz der vorkommenden Integrale, der Vertausch- 
barkeit von Integration und Differentiation usw., wird an dieser Stelle nicht eingegangen!. 


2. Der Fall n=1, a) Ein Tragfliigel erstrecke sich in Spannweitenrichtung lings der 
y-Achse von a bis b (Abb. 1). Die értliche Zirkulation sei I"(y), die von dem Tragfliigel am 
Ort der tragenden Linie hervorgerufene Abwindgeschwindigkeit werde mit w(y) bezeichnet. 
Unter den Vernachlassigungen der Theorie der tragenden Linie hangen dann nach L. Prandtl 2 
Abwindgeschwindigkeit und Zirkulation nach folgender Gleichung zusammen: | 


b 
ehat “dI'(n) dy 
BO a Oa Fa (1) 


a 


Dabei ist an der singularen Stelle 7 = y das Integral als Cauchyscher Hauptwert aufzufassen, 
was durch das Symbol f gekennzeichnet werden soll. 


Betrachtet man statt des einen Tragfliigels der A 
Abb. 1 ein regelmaBiges Tragfliigelgitter nach Abb. 2, ° 
wobei die einzelnen Tragfliigel jeweils um den Vektor 
h =h,+ ih, gegeneinander verschoben sind 3, so ist serra menccmneeroeeN 
(1) zu ersetzen durch 

COs. 
tS ‘ (n) dn 
hy) Age S ¢ dyn y—yn 
Se rates hy 
eee a h 
pet ie (yn + vh) dij 
pemeee ¢ dy y—n—vh* 


Nun ist aber ]"(yn+-vh) =I'(n), 
da die Zirkulationen der ein- 
zelnen Gitterfliigel iiberein- 
stimmen, und weiter gilt die 
Partialbruchzerlegung 


Abb. 1. Einzeltragfliigel. Abb. 2. Tragfliigelgitter. 


co 


1 1 
ctg wz = — > . Damit findet man also 


~~ — Y 


=-— co 


b 
fara) 
w(y) =¢ Gq 8 (¥ — 1) an, 


a 


oder mit der Schreibweise * k == 
, i 
dI(1 
(9) = Fy ay te ROY) a. (2) 


a 


Fir k — 0 geht (2) in die Gleichung (1) iiber. 
b) Um bei bekannter Zirkulation I(y) den Abwind w(y) in einfacher Weise berechnen zu 
kénnen, formt man (2) noch durch partielle Integration um in 


b 
(9) = Fy ay Mn) ets My — 0) dr. (2a) 


1 Man vel. K. Nickel, Math. Z. 58 (1953), S. 49. 

2 L. Prandtl u. A. Betz, Vier Abhandlungen zur Hydrodynamik und Aerodynamik. Neudruck aus den 
Verhandlungen des III. Internationalen Mathematiker-Kongresses zu Heidelberg und aus den Nach- 
richten der Gesellschaft der Wissenschaften zu Gottingen, Gottingen 1927, 

3 Dabei muB fiir h, = 0 offenbar die Einschrankung h = h, >b —a getroffen werden, wenn sich die 
Einzelfliigel nicht zum Teil iiberdecken sollen. 

4 Die Einfiihrung der GréBe k ist fiir das folgende zwar nicht notwendig, laBt jedoch die Parallelen 
zwischen Hinzelfliigel und Fligelgitter etwas deutlicher hervortreten. 
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Zerlegt man weiter w(y) und /’(y) in einen symmetrischen und einen antimetrischen Anteil 
+b +6 
lcs +y) = 0 (5 y) 
a+b re) a+b 
pam e 
I'(y) =T'(y)-+T_(y) entsprechend, 


so findet man nach einigen Zwischenrechnungen } mit der Identitat 


w(y) = w,(y) + wy) mit 


sink (y—" ai ) cos k(y care x * + sin k(n —" 5) cos k( ae = ‘) 
ctg k(y — 7) = Paine a+b 
sint k(y — ) sine & (n — 
2 2 
den Wert 
b I'_(n) cos k (1 — =) dy 
ih 2 
wy) =n (2b) 
4a dy ( a+ ") ( a+ *) 
nk{y — — sin k {4 ——— 
; 2 2 
und 
aes | spr awie I. ()4 
ea a (1) dy 
0.0) =Fr 7 cos k(y : \¢ = rt (2c) 
| : sin k(y— D )=sine(n— 9 )| 
Nun gilt ? fir m= 0,1T>. .. 
o AS] 
il Va? — v2 om 2 ’ 
1% ¢ u—v aed =, qa, ar? uma (3) 
—4 
mit 
i oe hb Swoncen uey ——= 
CRS es Urea sama 2 ies Saw aS 2) fir A>1 
Setzt man 
: b 
w= sin k(y—*F ) v= sin k(n“), a =sin ke, 
so folgt aus (3) 
b sin k(b — 7) sin k(n — a) sin™k (n—" A. * cos k (—) 
d a: 
oe sink (y— 95) *) — sin k(n —*) 4 
2 
Eee) (3a) 
2 
=— > qo, sin? ps2 > sin™ +1—24 Bere els ab 
7=0 A 
Durch (3a) wird nun folgender Ansatz nahegelegt: 
P{n) = Vein EO — 9) sin BG — a) sine + (y— 9"), (4a) 
a ae ; agen a+b b\ 
lie (7) Ysin k(b — n) sin k(n — a) sin? Kn =a ) eos K(n — mae ‘ (4b) 
mitm = 0,1,... Diese Zirkulationsverteilungen ergeben dann nach (3a) die Abwindgeschwin- 


digkeiten 


) D, ea@A— 2m — 2) ints RY poate At yoo a ‘), (5a) 


w_(y) =~ cos ees - 


) 
a 


1Vegl. J. Dérr, Ing.-Arch, 19 (1951), S. 66. 


2 Vel. A. Betz, Beitrage zur Tragfligeltheori : PS ce Fh ; 
eckigen Fliigels, Diss. ee i010. geltheorie, mit besonderer Beriicksichtigung des einfachen recht- 
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kx 2 
Ne S32 Pm I) sinh aint tt bf SE) | 
Ex (5b) 
k — 
ae gz, (2A — 2m — 2) sin?+ k? 7 © sin2m+ 2-24 foes | 


Fir k > 0 gehen (4a) und (4b) itber in den Ansatz von Prandtl! und Betz2 fiir den ein- 
fachen Tragfliigel. Die vorstehenden Formeln lassen sich auch derart umformen, daB die Poten- 
zen von sin k/ y — * 3 g 


durch trigonometrische Funktionen der Vielfachen des Arguments 
a+b : : 
K(y— =) ersetzt werden, doch ergibt sich dadurch keine wesentliche Vereinfachung. 


Fiir manche Zwecke giinstiger diirfte analog zur Trefftzschen * Umformung der folgende 
Weg sein: Man setzt 


sin k(y—*7~) = sin 9S“ cos s, sin k(n “5 ) = sin BF cos (6) 


und schreibt w(y) = w(s), I'(n) = I(t). Damit gehen (2b) und (2c) itber in 


b—a 
Lee 2 n 

4 k V os Fathom bey ON ee 
w_(s) = — : (7a) 

; b—a sin s ds ¥’ cos s —cost 

4asink 
2 0 
und 
k i a | ae intd | 
a ie : ae ji ORs ', (t) sin t dt . (Tb) 
nhewet sins ds cos s— cost Nee hea 
4 sin 6 1—sin? k 5 cos? t 


Wegen der bekannten 4 Forme] 


7% 


aiaegearrs dt =cosms (m = 1,2,...) wy) 


a4 cos s — cost 


0 
ist dann folgende Setzung zweckmaBig (m — 0,1,...): 


[_(t)=sin 2 mt (9a) 
und 
P(t) =f 1 — sin? k° =" cost t sin (2 m + 1)t. (9b) 
Dazu gehéren die Abwindgeschwindigkeiten 
= sin 2m 
w_(s) = ——— yi = (tae ie 5 * cos? s saa 6 : (10a) 
4sink 
und 
w,(s) = eee (1 FO ere 
4 sin ree ae need 
5) (10b) 
ae sin be cos (2m-+1)coss. 


Mit den vorstehend abgeleiteten Formeln lassen sich manche aerodynamischen Probleme 
bei Fliigelgittern lésen, insbesondere kann man damit die Verfahren zur Lésung der ,,zweiten 
Grundaufgabe“ von einem Einzelfliigel auf ein Fliigelgitter tibertragen. In den beiden Sonder- 
fallen des waagrechten und senkrechten Gitters (kreell oder rein imaginar) werden samtliche 


1 Siehe FuBnote 2 von S. 109. 

2 Siehe FuBnote 2 von S. 110. 

3 E. Trefftz, Z. angew. Math. Mech. 1 (1921), S. 206. 

4 Vel. K. Nickel, Math. Z. 53 (1950), S. 28, Formel (14a). 
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CréBen rein reell, nur ersetzen sich im zweiten Fall die trigonometrischen durch die ent- 
sprechenden Hyperbelfunktionen. Dabei ist jedoch die richtige Normierung der auftretenden 
Quadratwurzeln zu beachten, die man z. B. durch den Grenziibergang k—> 0 feststellen 
kann. Fiir beliebige Gitter (k komplex) ist auch w(y) = u(y) + 1 w, (y) komplex. Darin 
ist w,(y) als Abwindgeschwindigkeit zu deuten, der zugehérige Imaginarteil w, (y) stellt 
eine dazu senkrechte Geschwindigkeit dar, die in der ttblichen Theorie der tragenden Linie 
vernachlassigt wird. 

c) Ist die Abwindgeschwindigkeit w(y) bekannt und die Zirkulation [(y) gesucht, so 
stellt (2) eine lineare Fredholmsche Integralgleichung erster Art fiir dI'(y)/dy dar. Ihre Lésung 


lautet 


dI(y) ik 1 r k a ‘| 
dy % sin k(b — 7) sin k(n — a) re : a 


(11) 


b 
w(y) 
~ 1g Ae sin k(h — y) sin k(y — a) dy) 


Dabei ist die Gesamtzirkulation I eine beliebige Konstante mit der Bedeutung 


b 

r= | oO) dy = T(b) — Pa). 
; af 

Da die Zirkulation I(y) an den Tragfliigelenden verschwinden soll, da also ["(a) = I'(b) = 0 

sein soll, findet man J’ = 0. Damit ergibt sich durch Integration von (11) 


" b 
i l w(y) Sie a ee 
— Seer ree k(b— k(y — a) d : 
m sin k(b — t) sin k(t — a) (¢ sin k(y — t) Vsin k( y) sin k(y — a) >) dt 


a ‘a / 


My) == 


Hierin 1a8t sich die Integrationsreihenfolge vertauschen und das innere Integral geschlossen 
berechnen. Man findet so als Analogon zur Betzschen Umkehrformel 


. b — be 
sin? k 5 *cosk(y—* +) cos k(n —* 4") cost? 3 * sin k(y—° £2) in e(n—* 4) 
ce b—a b—« rr 
sink cos k sin k |y — | 
2 Z 
_ sin k(b — y) sin k(y — a) sin k(b — 7) sin k(n — a) d 
* baa b—a., bf 
sink 5 cos k 5 sin k|y — | 


b 
Durch nochmalige Integration ergibt sich der Gesamtauftrieb 2 A = i I(y) dy eines Fliigels 


aus dem Gitterverband. Auch hier 14Bt sich wieder die Integrationsreihenfolge umkehren und 
das innere Integral durch elementare Funktionen ausdriicken. Da diese Ausrechnung jedoch 
sehr mithsam ist, soll hier ein anderer Weg beschritten werden. Man faBt dazu A als induzierten 
Widerstand auf, den die Zirkulationsverteilung I'(y) auf eine zweite Verteilung "*(y) ausiibt. 
Dabei sei /*(y) so bestimmt, da die zugehérige Abwindgeschwindigkeit w*(y) = list. Nach 
einem bekannten Satz der Tragfliigeltheorie ? ist dann 


b b 
A= [T(y) 1+ dy = | w(y) P*(y) dy! (13) 
1 Siehe FuBnote 1 von S. 109. 
2 Ein konstanter Faktor, Luftdichte und Anblasgeschwindigkeit enthaltend, ist hlassi 
° Vgl. die in FuGBnote 2 von S. 109 zitierte epee S. 38. See eet kay ae 
_ 4 Diese Uberlegung ist offenbar nicht auf das betrachtete spezielle Fliigelgitter beschrinkt, sondern 
gilt ganz allgemein fiir beliebige Zusammenstellungen von Tragfliigeln. Weiter lassen sich auf dieselbe 
Weise auch Formeln fiir Rollmoment, Langsmoment bei Pfeilfliigeln usw. gewinnen. 


XXII. Band 1954. Nickel: Uber spezielle Systeme von Tragfligelgittern I. 113 


Es bleibt jetzt nur noch [’*(y) zu bestimmen. Mit der Formel 


b 
k /sin k(b — y) sin k(y — a) ite 
T ¢ sin k(y —n) dy = sin X( —*f4) (14) 
folgt aus (11) 
P ; a+b 
dI'*(y) eee sin k (1— 5} 
dy /sin k(b — 7) sin k( —a) 


und daraus durch Integration 


, cosk (v +) + /sin k(b — n) sin k(n — a) 
(PS = —] na 
(n) =+log , (15) 


Setzt man (15) in (13) ein, so kommt schlieBlich die gesuchte Auftriebsformel 


> as cos k (y—* 4) + Voi R(— ssi FU =a) 
A= Py) dy =>- w(y) log dy. (16) 


b—a 


2 


a e cos k 


Fiir w(y) = 114Bt sich das Integral auswerten (vgl. Ziff. 5) und wird zu 


b cox k(y—* 5) + sin k(b — y) sin k(y —a) 
- | log 


es cos k 


Oy = — “Flog cos ke > : (16a) 


2 


Pls 


b—a 


Die vorstehend abgeleiteten Formeln sind zwar allgemein giiltig, doch ergibt sich fiir 
beliebige komplexe Gitterkonstante k die Schwierigkeit, daB® w(y) bei reellem I'(y) komplex 
wird. Im allgemeinen wird nun von w(y) = w,(y) + i w,(y) nur der Realteil w,(y) gegeben 
sein, wahrend ja w,(y) fiir gewéhnlich vernachlassigt wird. Setzt man also w(y) = w,(y) z. B. 
in die Forme] (12) ein, so ergibt sich eine komplexe Funktion [(y), wahrend ihrer physika- 
lischen Bedeutung nach die Zirkulation eine reelle Gré®e ist. Um ein reelles '(y) zu bekom- 
men, miifte man zu w,(y) noch eine imaginare Funktion i w,(y) addieren, deren Bestimmung 
jedoch auf eine bis heute noch nicht geléste Integralgleichung fiihrt. Damit lassen sich die 
Formeln (12), (15) und (16) im allgemeinen nur fiir reelles oder rein imaginares k (Gitter waag- 
recht bzw. senkrecht) benutzen; es sei denn, daB mit w,(y) auch zufallig w,(y) bekannt ist. 

d) Ist die durch (15) dargestellte Zirkulationsverteilung ["*(y) reell, so ist sie nach einem 
Ergebnis von Munk! gerade diejenige, die bei gegebenem Auftrieb A geringsten induzierten 
Widerstand? 


W = [ (y) w(y) ay (17) 


besitzt. Sie soll fiir die beiden Grenzfalle des waagerechten und senkrechten Gitters noch etwas 
naher betrachtet werden: 
&) h==h,. Man setzt a=>—t, b=}, also b—a=1 und h,=I1+s, also 
= =Ta3 mit der Spaltbreite s (Abb. 3). Damit findet man aus (15) und (16a) die Zirku- 
lation mit konstanter Abwindgeschwindigkeit und dem Auftrieb A 


A 1 get : Dahl 
I'*(n) = log ah = : as (18a) 
ie, NOE POP 5) COE ar a) 


1M. Munk, Isoperimetrische Aufgaben aus der Theorie des Fluges, Diss. Gottingen 1919. 
2 Vel. FuBnote 2 von 5S. 112. 
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Der zugehdrige induzierte Widerstand berechnet sich nach (17) und (16a) zu (W—W,, fir 


$—> 0) 


Wit us (19a) 


ie 8(1 + s)? log cos sD 

Diese Ergebnisse stimmen (bis auf einen Faktor 2 bei [’*) tberein mit den Formeln, die von 
A, Betz! angegeben wurden. In Abb. 3 sind einige Zirkulationsverteilungen nach (18a) ein- 
getragen, fiir s > oo ergibt sich die bekannte elliptische Auftriebsverteilung. Der induzierte 
Widerstand wird durch die Gitteranordnung verringert und fallt von W,, auf Null ab, wenn 
s— 0 geht (Abb. 3). 

B) h=ih,. Schreibt man hier hy =s und setzt wieder a =— 3}, b =} (vgl. Abb.4), so 
lauten die entsprechenden Formeln 


x oj —n 
1 Best G7) Vem arc cos —————— (18b) 
" rt It 
Ss log Co} os Cp} ow 3 


und 


oo = fe (19b) 
= 8 s? log Co} 3 


Fallt s von groBen Werten auf Null ab, so geht die Zirkulationsverteilung wieder von der Halb- 
ellipse in eine Rechtecksverteilung iiber (Abb. 4), nur beginnen merkliche Anderungen diesmal 
schon erheblich frither als im Falle «). Wahrend dort I(y) fiir s = 0,1 noch recht nahe bei der 
Halbellipse s = oo verlauft, ist hier die Zirkulation bei s = 0,1 schon iiber einen groBen Teil 
des Fliigels praktisch konstant, und der Auftriebsabfall konzentriert sich auf einen kleinen 
Bereich am Fliigelende. Weiter ist jetzt der induzierte Widerstand bei konstantem Auftrieb 
stets héher als fiir den Einzelfliigel und wachst fiir s > 0 tiber alle Grenzen. In Abb. 4 ist daher 
noch W,,/W eingetragen. Diese Kurve kann gedeutet werden als induzierter Widerstand des 
betrachteten Fliigelgitters bei festgehaltener Abwindgeschwindigkeit w. 


3. Der Fall n >1. a) Betrachtet man n Tragfligel mit derselben Gitterkonstanten k, so 
ist Gleichung (2) zu ersetzen durch 


b, 
kof aly(n 
(Yu) Be ace K(¥,— My) an, (wu = 1,.- +5) (20) 
oder 
b, 

kd ~ 
wz, (Yu) Se ear, ¢ I’, (ny) ctg k(¥,—) dn, (uw =1,...,n) (20a) 

— 


v 


wenn sich ein beliebiger Fliigel des y-ten Gitters von a, bis b, erstreckt (Abb. 5). Welchen 
Fliigel man dazu herausgreift, ist offenbar gleichgiiltig. 


In diesem allgemeinen Falle ergeben sich wieder fiir reelles und imaginares k reelle Abwind- 
geschwindigkeiten w,(¥,), falls alle a, und b, auf einer Geraden liegen, doch sind dies offenbar 
nicht die einzigen derartigen Falle. Z.B. liefern auch bei der nebenstehend skizzierten Anord- 
nung (Abb. 6) alle (reellen) Zirkulationen J’, (n,) rein reelle Abwindgeschwindigkeiten w, (y,). 
Weitere solche Kombinationen lassen sich leicht angeben. Fir die nachfolgend abgeleiteten 
Formeln gelten daher die entsprechenden Uberlegungen wie im Falle n =1, nur sind jetzt 
noch mehr spezielle Gitterformen vorhanden, bei denen die dort genannten Einschrankungen 
wegfallen. 
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b) Zur Berechnung der Abwindgeschwindigkeiten w, (y,) setzt man entsprechend zu 2b, 


L(4) = (— y)/- 


(b, —n,) sin k(a, — sex 


1 


IT sink 


a 


n 


(21) 


| 


n ungerade 


=n aes 


© ‘y WoJUBISMOyIE4IL OWOII]S yrur usszSpaSnpy3vry, u woa maysig “Gg ‘qqy 


WaT]IeI UleT yor uIEq3TS[eSn{FSery, = 


T9MZ sne urazsig -9 "Gqy { 


‘udIUIT uspueseiy Jap a 
UQ We 19z10YSIpurmygoseSpurmqy 
+ 
sS 


‘JOY SIPUIMYOSISpUIMGy 193UeISUOY 1oq puKysIeprAy 
eWIIZNpur pun WeSUNIISWIIASqarayjny “109713 aSn pA sazqov1yUeg “Pp ‘qqy 


we me . wv 0 
Te 


. Lan i T 


s 


“‘IOYSIPUIMOSISpuIMqy I9}ueJsuoy 19q pueysIepr jy 
JOIVIZNPUL pun UeISuNeWIAsqoeyny ‘A9WS[SnTy soyqooisee ny “e ‘qqy 


Z L 


— 


-1 Sin k ny | 
m 
kn, lyo, coskn, + yox41sinky,] ngerade 


0 


& sin?* kn, [y2 4 cos kn, + vex- 


m 
y, + » sin?*—1 
‘ “u=l1 
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mit Konstanten y,, die nicht von der Tragfliigelnummer » abhangen. Mit Formeln entsprechend 
zu (3a) 1 findet man dann fiir alle n = 1 Sc 
1 


—|+m 


S) sint’*—l hy, (A, cos k y, + Azy41 sin k y,) (22) 


1 


2 
Ww, (¥,) =r A, ae 
mit Konstanten A, die noch von a,, b, und k abhangen. Sie lassen sichi. a. nicht durch elemen- 
tare Funktionen geschlossen ausdriicken. ; 

Die Formeln (21) und (22) entsprechen vollkommen denjenigen bei endlichen Tragfliigel- 
systemen”; fiir der Grenziibergang k > 0 ist nur zu beachten, daB dabei der hichste Koeffi- 
zient von (22) verschwindet. 

c) Bei bekannten Abwindgeschwindigkeiten w,(y,) ist (20) ein Integralgleichungssystem 
fiir dI’,(n,)/dy, . Seine Lésung ist unter sehr allgemeinen Voraussetzungen ¢ 


dI,(m) 4k (=a P 
i i = n(») 
j= IT sin k (by, — ny) sin k( ax — nv) 
Pr (23) 
_ Z _————— n ungerade| Lae : 
ee ng 1, (.) 88 (Fu — MP) — HT sin b(b, — y,)sin k(% — Yu) Ay 
sl ay, loee k(y,— "») nm gerade | is 
mit 
P,(m) =1B, cos k(n _ SS sin®—1 k 7 
=1 
n—3 ; 
2 
>» sin?” k n[ B.,, coskn + Bo,41 sin ky] (n ungerade) , 
a5 eee n-—2 


2 
2, + >) sin?*—lkn[B,, cosky + Bo,41sin kn] (n gerade) 
“=1 


und n beliebigen Konstanten B,. 

Wahrend sich im Falle n = 1 durch Integration geschlossene Formeln fiir /‘(y) und A 
angeben lieBen, kommt man jetzt bei der entsprechenden Umformung auf hyperelliptische 
, Integrale, die sich i. a. nicht durch elementare Funk- 

tionen geschlossen ausdriicken lassen. 
TTT DS Wie dies bei Systemen von Tragfliigeln durchgefiithrt 
wurde‘, 1aBt sich auch hier beweisen, da8B es zu beliebigen 
Abwindgeschwindigkeiten w,(y,) immer genau einen zu- 
gehoérigen ,,Satz‘* von Zirkulationsverteilungen I’, (7,) 
£ 5 gibt, der das Gleichungssystem (20) befriedigt. Man 
aN eal 4g erhalt die I’,(7,) aus (23); dabei sind die n unbekann- 


i b@)| [J ten Koeffizienten B, in P, eindeutig durch die Forde- 
A rung bestimmt, das die Zirkulationen J’,(y,) an den 
Ww Fliigelenden verschwinden sollen, daB also 
(Tae DS vA Day Sy 
ay by 


sein soll. Zum Beweis ist genau wie dort das Nicht- 
verschwinden einer Determinante zu zeigen. Da dieser 
Nachweis jedoch sehr langwierig ist, soll er hier nicht 
gefiihrt werden. 


NG ERY 4, Alternierende Fliigelgitter. a) Ein spezielles Sy- 


stem von Fliigelgittern, tiber das sich im einfachsten 
Fall noch weitere Aussagen machen lassen, erhalt man 


durch folgende Anordnung: 


Abb. 7, Alternierendes Tragfliigelgitter. 


1 Diese Formeln sind in der in FuBnote 1 von §,109 zitierten Arbeit angegeben’oder lassen sich durch 
einfache Umformungen daraus gewinnen. 2 Vgl. die in FuBnote 1 von S. 108 zitierte Arbeit, S. 370. 
3 Vgl. diein FuBnotel von $.109 zitierte Arbeit. 4 Vgl. diein FuBnote 1 von S. 108 zitierte Arbeit S. 366. 
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Zu jedem Gitter mit der Zirkulationsverteilung [’,(y,) im Intervall a, <y,<b, gebe es 
ein zugehériges, dessen Bezugsfliigel sich von a, — =a, — h/2 =a, —a/2k bis b, = b, — h/2 
=b,—a/2k erstrecken mége, und das auf dieser Strecke die Zirkulationen T, (7) 
=—TY, (n, — zi) besitzen soll (vgl. Abb. 7). Das zweite Gitter ,,halbiert“* also das erste 
und ist foe nentgegengesetzt“’, die Vereinigung dieser beiden soll als ,,alternierend“ be- 
zeichnet werden. 

Damit lautet (20) 


aL'y (1) dl (nv VARs 
ae -[s ~ dy 88 BY — Me) ana 3 gg He) ety (yy, — Tn) diy 


4 


J 


2 
sin 2 k(y,,— ny) 


und, mit ctg k(y, — ,) — etg k(y, — 4,) = etg k(y, —,) + tg k(y, —n,) = 
umgeformt, 


es dIy (ny) dn, Ly (1%) 
Wy (Yu) = Ant > ¢ dn, sin k’ (yu — mp) =; It a 7 ¢ sin k’ (yp, — my) dy» 


wenn man h’ = h/2, also k’ = 2 k setzt (Abb. 7). 
b) Durch einen Ansatz entsprechend zu (21) erhalt man hier die Beziehungen! (m= 0, 1,...) 


et.) =(— 1) \—it sin ki (b, — ,) sin k’ (a, — ,) x 


Yi + & sin?*—1 k’ n, [ye, cos k’ n, + yous. sin k’ ny] nungerade ; (25) 
x . #=1 
& sin? k’ n, [y., cos k’ n, + yoy 41 sin k’ n,] n gerade 
x=0 
und 
[a] +m 
W4(Y 4) = > sin?” k’ ¥u [Con cos k’ vy, + Co, 44 sin k’ vil - (26) 
“x=0 


c) Gleichung (24) hat bei bekannten w, (y,,) die Auflésung? 


dIy(yv) 4k’ (—1)’ 
dny It 


; 
Bie ie iyee TAG ian On (te) 


“=! L 


(27) 


ctg k’(y, —n,) n ungerade n 
| = | V — TIT sin k’(b,— y,,) sin k’(a,— y,) dy, 
eres n gerade A=1 
sin k’ (yp — nv) 


by 
+>) (= G wl) 


| 
mit 
n—1 
D,+ » sin?*—1 k’ y[D,, cos k’ n+ Dz,41s8in k’ ny] (nungerade) , 
Qn (7) =jn—2 *~ 
Ss sin?” k’ yn [D,, cos k’n + Dz,,4 1 sin k’ y] (n gerade) 


“u=0 


und n beliebigen Konstanten D,. Esla8t sich auch hier zeigen, daB man die Abwindgeschwin- 
digkeiten w,(y,) beliebig vorschreiben darf und dann durch (27) und die Forderung I’, (a,) 
=I", (b,)= 0 eindeutig bestimmte Zirkulationen I’, (7,) erhalt, die das Gleichungssystem (24) 
erfillen. 


1 Vgl. FuBnote 1 von S. 116. 
2 Vel. FuBnote 1 von S. 109. 
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d) Im einfachsten Falle n = 1 kann man noch weitere Ergebnisse erhalten. Gleichung (24) 
vereinfacht sich dann zu 


b 
k’ ( dI(n) dy 24 
wi =F5¢ dn sin k’(y —n) - Eid KO es rte — (ct) 


a 


was man mit der Identitat 


b b 
sin (y—2 5) cos #( — 25") 4 sink (n— 95 ) cos w (y—*F 


1 2 
UA het sin? K’ (y—" : *) — sin h’ (nF ' 
wieder aufspalten kann in 
b 
k’ a| , a+b T_(n) dy 
ce ead ea (y ae aa rab } (24b) 
a sin (y— 5 )—sin we ( 5 J 
und 
k d b cos (y= " 2) ay 
w,.(¥) aig - (7) I ca (24) 
a ar 


Setzt man hierin 


Ly) = Vsin k’(b — n) sin k’(n — a) cos K(n—* - * sin?™+1 ae - 5 (28a) 


und 
L(y) = Ysin k’ (6 k Poy oan 28b 
=a sin — 7) sin k’(y — a) sin 1 Se (28b) 
so ergeben sich durch Rechnungen, die wie in 2 b) verlaufen, die Formeln 
7 m } 
wiy)=> D)i(2A—2m—2) sind? Fo aint sits e(g 84) — 
(29a) 
k’ m+1 aa ‘ 
83h, a1—am— ants eh gt antes toon (gtd 
und 
- 
w,(y) == cos (x9) > aa (24 — 2 m— 1) sin?4 ped Deas o in2m—2A fy (nee t >). (29) 


Die Trefftzsche Umformung (6) ba 


w_(s) = As Sa rg ass yin a | (24d) 


mae sin $s cos s — cost bE=6 
/1—sin Lae cos? t 


Pe pes 


und 
‘ yi —sin? k’ ome cos? s bs 

AO k “ I’, (t) sint 24 

+ eee i leet sin s cos s — Cost oo) 

0 
Damit entsprechen sich die folgenden Zirkulationen (m = 0,1,...): 
: b — 
Et) =/i — sin? k’ “ cos?t sin 2 mt (30a) 


I’, (t) =sin(2m-+ 1)t (30b) 
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und die Abwindgeschwindigkeiten 


w_(s) = ee (1 Sy eee ae gs 
esi peas 2 sin s 
sin l ) 
3la 
ere po 4 9 ( 
—7 ain 5 Cos2mscoss 
Pe ees oe a St i icaette ee the Ls (31b) 
pies OL 2 sin s 
4 sin k’ ——— 
2 
Die Auflésung von (24a) bei bekanntem w(y) lautet ! 
ar ki ; 
= = |b f Py) sin by —2 +) ay 4 
dy % |/sin k’(b — 7) sin k’ (4 — a) a 2 
(27a) 


b 
-- ag w(y) etg k’(y — n) sin k’(b — y) sin k’(y — a) ay| : 


Ist I'(y) symmetrisch zu“ sa 


, also (ye ary) mit | hi Cae Shim 9) so 


b 
verschwindet dle (y) sin k’ [¥— 2 = * dy. In diesem Falle 14Bt sich (27a) noch einmal inte- 


grieren und gibt 


b sees Sey qs! a 7 ee 
4 2 2 2 
L(y) == | w, (y)log ; Se 
2sin ki -——* cos ¥ (y + b) sin fy — 9 
=" n — a — b) sin — |y— 
: 2 2 : 2 0 (32) 
__ Vsin k’(b — y) sin k’(y — a) sin k’(6 — 7) sin k(n —a) | dy. 
A eye cm b) in 24 
2s] 5) By Waa ae sp cers ce] 


Fir k’ + 0 geht (32) wieder in die Betzsche Umkehrformel iiber. 
a+b 
2 
Formel fiir den Auftrieb, entsprechend (16) beim einfachen Gitter, 1aBt sich hier nicht mit 
elementaren Funktionen angeben, da die Berechnung der Auftriebsverteilung [*(y) mit 

w*(y) = 1 ebenfalls auf elliptische Integrale fithrt. 
Wie in Ziff. 2 sind auch hier die Falle des waagerechten und des senkrechten Gitters die 
einzigen, in denen w(y) bei reellem [’(y) reell wird. 


Ist [‘(y) unsymmetrisch zu 


, so addiert sich zu (32) noch ein elliptisches Integral. Eine 


5. Anhang. Zur Berechnung des Integrals (16a) 
b cosk(y—2 +4) 4 yank ys My =a) 
4 2 
T=7 [log eae dy =4 


cos k 


Lysine (y—244) 
d 


/sin k(b — y) sin k(y — a) 


a 


act 


‘cos k (y— 2 +) + Vain k(b — y) sin k(y — a) 


b 
=_ [le Gah 
: cos k (y— 5} ) EH w 


dy 


betrachtet man, wie es ahnlich schon A. Betz? getan hat, den allgemeineren Ausdruck 


5 b ccosk(y—*2) 4 yom e(b— yen Fy —a) F 
(0) = 7 [les Nf 
b 
a ccos k(y— °F ) ie kO ym EO 


1 Vel. FuBnote 1 von S. 109. 
2 Vel. FuBnote 1 von S. 114. 
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mit ¢> 1, aus dem man riickwarts I =lim I (e) gewinnen kann. Durch Differentiation nach 


e>l 
¢ unter dem Integralzeichen findet man 
3 b—a : a+b a+ ) 
b : 2 a 2 pans SE UAE =—— 
dl (e) ae sin k 5 sin? k (y D ) eos k (y 9 : 
de k a Oe 


und dieses Integral 148t sich leicht ausrechnen, wenn man sin ky — = Z ‘ als neue Integra- 
is 


tionsverinderliche einfiihrt. Man findet dann mit « = sin k 5 £ 
dI(e) _ 4 Tee 1 — a? 
de k? (1 —e?) e? — a? 


Wie man durch Ableiten bestatigt, gehért dazu das Integral 


l+e C 
Ve2 — a? +)/1 — a? ; 
LaBt man ¢ > co gehen, so ergibt sich C = 0. Also gilt schlieBlich 


4 


4 
I(c) = 43 log 


: Al TGR l+e 4m 1 
IT =lim I(¢) = — lim lo = ——lo 
Ai (©) Be ey e Vez — a? + 1 — a? ke 2 V1 — a 


(Eingegangen am 6. Juli 1953.) 


Anschrift des Verfassers: Dr. Karl Nickel, Instituto Aerotecnico, Cérdoba, Argentinien. 
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Einige Anwendungen eines nichtlinearen Elastizititsgesetzes*. 
Von F, Jindra, 


1. Einleitung. Die klassische Elastizitatstheorie stiitzt sich durchaus auf die Giltigkeit des 
Hookeschen Gesetzes, das einen linearen Zusammenhang zwischen den Spannungen und Ver- 
zerrungen annimmt. Dabei diirfen die Verzerrungen als unendlich klein behandelt werden. 
Viele Werkstoffe zeigen aber schon bei diesen kleinen Verzerrungen ein nichtlineares Verhalten, 
so daf das Bediirfnis nach einem nichtlinearen Gesetz bestand. H. Kauderer! hat daher ein 
allgemeines nichtlineares Elastizitatsgesetz fiir kleine Verzerrungen angegeben, das simtliche 
elastischen Probleme umfaBt und die klassische Theorie als Grenzfall enthalt. 

Die vorliegende Arbeit beschaftigt sich mit den Anwendungen dieses nichtlinearen Elasti- 
zitatsgesetzes auf einige elastomechanische Probleme. In einem einleitenden Abschnitt 
(Ziff. 2) werden die Grundgleichungen der nichtlinearen Elastizitatstheorie der kleinen Ver- 
zerrungen zusammengestellt. Im weiteren verweisen wir dabei auf die Arbeit von H. Kauderer. 
Wir behandeln zuerst den ebenen Spannungszustand durch Einfiithrung der Airyschen Span- 
nungsfunktion, wobei wir einige einfache exakte Lésungen angeben und dann die Biegung des 
Rechteckstreifens durch eine Einzelkraft mit der Methode der Stérungsrechnung behandeln 
(Ziff. 3). Im Anschlu8 daran stellen wir die Differentialgleichung der Airyschen Spannungs- 
funktion in Polarkoordinaten auf, teilen einige exakte Lésungen mit und behandeln durch 
Reihenentwicklungen als weitere Beispiele die Kreisringscheibe unter AuBendruck und den 
ebenen Keil mit einem Moment an der Spitze (Ziff. 4). Sodann folgen die rotierenden Scheiben, 
wobei wir zuerst die allgemeinen Gleichungen fiir die Spannungen herleiten und dann die 
Scheibe gleicher Festigkeit exakt, die Scheibe gleicher Dicke aber wieder mit Reihenentwick- 
lungen untersuchen (Ziff. 5). SchlieBlich stellen wir noch die Grundgleichungen der Torsion 
des Umdrehungskérpers auf und behandeln mit der Methode des Minimums der Fehlerquadrate 
als Beispiel den Kreiskegel (Ziff. 6). 


2. Das nichtlineare Elastizitatsgesetz fiir kleine Verzerrungen. In kartesischen Koordinaten 
x, y, zist der Spannungszustand vollstandig durch die drei Normalspannungen o,, o,, 0, und 
die drei Schubspannungen t,,, ty,, ts gegeben. Sind X, Y, Z die Komponenten der auf die 
Volumeinheit wirkenden Krdafte, so gelten die Gleichgewichtsbedingungen 


00x Oty x OTz x 
eee, ae 


oy Oz 
OTx y doy OT y \ se al 
Ox q dy is dz Bice (2.1) 
OT xs OTyz \ 00z 


ae dy | Oz Sag eae 

Unter Beschrankung auf unendlich kleine Verzerrungen, d. h. Verzerrungen, bei denen GréBen 
von der zweiten Ordnung in den Verschiebungskomponenten wu, v, w und deren Ableitungen 
gegeniiber den linearen Groen vernachlassigt werden kénnen, wird der Verzerrungszustand 
durch die drei Dehnungen 

meen ae ae. (2.2) 


oad OK = Bye oe a 


und die drei Winkelanderungen 


du dv __ ov Ow _ dw | ou 23 

Eat ree Fe Ge a ay aia Pal (2.3) 

beschrieben. Das Elastizitatsgesetz gibt den Zusammenhang zwischen dem Spannungs- und 
dem Verzerrungszustand. Die Spannungs-Dehnungsbeziehungen des nichtlinearen Elastizitiats- 


* Auszug aus der Stuttgarter Dissertation (Berichter: Prof. Grammel, Mitberichter: Prof. Siebel). 
1 H. Kauderer, Ing.-Arch, 17 (1949) 5. 450. 
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gesetzes fiir kleine Verzerrungen lauten in der allgemeinen Form 


0, = 3 K x(é) & + 2 6 (ys) (& — £0) » | 
dy = 3K x(é ) & + 2 G y (yo) (&y — €0) + 
5, = 3K (69) & + 26 y(ys) (& — &) 
Cey= Cy (Y2) Yay > Tyz= EY (YE) yes Tex = © Y (YO) Yeo - 
Dabei bedeuten K und G den Kompressionsmodul und den Schubmodul beim Hookeschen Ge- 
setz, €), und 5 die Dehnungsgréfen 


oy = (ee + ey + &)> | 


(2.4) 


(2.5) 


ye = 2 Olle — a)? Fe — 00? +(e a0) + hy + het ve | 


welche die Invarianten der Tensoren der Volumdehnung und der Gehaltsanderung darstellen. 
Die Dehnungsfunktion x(é) und die Scherungsfunktion y(yj) werden als 
Potenzreihen der Dehnungsgréfen (2.5) angesetzt : 


2% (€) =1 + e+ eet x%segt ''> | (2.6) 
y(y2) =1+ revit revit roy to-°-) 
Die Gleichungen (2.4) lassen sich umkehren. Man erhalt 


I 1 
b= 3 RMS) 0 + 7G Bl) (% — %) 


1 i 
by = 35g (50) Oo + aq B40) (oy — 9p) » (2.7) 
1 J : 
E=3K k (So) % +36 8 (t) (02 — %) » 
1 2 1 2 — 1 o Fe 
Pry = & B(to) Try » Pyz=G Bt) Tye» Pex = G B(to) Tex » 
wo die dimensionslosen Spannungsgréfen sy und t, durch die invarianten Beziehungen 
o 1 
eae Oo=—z (Gx + Oy + 02); | (2.8) 
2a 
p= > %= 3 r [(G, — 05)? + (oy — Gp)? + (6, — 9)"] + tay 4 Tye ty 


gegeben sind. Die Kompressionsfunktion k(s) und die Schubfunktion g(i) 
nehmen als Potenzreihen der SpannungsgréBen (2.8) die Gestalt an 


k (sy) = 1+ ky sy + hy 52 ++ kg s8 +-°--, | 

g(§) =1l+He+astseot'- } 
Mit Hilfe der drei Gleichgewichtsbedingungen (2.1), der sechs VerzerrungsgréBen (2.2), (2.3) 
und der sechs Spannungs-Dehnungsbeziehungen (2.4) oder (2.7) ist es méglich, die inneren 
Krafte und die Formanderungen eines dem nichtlinearen Elastizitatsgesetz gehorchenden, in 
bestimmter Weise belasteten und gestiitzten Koérpers zu berechnen. Die Aufgabe der Ermitt- 
lung des wahren Spannungsverlaufs bei willkiirlich vorgegebenen Kompressions- und Schub- 
funktionen k(s,) und g(t?) bringt uniiberwindliche mathematische Schwierigkeiten mit sich. 
Man vereinfacht die Aufgabe, indem man sich auf bestimmte, einfache Formen dieser Funk- 
tionen beschrankt. Es wird dann angenommen, dafB die Funktionen k(sy) und g (#2) in der Form 


K(s) =1, s)=1lt+ahq - (2.10) 


darstellbar sind und g(t) nur wenig vom Wert Eins abweicht. Diese Annahme kann bei vielen 
Stoffen als erfiillt angesehen werden }. 


(2.9) 


3. Der ebene Spannungszustand in kartesischen Koordinaten. Wir behandeln einen platten- 
formigen Kérper, bei dem die iiberall gleiche Dicke gegen die Abmessungen der -Grundflache 
gering ist. Die Scheibe wird nur an ihrem Rand durch Krafte belastet, die in der Ebene der 


1 Vel. H. Kauderer, a. a. O. S. 468. 
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Scheibe wirken und iiber die Dicke der Scheibe gleichmaBig verteilt sind. Wir sprechen dann 
von einem ebenen Spannungszustand !, 


Wir legen ein rechtwinkliges (x, y)-Koordinatensystem in die Mittelebene der Scheibe. Fiir 
den ebenen Spannungszustand gilt also 
O; =x =TMy = 0. (3.1) 
Bei fehlenden Volumkraften (X = Y=Z = 0) und mit (3.1) erhalten wir aus (2.1) die Gleich- 
gewichtsbedingungen des ebenen Spannungszustandes 


00x Oty x — OTx y doy 

ch ew Pea Pa er ogi (3.2) 
Wegen ty, =T,x = 0 sind auch die Winkel anderungen py. = y.x = 0. Die VerzerrungsgréBen 
(2.2) und (2.3) reduzieren sich dann auf 


_ Ou Lop dw du , dv 
saad A aaa mee PO Pay ae ae (3.3) 


1 1 
a ae k (9) 0 + 26 & (th) (Oz — oO), 


i 
oy aie k(s9) 09 + oval 8 (#9) (Gy — 99) » (3.4) 
1 
Pry = & 8(t) try | 
mit den Abkiirzungen 
1 2 2 
$0 = 9K (Ox + Fy) » Neeayar (Gp. 0, —=07 07", © Tey) (3.5) 
bestehen. Aus (3.3) leiten wir die Vertraglichkeitsbedingung des Verzerrungszustandes 
2 2 2. 
Or ex I Oey _ Pps y (3.6) 


dy? | Qx2 Ax dy 
ab. In diese Gleichung fithren wir die Ausdriicke (3.4) fiir die Verzerrungskomponenten ein und 
erhalten 
a2 


A lax Moo) — ag 8(4) 06} +3 gys (BC) On] + 5 pun (BU) 4] = gy BUH) Fey] 


2 
+53 den Laplaceschen Operator bedeutet. Die rechte Seite der Gleichung 14Bt 


2 
Pe 


7 Os" 
sich mit Hilfe von (3.2) umschreiben. Es ergibt sich 
i_Z fia | 1 0°g (16) 0°g (t5) 0g (t5)] __ 
a aK *(%) +36 & (1) (ox + oy) aC. Ox Ox + Oy ay? + 2 Try ax dy|~ °° (3.7) 
Da die Gleichungen (3.2) in ihrer Form mit den Gleichungen der linearen Theorie iiberein- 
stimmen, so kénnen die Spannungen auch hier aus einer einzigen Funktion F(x, y), der sog. 
Airyschen Spannungsfunktion hergeleitet werden, welche mit den Spannungs- 
komponenten folgendermafen zusammenhiangt 
or Sy See 
Ox Oy ~ 


(3.8) 


pee Sina 


Mit diesem Ansatz werden die beiden Gleichungen (3.2) identisch befriedigt, wahrend (3.7) 
iibergeht in 


al I 1 [OF 0g(t OF g(t? OF dg (t? 
Alog Ms) +36 88) Fl —se |e oat il ae OO 


| Oy? Ox Ox? dy? Ox Oy Ox dy 
mit 
es. AS af @F\t . OF BF 
ae ae =7e|(4F) ts 3(30 a5) 8 Gat aye (3.10) 


Dies ist eine quasilineare partielle Differentialgleichung vierter Ordnung, und die allgemeinste 


1 Siehe z. B. C. B. Biezeno und R. Grammel, Technische Dynamik, 2. Aufl., Bd. I, $.111, Berlin 1953, 
gx 
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Aufgabe des ebenen Spannungszustandes lauft auf die Lésung dieser Gleichung bei vorge- 
gebenen Randbedingungen hinaus. 

a) Einfache Lésungen fir die Rechteckscheibe. Wir betrachten zu- 
erst einige einfache Spannungszustande, die eine exakte Lésung der Differential gleichung (3.9) 
zulassen. 

Der in Abb. 1 dargestellte rechteckige Streifen von beliebiger Lange und gleicher Dicke, die 
gleich Kins gesetzt ist, wird langs der Querseiten durch gleichmaBig verteilte Lasten p auf Zug 
beansprucht. Die zugehérige Spannungsfunktion 


F= = y? (3.11) 
liefert nach (3.8) die Spannungskomponenten 
eich ee ~~ Ons er ee 3.12) ~ 
Ge Srp athe Oy =e a aay GSe 
und nach (3.10) 
1 2, 
So one te 0 =oaP- (3.13) 
¥ 
ha fu 
= aa 
Abb.1. Rechteckscheibe unter Zug. Abb. 2. Uberlagerung zweier Zugbelastungen. 


Wie man leicht sieht, erfiillt (3.11) mit (3.13) sowohl die Differentialgleichung (3.9) als auch die 
Randbedingungen. 

Durch Uberlagerung zweier einachsiger Spannungszustande in zwei zueinander senkrechten 
Richtungen erhalt man 


J 2 

ee (3.14) 

woraus 
O, =P, ~<“Oy =, Txy = 9 (3.15) 

und 
1 

Abb. 3. Rechteckscheibe unter Schub. $9 = OK (p + q) ° i = ae (p: aa ¢” ——— q) (3.16) 
folgt. Die Funktion (3.14) mit (3.16) geniigt somit der Gleichung (3.9), die Randbelastung der 


Scheibe ist in Abb. 2 dargestellt. 
Zu der Spannungsfunktion 
Forxy (3.17) 


gehért ein besonders einfacher Spannungszustand. Mit Hilfe von (3.8) berechnet man die 
Spannungskomponenten 


Or 
Oy, = 0). Oy — 0, Txy = ~ Ox dy =— (3.18) 
und ; 
2 
So = 0 2 ie = 3G 2 (3.19) 


es entspricht ihr der durch Abb. 3 wiedergegebene Zustand der reinen Schubbeanspruchung 
einer Rechteckscheibe. 

Auch die Uberlagerung einer Schubbelastung mit einer gleichmaBigen Zug- oder Druck- 
beanspruchung 1aBt sich einfach behandeln. 
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by Die Biegung des Rechteckstreifens durch eine Einzellast. 
Als weitere Anwendung der Methode der Airyschen Spannungsfunktion behandeln wir die 
Biegung des Balkens als ein ebenes Problem. Wir beschranken uns hier auf den Sonderfall, daB 
die Kompressions- und Schubfunktion in der Form (2.10) darstellbar sind. Mit dieser An- 
nahme wird aus (3.9) 


il 1 I OF O71 OF 071? OF 4? 
A 2 82 0 0 t 2: 
lieve eee +36 gt) AF] a6 (39 Ox? “i Ox Oy? 5p Oy Ox 35) 


Fae . iL 1 
Nach Division mit (ox | 36 und nach Ausfiithrung der Differentiationen entsteht 


Fense t 2 Fesyy + Fyyyy 
eta Byg a vey = 2 Fy By) Bea c 
air (Fo. + 2 Fy, + 2 Fly — 2 F 2 Fyy) Fyyyy 
+- BaP © Ley) Lx nix 
+ 2 (2 Fyy — Faz) (2 Fay Fayyy + Feyy) | Gy) 
FAO ey) Foye (8 Hy — Fy x) Ee 
= ae ae ly SRS eo 


SY lee ORY 
eae Pee sy 5 — Peay lays) 
Pe Fy (Pave Pary 1 2 Pony Peyy + Fayy Fey y)} =0 


mit der Stoffzahl 


SY 


ee Seek Lee & K 
A pai (aed 30) (3K + G6) G° 2) 


Wir entwickeln zuerst eine Naherungsmethode zur Behandlung der Differentialgleichung 
(3.20). Der Ansatz der Stérungsrechnung besteht darin, daB man die Lésungen von (3.20) als 
Potenzreihe eines kleinen Parameters zu bestimmen sucht. Die Spannungsfunktion F(x, y) 
wird hier als Potenzreihe in A angesetzt: 


B(x, y) = FO) (x, y) + A FO (x, y) + 2 FO) (a, x) foe. (3.22) 
Durch Einsetzen der Reihe (3.22) in (3.20), Entwickeln nach Potenzen von A und Nullsetzen 
der Koeffizienten der aufeinanderfolgenden Potenzen von / entsteht eine Folge von linearen, 
inhomogenen Differentialgleichungen, die man nacheinander integrieren kann. 
Die von / freien Glieder ergeben die Differentialgleichung der ,,ungestérten“* Aufgabe, die 
Bipotentialgleichung 


AAF = 0. (3.23) 

Der Koeffizient von /! liefert als Differentialgleichung fiir F(): 
( ( (0) 
AAF® + FO! FU),, + FO) Fe 


Speke fel 
— 2(Fie + Fyy'— 3 FDL Fyy — 6 Fy) Fotyy 
+ 2(2 FLL — Fy) (2 FS, Foley + Fe) 
SAB ED ICA DT ara) | (3.24) 


+ 2(3 FY — FS) FOS + 28 FY) — Fy FY: 
$2RY PY FO, +2 PO, FO, 
oer sae Baye i Peng Vaya) a () 
0 
a= 12 PE ey Lae 4 ple Ee aie Beyy Py yy) =0. 
Durch Nullsetzen des Faktors von A? folgt 
AAF® + (2 FRY + Fy + 2 Fay — 2 FY, Fyy) Fates | 
op (FEY p 2 Fyy 2 FOP — 2 BS) BP) Evy yy + 2 (Fee Fey + 8 Fey) Petey 
+ 4(2 FS, — Fy) Fey Foley + 4(2 yy — Fes) FE, Feyyy (3.25) 
42 BO FO) ro FO), 1 4(2 Fo) — FO) Fe Fo) 
si 4(3 fee a Ee) Be, Ee ie 4(3 eh re. Ey) 1 ie joes 
+ 2(FY PY, — POP, + 6S Pl, Fe, 
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+ 2(FEL IY) yt 6 Fay Fe be + 13 Fey Feyy) Petey 
Lo( HO) BL) Obey ae 1S FOE ae 
ELE 2( Fe) Fe — Fe FOO 3 6 1) F,,) Ele 
see y(2 Fo, =e Fe) ,.) FO) + 3 FY, Fi Fe) 22 Fo, ey he, 
sls 9 pos ts 3 Lah a ie, ae Fiery "FY, oy Eo Fo) 3 95) 
+ 2[(2 HANS el oe a) Se) a) oe 3 FY ley OL 2 (2 FO yy — sew ge , 
— FO? — FOR, +3 FO, A FOL RO) ee Ry 


i 


Hao Ie ae Fo y+ 2(2 ie iss Fp) Fo, S23 Fo, cae Fey) oe 
Se hen DB leash 7 , 
0 0 
+ 6(F§ sire Eo aE eke a ayy Eye) ee J 


Die Spannungskomponenten stellen sich dann als oe in A dar; denn durch Differen- 
tiation von (3.22) erhalt man die ato 


ea eS o!®) gone) Regs 
Cas 
Oy = ae + A oy L 2 ye a $ (3.26) 
_ OPF 7) (2) | 
tae yA eye tae oe 


3 ris a2 FC) Aye o2FQ) ome o2F (2) 
dy? 2 x% ay? 9 x dy? ) 


die zweiten partiellen Ableitungen der Naherungsfunktionen F, F“), F® usw. bedeuten. 
In gleicher Form lassen sich auch die Invarianten (3.5) darstellen; wir schreiben 


55 = xz (of + dof + A? of? +:+°) 


mit 

oy Sag See oy) 

of) == (ol) + of), 

si) = (of of. (3.27) 
und 

D5 = (OD +4 GY 4 22 G1. ++) 
durch die Festsetzung 
DB") = oP* + of)? — 6) of) + 3 7)2, | (3.28) 


@ —2 o) o) =. oy) gy ol (0) of) — o) Ube 126 oe To) ’ | 


Sobald die Spannungskomponenten (3.26) und die Invarianten (3.27) und (3.28) bekannt sind, 
kann man die Verschiebungen ermitteln. Mit der Annahme der Funktionen k(sy) und g(t) in 


oa te (2.10) lauten die Spannungs-Dehnungsgleichungen (3.4) des ebenen Spannungszu- 
standes 


ee = sgt yg + go) 0%, | 
oy = 3p Ot oq (lt hei) 0%, (3.29) 
Yuy = (TisiBe to heey 

mit Ox = (G20), Oy = (Gy —%). 


Zur Auffindung der Verschiebungen werden jetzt sowohl fiir die Verschiebungen u, v als auch 
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fiir die Verzerrungen e,, éy und Wx y Ansatze von der Form (3.22) eingefithrt, also 


ee ee eee) 


vy = yp) = Avy aL? p() celta f (3.30) 
und 
ee =H M4 AME BM Ho, ) 
Ov rer 
by =a = by) + A ey + Mey) poe, (3.31) 
Ou dv 
Yay =F + og = Pex t Aye + Pye to 


gesetzt. In (3.29) driicken wir nun alle Spannungs- und Verzerrungsgré®en durch die ent- 
sprechenden Potenzreihen in J aus. Beachtet man die Abkiirzung (3.21), so ]4Bt sich die rechte 
Seite nach Potenzen von 4 ordnen, und durch Vergleich der Koeffizienten gleicher Potenzen 
auf beiden Seiten erhalten wir die Koeffizienten der Entwicklungen in (3.31). 

Fir die nullte Naherung erhalten wir 


= 1 0) 

ef = 3K 00 tag 0”) | 
1 1 | 

(0) (0) “( 

ee 5g oy (3.32) 
1 

( ( 

Pry = EG tey | 


Sie stellt die .,ungestérten‘* Verzerrungskomponenten beim Hookeschen Gesetz (A = 0) dar. 
Durch Vergleich der Koeffizienten von / folgt 


1 il 1 1 
Coe (1) (1) > (0) G0) _) 
bx ure 26 °* (ox 36) 
1 1 1 1 
Geass (1 4( fi @ 
23 3K" 26%" (ox 36)" o> (3.33) 


3K 0K 3G 
De 1 1 ie: } 
ep = 5K 13g (ox zc) (Cy ea 0, OO), (3.34) 


(2) _ 1 _@) 1 1 (1) ( 


potes te G0) oe o/ ( 1 I iow DO + of) PW), | 
| 
) 


Aus den Verzerrungskomponenten werden nun die 
Verschiebungen u,v selbst ermittelt. Die Integration 
von (3.32) bis (3.34) bereitet fiir einfach zusammen- 
hangende Bereiche im allgemeinen keine Schwierigkeiten. 

Als Beispiel wollen wir die Spannungen in einem 
Balken von der Héhe 26 und der Lange / berechnen, 
der an einem Ende eingespannt und.am anderen freien 
Ende durch eine Last Q nach Abb. 4 belastet wird. Die 
als klein anzusehende Dicke h des Balkens senkrecht 
zur Zeichenebene ist gleich Eins gesetzt. 

In der linearen Elastizitatstheorie wird dieser Spannungszustand durch die Spannungs- 
funktion null/ter Ordnung + 


Abb. 4. Hinseitig eingespannter Balken. 


FO Le x(ys—3bty) mit ¢=—+ (3.35) 


1 Vgl. S. Timoshenko, Theory of Elasticity, 1. Aufl., S. 33, New York 1934. 
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beschrieben. Nach (3.8) erhalten wir hieraus die Spannungskomponenten 


oo) = oy a0 mre, = ey te (3.36) 


die den Randbedingungen der Aufgabe geniigen. Denn 


fir_ y = -—b ist Gay a0 


fir « =0 ist OL Or; 
4b 
auRerdem ist dort die Resultierende der Schubspannungen im Querschnitt J 7) dy =— Q. 
Fiir die erste Naherungsfunktion erhalt man nach (3.24) die Gleichung 
AAF) = — 2 c3(2 nw y + 8 x y? — 5 Bx y) (3.37) 
mit den sog. homogenen Randbedingungen: 
in ‘y= tb muB “oS = 7 = 0. 
4b 
in a —=0 mu go.) 05 J ty dy = 0 | (3.38) 
sein, Durch Einsetzen bestatigt man, daB sich ein partikulares Integral von (3.37) in der Form 
Fo Smt no 2 xy? — 25 b2 x y8) 


darstellen 14Bt. Zu dieser Partikularlésung nimmt man eine Lésung der homogenen Gleichung 
AAF = 0 hinzu, welche aus bekannten Pol ynomlésungen der Bipotentialgleichung * zusam- 
mengesetzt ist: 


1 1 1 1 
HO ie(z it ero: “ 9) + bisa x8 y+ dys xy? +b xy. 
Die hierin auftretenden Festwerte d,, b,, d, und b, sind derart zu bestimmen, da die Funktion 
FQ) = FO + FP) simtliche Randbedingungen (3.38) befriedigt. Nach kurzer Rechnung findet 
man die erste Naherungsfunktion 
3 
one Sule xy A xy? — 4b aP y8— bam ys 2d ary oe bt y8— 2 iF x) - (3:39) 


Die Verbesserungsglieder erster Ordnung der Spannungen errechnen sich dann zu 


of = —© gy (10 xy? — 6b? x? + 6 yt — 13 Bey? + OF 
sees a) Ms Bb Se ege a) | 
J 


3 
1 = 5 (30 2 x4 — 36 b? a? y2 + 6 BA a? $d 96 — 13 Be yt 4 200 pt ye vs). 


Nunmehr kann die Differentialgleichung fiir die zweite Naherungsfunktion F@) nach (3.25) 
angeschrieben werden. Sie lautet 


5 
AAF®) =] (400 x y3 +. 1858 x3 y5 + 944 x y? — 72 b? x* y — 2000 Bb? x3 y3 — 2298 b? = 
(3.41) 
At 12486 eae 57434 Paleo) | 


35 35 


mit den mit dem Zeiger (2) geschriebenen homogenen Randbedingungen (3.38). Man findet zu- 
nachst leicht ein partikulares Integral dieser Gleichung. Sodann fiigt man hierzu eine Lésung 
der homogenen Gleichung, welche eine geniigende Zahl von verfiigbaren Konstanten enthalt, 
um die Randbedingungen des Problems befriedigen zu kénnen. So erhalt man folgende 


1Vegl. S. Timoshenko, a. a. O. S. 27. 
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Funktion: 
BLS vig ( 10a 9b #98 5, a yt — 98 52 a8 ys — 38 be a8 y? — 3 pea ys 
— bE ad 9 $08 bt go ys 4 FNP pty? 
PE TBE tt I = 
tM ny pa). 


Hieraus leiten sich nach (3.8) die Korrektionen zweiter Ordnung of), of) und t£) fiir die Span- 
nungen ab, die wir aber hier nicht anschreiben wollen. 

Die errechnete Spannungsverteilung ist nur dann streng richtig, wenn die Belastung am 
freien Ende durch stetig verteilte Krafte, entsprechend dem ermittelten Verlauf der Schub- 
spannungen 7t,y = 7%) + Aro) + 22 72, +-++, auf den Endquerschnitt iibertragen wird. 
Greift die 4uBere Kraft Q nicht in der vorgeschriebenen Weise an, so wird der Spannungszu- 
stand nach dem St. Venantschen Prinzip nur in der unmittelbaren Umgebung des freien Endes 
beeinfluBbt. 

Wir behandeln noch die Forminderung des Balkens. Dabei nehmen wir an, der Mittel punkt 
des Endquerschnitts sei so festgehalten, daB folgende geometrische Randbedingungen gelten: 


Ov 


0 
Ox 


10 PUL eh ey (3.43) 


Mit den in (3.36) gerechneten Spannungen erhalten wir nach (3.32) die Verzerrungskompo- 
nenten fiir die Naherungslésung nullter Ordnung 


dul) 1 1 dv) 1 1 

(0) _ (0) j= ae 

ex Ox =(sx+a¢)¢*7: 2) dy (ax sg) ony 
du) = av) c 

CO) mea 2 2 


woraus die Verschiebungen nullter Ordnung! folgen: 


Cy a. 3 Ise 5 ; 
(esp ane ack BP TP sages leaf SEE CX a ete 2 2 
re 3 Py) sea ("9 ao Py +3bey), | ae 
x 2 ec eae ae BT 2 
= oR (#9 = 1 J? % a )-sea ert geht ah), 
Nach (3.33) findet man leicht die Verbesserungsglieder erster Ordnung der Verzerrungs- 
komponenten in der Form 


e) i aK 60° xy? — 24 b? ny — 16 B xy? + 2 bt x yl 
gee 12 x y> — 2482 wy + 20 xy — FP bt ay), 

ef) =~ = sK50(% x? y® 4 51 xy — 240% x9 y — 106 B® x y? 4 A077 ps 9) 
—30 wl 15 xy> + 12 B? my — btw yt +e bey), 

yo) = a s See a yt 1 45 y8— 60 b2 a2 y2— 135 B® y4-4 135 bt y2— 4558) 


1 3 
Sera 30 x? yt 4 33 y+ 48 B? x? y? 


me 4 42 2 
— 96 b? y 18 b4 x + ge Oty 35 


B62 tee aes PAVAL6 v) 


1 Vel. S. Timoshenko, a. a. O. S. 36. 
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Die Integration dieser Gleichungen liefert unter Beriicksichtigung von (3.43) die Verbesse- 
rungen erster Ordnung fiir die Verschiebungen 


wi) = — LO (3 92 ys — 6 Be xt y — BBY ay? 4 22 BE ay 
9K 60 ee y ve 
73 217 2591 51 
Tayioe TH bey. | A yas oy DoT) 
Mon — 6 Be xt y 4 10 be a? y? — "08 8 a? y 
61 4 1116 366 
oF yt — 2 oe ys + ey? 5 taney 
le Li 53 Loi, aes 
wt ee gots ot ay eee eee ta yi t+ —— bt ay? 
6 17 4 51 4 ]2 2 5 34 4]3 
sags Cem yeh iar Osa as Det ne ag bl 
1 8/5 579 
See ula te Oe ale le tea 


6 120 oe 366 14 Danes EDA are 
fe OP ae DE OD ee Pa = bl aaa - | 


In gleicher Weise berechnet man die Korrektionsglieder zweiter Ordnung, deren Richtigkeit 
leicht nachgepriift werden kann: 


7) aKa (2 xf y? sree ee b? x4 y® ke x2 v7 — bt yt 129b* x4? 
a es ee Peecore ti aD see yey 40127 oe, ba? y 
A ae Ye 955 by? bat bt y? ats be 5 sn by? 
egos) Oe, OD segs Ea 
assis AS x4 ype x y9— 21 Beat ys 1 Baty? 72 bt xt y 4 165 bt ty? 
= Ti cae aay es oe a sane re 
all a yi ot yy? + oe ty? — ae ecoe Be by? 
a, Foye ny a 6 hy — ss ie 
y(2) — xx as(14 x 6 ae x3 8 ma ae b2 x9 y4 413 b? x3 ye — 2059 pe, ye ee“ 
8 
a = ps x yy? “ye, yt ee ee a be x3 y2 
5051, yt oe yrs eee ee 
— SU +, BH pe EEE ae ory BE + aay | 
at sean ie x9? — 2 oe yO 1488 x ¥8 += Bee y® + 2 bt ad ye 907 4 xs ya 
+ bt ais re B® x3 y? — a8 bo x y4 + ee be x y? 


24 2612 18551 
oo a ee 6 45 
T, 33 375 & * 2000 


55653 144 Ss 18551 
Ze BSP x + = ba — ps — 8 13 
2000 a b 1000 st 


O8 a? — 2 DAD a + OE" Bt 


— a 
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Zahlenbeis piel: EsseiQ = 33,3 kg, b =5cm,! = 50cm; Werkstoff: Kupfer + mit 
K = 1,33- 106 kg/em’, G = 0,47- 108 kg/cm’, g. = 7,26+ 108, also nach (3.21) A = 0,98- 10 
em4/kg*, Es soll die Normalspannung o, an der Stelle x = 1, y = — b ermittelt werden. Wir 
erhalten 


oz) = 100,0 kg/em®, Aol) = — 2,60 kg/em®, 220%) = + 0,42 kg/em?. 


Mit diesen Werten ergibt die zweite Naherungslésung der Spannungen gemab (3.26) den Rand- 
wert (fir x =—1, y = — BD) of) = 97,82 kg/em?. Die maximale Normalspannung verringert 
sich demnach gegeniiber dem Wert der linearen Theorie um 2,18 %. 


Nun ermitteln wir noch die Durchbiegung fiir das belastete Ende (x= 0; y = 0). Wir 
erhalten 


v = 13,22-10-3em, Avi!) = + 0,31-10-3cm, 22 v@) = — 0,47- 10-5 cm. 


Mit diesen Ergebnissen liefert die zweite Naherungslésung fiir die Durchbiegung an der Stelle 
x = 0, y =0 den Wert vo!) = 13,52- 10-3 em; also betrigt der Unterschied gegeniiber dem 
Wert der linearen Theorie bei gleicher Belastung + 2,28 %. Wir werden aber in spiteren Bei- 
spielen einen zahlenmafig viel gréBeren Einflu8 der Nichtlinearitat des Elastizitatsgesetzes 
finden. 

In gleicher Weise mit Hilfe der Methode der Stérungsrechnung 1] aft sich der beiderseits auf- 
liegende, gleichmaBig belastete Balken behandeln. (In der Dissertation sind die Verbesse- 
rungen erster Ordnung angegeben.) 


4, Der ebene Spannungszustand in Polarkoordinaten, Fiir viele Aufgaben ist es bekanntlich 
vorteilhaft, ebene Polarkoordinaten r, m an Stelle der kartesischen Koordinaten x, y zu ver- 
wenden. Der ebene Spannungszustand einer Scheibe im Punkt (r, y) ist dann durch die Nor- 
malspannung g, in Richtung des Radius r, die senkrecht dazu stehende Normalspannung o, 
und die Schubspannung f,,, bestimmt. 

Da die Darstellung der Spannungen durch eine Spannungsfunktion unabhangig vom Koor- 
dinatensystem ist, kénnen wir diese Spannungen aus einer Spannungsfunktion F(r, y) folgen- 
dermafen ableiten: 


_1@F , 10F 

eae ap? Pere | 

OF 

a (4.1) 
Ea | 

= lee : 


In (3.9) haben wir die Differential gleichung der Airyschen Spannungsfunktion in rechtwink- 
ligen Koordinaten hergeleitet. Mit Beniitzung der bekannten Transformationsformeln? auf 
das Polarkoordinatensystem umgeschrieben lautet sie 

1 [/1@F , 1 AF\@g(2) , @F/1 Og(12) , 1 agit) 
Ge Og? | r Or} Ore? | Or®\r2 ag? | r ar | 42 
1 @F 1 4F\/1 og(3) 1 agi8)\] _ 9 
( Ordp 2 Op)\r ar dy r? 0p 


mit 
Il eee! ee C2 a So ee oF 
80 = 9K (Or tO) =9n4F =5K( Ge bat age Tar) 
2 ‘ 
6 = 9a (Or + 9 — 9, Oy + 3 tre) (4.3) 


il elm IRORN? Or fA oe a! oF 
= 1G (44 a(7 Pp ay D ar) Or? (3 OG? T x ar)| 0 


1 Siehe H. Kauderer, a. a. O. S. 468. 
2 Vel. etwa K. Knopp, Einfiihrung in die héhere Mathematik, Bd. II, $.360, 8. Aufl., Stuttgart 1947. 
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Beschrinkt man sich auf den schon im vorigen Abschnitt behandelten Sonderfall, daB die 
Kompressions- und Schubfunktion in der Form (2.10) darstellbar sind, so vereinfacht sich (4.3) zu 


AAF+A [Free (Fe eh: 7 ee & cs ae 2(7 i F,) | 
eid), le ig fat te) pee 8(- ee F,) 
eT ee (F = =F, — Fe) jae 2(- abe F,,) 
anda 
+5 Fere(> Fea Fe) 


4 1 1 2 2 
Belew G os F,) (Fe —in—-+ es) 
1 


1 1 
Ue UG eer tamer F 9] 


i 2 Fre (2 Fr ~Lh—| Fe) + Feo(3 EF. —iR=G He 

3 Fyo(Fr > F, — 5 Feee(Fe Fhe i) 

+7 Foor Fee a Fr tae Fog) + zt Fore Pope Fo 

pA ci id Me (ape F,) — Free Foray i F,) 

$2 Byg (Ly DE) + Bey Real Dy oF ; 


2 2 2 ‘ 1 3 2 
FFP ih a Pe Pg tah tsk Fp + Fy 


3 


fl 1 6 
Hs eat eae 


2 / 14 15 9 9 4 4 
— at Pele Ne eed Ge eV cee oe Oe ea oli 


2 : 2 6 6 34 74 40 

+ 33 Fegg| Fhe > Fre Fr— SEF, —% Fo — = Flot a Fig Fy 7 By 
2 4 6 6 10 12 16 

+E peep Berkey Fer By pak Feg t oly —ix ig log tz Fo Pye 
1 3 2 4 oe 6 

+ a ie +a Fy, Fp —3 i, Foot Farle ta Fr Ue 


60 18 2 26 2 
—-F,, F,,F,—=F, Fy +4 Met zh Rta lily, 


32 52 10 0 Whee 
+55 Foti hehe ta hie — a EM — a Fre My Poe 
es 8 
a ietee a Eee 0 ! 
mit der Stoffzahl 
ms ER 1 Te fa dS . 
4 ge ‘(ox +3¢) BK+6e° >) 


a) Einfache Lésungen in Polarkoordinaten. Wir wollen uns zuerst 
den einfachsten exakten Lésungen der Differentialgleichung (4.2) der Airyschen Spannungs- 
funktion zuwenden. 


Der Spannungszustand, der zu der Spannungsfunktion 


He G yA (4.6) 
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gehért, ist von y unabhangig. Gemaf® (4.1) erhalt man die Spannungskomponenten 


wees a AO, ONeee t—  Ce Tre = 9 (4.7) 
und nach (4.3) die Invarianten : 
1 2 
%=oKAC: G7? O- (4.8) 


Die Funktion (4.6) mit (4.8) geniigt somit der Differentialgleichung (4.2). Sie stellt den Span- 
nungszustand einer in ihrer Ebene auf allseitig gleichen Zug von der Grife 2 C beanspruchten 


Scheibe dar. 


Wir untersuchen den Spannungszustand, der durch die Spannungsfunktion 


F=Co (4.9) 
beschrieben wird. Die zugehérigen Spannungskomponenten sind 
Cc 
= Oy al) c Tro ey (4.10) 


und damit kommt nach (4.3) 


%=0, #242. (4.11) 


2 
Wegen a = 4 a =0 ist die Differentialgleichung (4.2) erfiillt. Der durch (4.9) 
gegebene Spannungszustand ist von der Richtung y unabhangig und stellt die Bean- 
spruchung einer Vollebene durch ein Moment M nach Abb, 5 dar. Die Konstante C berechnet 
sich aus dem Gleichgewicht zwischen M und den Schub- 
spannungen T,,, an einem herausgeschnittenen konzentrischen 
Kreis von beliebigem Radius r, namlich M = 2 qr? Tre > 


wora Cla ———_ 

b) Die Kreisscheibe We \\ 
mnter gleichmaBigem 
AuBendruck. Wir behandeln | 


jetzt den Spannungszustand in 
einem geschlossenen, scheiben- \ 

férmigen Kreisring, der nur an SS 

seinem AuBenrand durch eine gleich- ss 
eee ute Normalbelactung tev emar | ABR E. Kreisringscheibe unter AuSendruck, 
nach Abb.6 beansprucht wird. 

Die Halbmesser seien mit r, und r, bezeichnet, die gegen die Abmessungen in der Mittelebene 
kleine Dicke h = 1 gesetzt, und der gleichmaBige Druck am Scheibenrand r = r, besitze die 
GréBe p. Es gelten also die Randbedingungen 


(GC; Jee fy = oF (0,), =I, > Ps (4.12) 


Die Belastung ruft einen drehsymmetrischen Spannungszustand hervor. Da dieser Span- 
nungszustand vom Winkel y unabhangig ist, verschwinden die partiellen Ableitungen der 
Spannungsfunktion nach 9; fiir die Spannungskomponenten erhalten wir die Bezeichnungen 


Eelear aE 


eign Pt art 


SS 


%, =0, (4.13) 


und die Differential gleichung (4.4) nimmt nach kurzer Umformung die Form an 


(Ar? FP’? — 2 Ar FY FA? 4 7) 72 FV | 
4 24(2r FY — F) 3 FV 4 202 BY — BY) BF 2 (4.14) 
4+ (AF? —34r FoF 4 2A0F!— 2?) F’ 47 F =0, 


wobei Striche Ableitungen nach r bedeuten. 
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Man kann nun ein erstes Integral dieser Gleichung ermitteln. Es lautet 
(2Ar2 Fe? —2Ar FF 4) FP? 1?) FY | 
(Ar FR fh Fee) a 

mit einer willkiirlichen Konstanten C. Man bestatigt es leicht, indem man das Integral (4.15) 
nach r ableitet und die Konstante C zwischen den beiden Gleichungen eliminiert. Analog zu 
der linearen Theorie stellt (4.15) fiir den Fall C+: 0 die Differentialgleichung der reinen Biegung 
eines Kreisringsektors und fiir den Sonderfall C = 0 diejenige einer durch gleichmaBigen 
Au®en- oder Innendruck beanspruchten Kreisringscheibe dar. Um dies fir den Sonderfall 
C = 0 zu beweisen, ist es nétig, auf die Formanderungen einzugehen. 

Wir haben also fiir C = 0 die nichtlineare Differentialgleichung zweiter Ordnung 

(2QArZ2—2ArAZ+ALZ2+9r)yrZ” \ 
DOTA Amy Ae WA V5 ==) |) 
zu lésen, wobei zur Abkiirzung Z(r) = F’(r) gesetzt ist. Wir wollen die Integralkurve von 
(4.16) in der Umgebung der Stelle r = rp durch eine Potenzreihe darstellen. Zu diesem Zweck 
setzen wir in (4.16) r = @ + ry und machen fiir Z(r) = Z(r — r) = Z(@) den Ansatz 
Z(e) = 40+ 0? + a0? Tes: (4.17) 

wobei wir das erste Glied a, = 0 setzen und somit von vornherein die erste Randbedingung 
(4.12) erfiillen. 

Durch Einsetzen in (4.16), Ordnen nach Potenzen von g und Nullsetzen der Koeffizienten 
gleichhoher Potenzen von g bestimmen sich die Koeffizienten a, dieser Potenzreihe rekursiv 
aus folgendem Gleichungssystem: 


27(2A0a?+ lha,ta=0, 

6 r2(2A a? + 1) a, + l6ArzZ a, a3 + 8r (Aaj + 1l)a,+Aai+a,=0, 
12 r3(2 A a? + 1) a, + 3 rZ(24A 1 a, a, + 8A at + 7) a3 

+ 16A73 a8 + 386A r2 a, af + (11 Aaj + 9)a,=0, 


20 r3(2 Aa? + 1) a, + 8 r2(16A ry a, ag + OAat + 5) a, 
+ 72 Ar8 a, a2 + 1)(96 A 12 af + 164A 1 a, a, + 25 A af + 23) a, 
+ 40 Ar? a3 + 36A 1 a, af + 3A aj+ 1lha,=0, 


30 73(2 Aa? + 1) a, + 5 re (40A ry a, ag + 16/ af + 13) a, 

+ 15 (240 A 12 a, a; + 160A r§ a3 + 292 A ry a, ap + 45 A af + 43) ay 

+ 12472 (15 ry a, + 14 4,) a2 + 2(122 A r§ a3 + 715A 1 a, a, + 4/ af + 4) a, 
+ 38A7r, 08+ llAa, af =0, 


(4.15) 


(4.16) 


usw. 


Die Auflésung dieser Gleichungen ergibt 


i I a, 
¢ 2r,24a? +)’ 
1 a, 
ay Grae aie Cie, ome 
1 a 
4 =F at pays (16 ¥ al? +. 32 A4 af + 28 A3 af + 8 A? af — 3), 
1 a 
% = — Toys Hae Tap (928M? alt + 3184 18 aj?-+ 4960 7° aj? + 4248 4 a’ (4.18) 
+ 2410 A3 af + 803 A? at + 225 A a® — 60) , 
1 a 
As = a3 7 WTat ip (13824 2° al® + 67392 48 a}? + 146816 17 aft + 183152 28 a? 


+ 149360 /5 al® + 81324 A a8 + 34464 23 af + 10497 22 af 
+ 4005 A a? — 360) , 


usw. 


XXII. Band 1954, Jindra: Einige Anwendungen eines nichtlinearen Elastizitatsgesetzes. 135 


Die Lésung ist nun der zweiten Randbedingung (4.12) anzupassen. Man erhiilt 
1 dF 1 


a Ta ar F 


(0,)rar = “(= %) =— p 


a= Ta 
= Tq 


oder gemaB (4.17) 


Fe Nae oh — PP; (4.19) 
“= 
womit eine Bestimmungsgleichung fiir die Konstante a, gegeben ist. 
Die Spannungskomponenten go, und o, findet man endlich nach (4.13) und (4.17) zu 


é ey (4.20) 


Zahlenbeispiel: Es sei p = 152,78 kg/cm? ; r,/ry = 1,2 ; Werkstoff: Kupfer } 
mit K = 1,33- 10¢ kg/em?, G = 0,47- 106 kg/cm? , go = 7,26: 10°; also A = 0,98- 10-5 
em4/kg?. Nach dem Newtonschen Verfahren bestimmt man die Wurzel von (4.19) zu 


4, = — 923,9 kg/em?. Aus (4.20) erhdlt man die Ringspannung o, am Innenrand (6,), — ;, 
= — 923,9kg/cm®, und am AufSenrand (0,,)r=r, = — 907,8 kg/em?. Die entsprechenden 
Werte der linearen Theorie sind 

(03), =7, = — 1000,0 kg/cm? , (O>)r =r, = — 847,2 kg/cm? . 


Die Spannungsverteilung der nichtlinearen Theorie weicht also am Innenrand um — HOLE 
und am Auffenrand um + 7,15 % ab. 

c) Der ebene Keil mit einem Moment 

-an der Spitze. Wir beschaftigen uns hier mit der 
Spannungsuntersuchung in einer von zwei sich schnei- 

‘denden Ebenen begrenzten Scheibe, welche sich ins 
Unendliche erstreckt, einen halben Offnungswinkel 
B <a/2 besitzt und an der Spitze durch ein Moment 
M nach Abb. 7 belastet ist. Die Scheibe ist in r = 00 
gestiitzt zu denken. 

Als Naherungsverfahren zur Lésung von (4.4) wenden 
wir die in Ziff. 3 bewahrte Methode der Stérungsrechnung 
an. Wie beim eingespannten Balken gehen wir auch ALE Gu Bhoner Keiltunif ciucad Moments 
bier von der bekannten strengen Lésung nullter Ordnung 
aus. Die Spannungsfunktion F, die der Bipotentialgleichung geniigt, lautet? 


(0) — sin2Q _ 9 ; < M 3 
F (cep YP mit c (2h —te26)° (4.21) 
Die Spannungen haben folgende Werte: 
(0) ss A csin 2p (0) __ ¢ (0) __2¢/cos 2 
O; r2 cos 2p 9 Og 0 o) Trop 72 le 2B ‘SE (4.22) 


Die Randbedingungen sind streng e-fiillt; denn aus (4.22) folgt of = m= 0 furg = +B. 


Mit dem Ansatz der Stérungsrechnung F = F(® + . FO +..-erhalten wir aus (4.4) 
fiir die erste Naherungsfunktion F 


2c)? sin 2 
AAFW) = 128 a (4.23) 


und die Randbedingungen 
gia 0 fir p= 4B. (4.24 
Eine partikulare Lisung ]4Bt sich sofort angeben. Sie lautet 


ubioee (2c)? sin2 
P53 cos 26 _-r4 


1 Siehe H. Kauderer, a. a. O. S. 468. 
2 Vel. L. Féppl, Drang und Zwang, Bd. III, S. 34, Miinchen 1947. 
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Durch Uberlagerung einer Lisung der homogenen Gleichung kénnen die Randbedingungen 
(4.24) befriedigt werden. So finden wir die erste Naherungslésung in der Form 

Fo = {C0 (sin 2 9 — Bsin 4g — Dsin 69) (4.25) 
mit 

2sin 8B—4sin 48 
~ gin 108 —5sin 2p ’ 

sin 6 B —3sin 28 
sin 10 fp —5sin2f - 


Die Verbesserungen erster Ordnung der Spannungen nehmen hiermit die Werte an 


DS 


4 (2c)? 1 : ; ) 
(1) Ss - 
Oo; = 3 cos2fhre (2 sin 2p 5 Bsin 4 10 Dsin 69) , 
20 (2c)? 1. . 
Ofte ore 2 a ema (4.26) 


ORO aa = 
\re= 3 cos 2fr8 (cos 2p — 2 Beos 4 3 Dceos 6¢). 
Somit haben wir die Spannungsverteilung in erster Naherung fiir den ebenen Keil bestimmt, 
der durch ein Moment an der Spitze belastet ist. Als Anwendung behandeln wir die Halbebene 


mit einem Randmoment nach Abb. 8. Mit B = 2/2 erhalten wir die 


ie Naherungsfunktionen 
j LY, FO =—™ (sin2p +29); | 
; 2 I 
I VING Te ye leas (4.27) 
rf G FY =—3 (4) (sin 29+ sindg +z sin69), 


| woraus die Spannungen nach (4.1) leicht abgeleitet werden kénnen. 
Mit Hilfe der Methode der Stérungsrechnung kénnen zahlreiche | 

airs Se ee Beispiele behandelt werden, fiir welche geschlossene Lisungen in- 
der linearen Theorie bekannt sind. (In der Dissertation wurden 

die Biegung des Kreisringsektors durch eine Einzelkraft, der ebene Keil unter stetiger Riicken- 
last und die kreisgelochte Scheibe unter Zug untersucht.) 


5. Die rotierende Scheibe. Als weitere Anwendung der nichtlinearen Theorie zeigen wir die 
Berechnung rotierender Scheiben. Sehen wir den Kérper in dem frither erlauterten Sinn als” 
eine ,,Scheibe“ an, so geniigt es fiir geringe Scheibendicken, das ebene Problem zu lésen. 

Fiir eine mit der gleichbleibenden Winkel geschwindigkeit w umlaufende, achsensymme- 
trische Scheibe von der Massendichte @ gilt, wie das quasistatische Kraftegleichgewicht an 
irgendeinem Scheibenelement im Abstand r von der Drehachse zeigt 1, die folgende Beziehung 
zwischen Radialspannung o,, Tangentialspannung o, und dortiger Scheibendicke y: 

d 


FRAUC Re) ae pea Se eM: (5.1) 


Fir die radiale Verschiebung u hat man die Spannungs-Dehnungsgleichungen des nicht- 
linearen Elastizitatsgesetzes (2.7) 


du 
e = at = 5a Fle) (a, + 0,) +a al) (2 0, — 9) » | 


gp = tt = gg B50) (61 + 09) + gaz BUR) (2% — 97) | 
mit 
= 1 2 2 2 2 
9 OK (Gai. yes is = 1e (Oni Oot 0; Oy) : (5.3) 
Durch Entfernen von u aus den beiden Gleichungen (5.2) ergibt sich die Vertraglichkeits- 
bedingung 


1 d 
wt [mls (0, + %9)] — ge Te [Bl8) (or — 2 6) — gq Bl) (6 —9,) =0. (5-4) 


1 Vel. C. B. Biezeno und R. Grammel, Technische Dynamik, 2. Aufl., Bd. IT, S. 5. 
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Die simultanen Gleichungen (5.1) und (5.4) lésen das ebene Spannungsproblem der 
rotierenden Scheibe bei gegebener Profilform y (r) vollstandig. 

a) Die Scheibe gleicher Festigkeit. Wir suchen hier diejenige Profil kurve 
y(r), bei der das Material an jeder Stelle die gleichen Festwerte fiir die Hauptspannungen 
besitzt. Mit dem gemeinsamen Festwert 
0, =O, =0 5.9 
wird aus (5.3) : ae 
ee 2 0 


== Konst.., 0-9@ 


= Kkonet., 


und durch diese Werte wird (5.4) befriedigt. Die Gleichung (5.1) gibt dann die Profilkurve 
wieder 

ay oa 

aes 
oder integriert 

ow? 

bast aad 
Die Integrationskonstante yg stellt die Scheibendicke fiir r = 0 dar. Die Lésung ist also von 
-der Form des Elastizitatsgesetzes unabhangig. 


= (5.6) 


Nach (5.2) sind die radiale und tangentiale Dehnung konstant und einander gleich: 


at, Ee Lal 2ZON eo > 207\ o 
c= 8 = % =i aRl oR e(se)eo- (6.7) 
Die Scheibe gleicher Spannung stellt daher auch eine Scheibe gleicher Dehnung dar. Die 
radiale Verschiebung 
20\20 207\ o 
gine Hox) 9K 8(sa) ec ; 
nimint gleichmaBig mit dem Radius zu und erreicht am AuBenrand den gréBten Wert. 


b) Die Scheibe gleicher Dicke. Da die Scheibendicke y unveranderlich ist, 
‘fallt sie aus (5.1) heraus, und man erhalt das simultane System 


d Sinn yee 
peer 6, ee OT 0, | ss 
1 d 1 d il 5 
KT gp Lh (0) (9; + 9,)1— ge ¥ z- [a (86) (0 — 2 9,)] — xq 8 (4) (0, — 9,) = 0. | 
Der ersten dieser Gleichungen geniigt man, wenn man mit einer Funktion Z(r) 
ate _ @ 2 72 
Camis Cee arr S (5.9) 


-ssetzt. Beschranken wir uns wieder nur auf den Sonderfall (2.10), so ergibt die zweite Glei- 
chung (5.8) fiir Z(r) die nichtlineare Differential gleichung zweiter Ordnung 
[2Ar?Z’2@— 2ArZZ’4+AZ?+Ahow* r* Z’— 2) ow? r? Z+ 2A 7? (gw?r?)?+ 17] 7 Z” l 

+fArZZ/’—12?2+ 5iow'rZ’ —4how' re Z+ 10Ar? (9 w? r?)? + 3] 2 (5.10) 

+ [2Agu%r® Z—5Ar(qw? 1%)? —r]Z +51 12(0.w? 1)? + (3+) ow rt= 0, 
-wenn Striche Ableitungen nach r bedeuten und folgende Abkiirzungen fiir die Stoffzahlen 
eingefiihrt werden: . 

Dae, 82 ( 1 1 a 8 K 


9@ ‘lok +36] ~@K1O@’ 
1 ei ela ls Ka 
=(ox sc)’ (ox | i, 23 ko 
-wobei y die Poissonsche Zah] der linearen Theorie bezeichnet. 

Da A der einzige Parameter von (5.10) ist, der die Abweichung vom linearen Elastizitats- 
gesetz enthalt, so wird man eine Liésung dieser Gleichung in der Form einer Potenzreihen- 
entwicklung nach positiven Potenzen von A ansetzen. Eine notwendige Voraussetzung fiir die 
_Anwendbarkeit des Verfahrens ist, daB der gewahlte Parameter auf hinreichend kleine Werte 

10 


(5.11) 
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beschrinkt bleibt, um die Konvergenz des Verfahrens zu verbiirgen. Wir machen also den 


Ansatz 


Z(r) = ZO (r) + AZO (r) + A? ZO (r) + --- (5.12) 
Geht man mit (5.12) in (5.10) ein, so erhalt man durch Vergleich gleichhoher Potenzen 
von A die folgenden Differentialgleichungen: 


fiir Z()(r) : 


r2 Z(0)"/ + r Z(0)/ — Z(0) +: (3 + y) 0 wr=0, (5.13) 
far Z") (r) : 
— [2 r? Z0)72 2 7 ZO)" Zo) 4. ZO)2 4. 4 0 w? r4 Z(0)’ 
nae). 0 w? r? ZO) + 2 r2 (o w2 r2)2] Z0)”" (5.14) 


=: |Z zo — + 202 + 59 wr? ZO” — 49a? 1? ZO + 10 r(gu? ry Z00" 


Se 


— [29 wr Z — 5(9 w? r?)?] ZO , 
fiir Z(?)(r) : 
2 ZA Zl — Ze?) ) 
= — [2 r? ZO — 27 ZO" ZO) + ZO? + 4 w? rt Z(0)’ 
— 29m? r? ZO + 2 r2(9 w? r?)?] ZH)” 
= 4 72 Z00" ZO 2 ZO” ZO + 49 wrt ZO” 4 2 ZO’ ZO == 702 


Syl 
+ 10 9 w?r? ZO’ — 49m? r? ZO + 107 (0 w? ry ZO)" ( ) 


uf fee 27 ZO” ZO’ 4.2 20 ZO — Four ZO" 4 Zor? 2+ Zio Zo 
x 
—4ow?r? ZO’ +49 0? r ZO — 5(oa? ry ZY) , 


Diese Gleichungen sind linear und nichthomogen, ihre Integration ist mit elementaren 
Mitteln ausfiihrbar. 

Die Funktion Z geniigt der Differentialgleichung (5.13), die sich fiir die lineare Theorie 
mit A = Oergibt. Die Lésung lautet 


Z) = Ayr+ Bye —agw?r (5.16) 


mit der Abkiirzung 
3+ 
eae 


Fithrt man (5.16) in (5.14) ein, so folgt nach kurzer Zwischenrechnung fiir die Funktion 
Z) die Differentialgleichung 


r2 Ziajy a PLANE —= ZO) 
= — (5 — 37a + 95 a? — 84 a) (9 w?)? r7 — (5 — 30a + 44 a?) A, (0 w?)? PP? 
— (3 — 8a) ABg wr? + (11 — 50 a + 60 a%) By (9 w*)? r? + 2(1— 4a) Ay Bow? 
1 1 1 
+ (1— 4a) By w? — + 4 A, B= 12 By. 


Die Lésung dieser Gleichung nimmt die Form an 


ZO) = Ayr + By— 


Bair 5 
— (5 — 374 + 95 a? — 84 a?) (9 w)® 7 — (5— 30 a + 44.0%) 4, (ow)? 
Sa br) 
— (38a) A2o uw? + (11 — 50a + 60 a2) B, (Qu?) Oe 
2g 1 a_i 


2 1 
+ (1— 4a) A, B,ow? rinr — (1 — 4a) Bog w* = Inr+ A, Boo — Oa ,8° | 
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Nunmehr kann die Differentialgleichung fiir di i 
g fur die Funktion Z(2) nach (5.15 
werden. Als allgemeine Lésung findet man dann S ae eae 


1 
22) = A,r + B,— } 
1 
Se 5760240 4,|(5— 300-4 44 a2) (0 2)? r+. 6(3 — 8 a) A,o wr? +12 Bt 
A(t a) Bo ot rin r| + 720 B,|(11 — 50a +60 a?) (9 2)? 78 


+ 84, By — 12 BRS + 8(1— 4a) A,ow% rin — 8(1— 4a) Bg win r 

+ (745 — 9263 a + 46805 a? — 120026 a® + 156100 a4 — 82320 a) (9 2) r™ 

+ 3(740 — 7971 a + 32401 a? — 58992 a? + 40608 a4) 4,(o w®)! r? 

— 5(1946 — 18885 + 69831 a? — 116480 a® + 73920 a4) B,(o w2)* 1? 

+5 (803 — 6626 a + 18428 a? — 17216 a) A2(9 w?)3 r? 

— 20(1069 — 8464 a + 22378 a? — 19864 a) A, B,(o w2)3 r+ 60(65 — 372 a 

+ 528 a2) A3(0 w)? r5 + 30(2087 — 14884 a + 36122 a? — 29880 a3) B2 (ow)? r? 

— 60(479 — 2532 a + 3440 a?) A? By(g w®)? r® + 720(3 — 8 a) Ado w? r3 (ote) 
+ 720(54 — 160 a) Ay B3 9 w+ — 2880 48 B24. — 120(161 — 404 a) BS g wr 


1 ] I 
+ 6720 A? By — 11280 A, Bos + 8496 Bo 
— 240(5 — 50 a+ 164 a — 176 a’) A, B,(o w?)3 ln r 
— 360(11 — 94a + 260 a? — 240 a’) B2(9 w?)? r3 In r 
— 720(6 — 40 a + 64 a?) A?B, (0 w?)? Fein r 
+ 1440(71 — 384 a + 520 a2) A, B2(9 wv)? rlnr 
— 2880(5 — 16 a) Ag B, gw? rin + 1440(11 — 32 a) 42 Beg w—Inr 
— 1440(27— 148 a + 200 a2) B8 (ow)? In r + 4320(1— 4a) Big wt nr 
Lf 
+ 1440(1 — 4 a)? A, B3(o w?)? r(In r)? + 1440(1 — 4 a)? B3(@ ws (In a : 
J 
Nachdem die Funktionen Z(r), Z)(r) und Z(?)(r) bekannt sind, kénnen die radiale und 
tangentiale Spannung leicht ermittelt werden. Aus der zweiten Naherungslisung 
ZY) = ZO + 4 ZO) 4 7? Z?) 


errechnet man die Spannungen gemaf (5.9) zu 


Ae il zm — + (Z0) AZ) Af? ZO), | 
Lie Tr. 

UI (11) (0) (a) (2) cue 

Boe a 2p dZ dZ dZ 


g ieee 0.0 Lif See ear wae 


Die in Z)(r) auftretenden Integrationskonstanten berechnet man aus den vorgegebenen 
Randbedingungen der Aufgabe. Wenn die Radialspannung am Innenrand r = 1% mit dem 
Wert o, = 0,, und am Aufenrand r =r, mit o, = 0,, vorgeschrieben ist, so unterwerfen wir 


die Funktionen Z(r), Z(r) und Z@)(r) den Randbedingungen 


of) = Or, fic r=T,. | 
aay fir r=1Tr,. | (5.20) 
Cao” =O fir r=7, und or =7,, 


woraus die Integrationskonstanten A,, By; A,, B,; bzw. A,, B, eindeutig bestimmt werden. 
Die Durchfithrung dieser elementaren Rechnung soll aber hier wegen der ziemlich um- 
standlichen Formeln nicht angegeben werden. 


6. Die Torsion des Umdrehungskérpers. Wir fiihren in bezug auf die Achse des Umdrehungs- 
kérpers Zylinderkoordinaten z, r, @ ein, wobei wegen der vorausgesetzten Drehsymmetrie 
samtliche Spannungen und Verzerrungen vom Winkel y unabhangig sind. Werden die Span- 

Og 
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an der Stelle z, r, g in diesen Koordinaten mit 0,, 07, Og» Trg» Tps> Var und 


nungskomponenten 
bezeichnet, so sind 


die entsprechenden Verzerrungskomponenten mit €,, &» Ey» Veg? Wyz> Yer 
fiir das reine Torsionsproblem nur die GréBen T,9, Tys' bzw. Yo> Yos beizubehalten und alle 


anderen gleich Null zu setzen 1: 


O, = 0; =O, = Gr = 0. (6.1) 
Dann bleibt als einzige Gleichgewichtsbedingung iibrig 
Otro OTpz 2Tro ae 
or z Oz r Ue (6.2) 


wobei die Volumkrifte vernachliassigt werden. Wenn t, die Oberflachenspannung in tangen- 
tialer Richtung und n die Richtung der nach auBen gerichteten Oberflachennormalen bezeichnet, 
so reduzieren sich die Randbedingungen auf die Beziehung 


Tre cos(T, N) + Tos cos(z,n) =t,. (6.3) 
Die beiden Winkelanderungen »,, und p,, lassen sich durch die Verschiebungskomponente v 
in tangentialer Richtung ausdriicken: 
__ Ov v __ Ov 
Wr p "A ee Yos =3, * (6.4) 
Mit der HilfsgréBe # = v/r gehen die Gleichungen (6.4) tiber in 
0b oo 


SS i — ai 
Por on id Oz” 


(6.5) 


woraus die Vertraglichkeitsbedingung folgt 


@ fil @ fil 
or ie vos) ~ %e & ver) au (8:9) 
Fir die Aufstellung der Gleichungen braucht man noch den Zusammenhang zwischen den 


Spannungs- und den VerzerrungsgréBen. Das nichtlineare Elastizitatsgesetz liefert nach 
(2.4) bzw. (2.7) die Beziehungen 


. tin = G y(y>) Yor > Tox = GY (yo) Vox (6.7) 
mit 
Dg 
Wo iy (Por + Yor) 
oder 
at 2 1 2 
Yor =E 8 (8) Tor » Vos = EG &(to) Gop (6.8) 
mit 
2 


1) = 3 Ga (Tor + Tos) - 


Bei dem hier definierten Problem ergeben sich nun zwei verschiedene Formulierungen fiir 
die Integration, von denen die eine fiir vorgegebene Randspannungen und die andere fir 
vorgegebene Oberflachenverschiebungen zur Liésung geeignet erscheint. 

a) Vorgeschriebene Oberflachenspan Di i | 
BAO i take P nungen. Die Gleichung (6.2) — 


Gh A 
F (gr) “5 (79 Zp) = 0 (6.9) 


schreiben. Diese Gleichung wird identisch befriedigt, wenn die Schubspannungen T,, und Tt), 
aus einer Spannungsfunktion @(r, z) durch den Ansatz ™ "7 


eee: : ab | 
Tr Tor = Roniar > r? Tox = a an (6.10) | 
: 
hergeleitet werden. Wegen (6.6) und (6.8) geniigt ® der Differential gleichung 
Oil elo Otee 1 a@ | 
5 608) a | +e 8 sa) 0 (6-11) 


1 Vel. C. B. Biezeno und R. Grammel, Technische Dynamik, 2. Aufl. Bd. I, S. 125 
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2 — 2 1{/aw?, /ao\: 
0 3G r4 (Fr) 1 (32) |: 
Fiithrt man die Differentiationen aus, so ergibt si i it 
man , ; gibt sich mit der Abkiirzung g’ (i?) = d g(t?)/d(t? 
folgende nichtlineare Differential gleichung zweiter Ordnung fiir die at omuiieatonldn Bes 


, (A@ 3 AD o2@D 
Yi (peer in i aoe 
8(t) ( Or? r or Oz? 


mit 


2 2 (/dD\2 2D DAD PD OD\? PD 2 2 (6.12) 
ae ag 2 aD /a@\? _ /a® 
+ 8) 5 Gz (ze) at" Ge Oy arae t (F:) dz? me ite) + (Fe) | is | 
Die Randbedingung (6.3) geht mit (6.10) iiber in 
1 d@ 
eae am Rand, (6.13) 


wobei s die Bogenlange der Meridiankurve bedeutet. 
b) Vorgeschriebene Oberf] achenverschiebungen. Setzt man (6.5) 


und (6.7) in (6.9) ein, so folgt fiir 9 die Gleichung 
on Ware a ao 
a (v8) 5, | +e [78 Fe] =o (6.14) 


3 o((2)+ (2) 


/ 


mit 


Nach Ausfiihrung der Differentiationen entsteht die nichtlineare partielle Differentialgleichung 
zweiter Ordnung 


ao 3 ao oo 
2 

(v8) (5p =e =e 7) 

Peay 2\(00 Oe 7 5 08 88 OO | (287099 | 1 od //e8\? « (o8\7)) 

as 3 Y’ (Po) ar (5, or2 ee Oz Or Oz Bice) poroaealeanacell = 


mit der Abkiirzung 


d (9) 

Diese Gleichung ist dann fiir die Integration heranzuziehen, wenn der Oberflachenwert 
der Verschiebung v gegeben ist. Der Wert der Funktion # an der Oberflache des Rotations- 
kérpers ist dann ohne weiteres bekannt, und man kann 
(6.15) bei gegebenen Randwerten von @ integrieren. 

Als Beis piel behandeln wir die kegelférmige Welle 
Abb. 9. Den Koordinatenursprung legen wir in die zu dem 
kegelférmigen Wellenstiick gehérende Kegelspitze und be- 
zeichnen mit « den halben Offnungswinkel des Kegels. Die 
Welle werde an den Endquerschnitten (z = 1, und z = 1.) 
durch Torsionsmomente W verdreht, wahrend am Mantel 
keine duBeren Krafte angreifen. 

Unsere Aufgabe besteht nun darin, eine Lésung der Abb. 9) Kegelisanine Welle: 
quasilinearen partiellen Differential gleichung zweiter Ord- 
nung (6.12) zu finden, die auf dem Rand des Meridianschnittes einen festen Wert besitzt. 
Da diese Aufgabe infolge der komplizierten Form dieser Differentialgleichung auBerordent- 
lich schwierig ist, soll hier eine Naherungsmethode fiir den Sonderfall entwickelt werden, 
daB die Schubfunktion g(i?) in der Form g(t2) = 1 + g, j darstellbar ist. Mit dieser Annahme 
vereinfacht sich (6.12) zu 


aD 30D , PD 1 {[ 4 /aP\? dD\?) 0? @D 7 [/ad\? 0D\?| 0D 
b= eee sarees ee (3) | a + (3) Sloe lee 


oDaD AD a@\2 aD\? aD) rus 
a Or Oz Or Oz (er) + 35) | Oz? | oy 


(6.16) 
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mit der Stoffzahl 
2 
=5a- (6.17) 


E 


Als Naherungsverfahren zur Lisung der Randwertaufgabe wenden wir die Fehlerquadrat- 
methode an. Fiir unser Problem wahlen wir den nur zweigliedrigen Ansatz 


P=pteg- (6.18) 
Die Funktion 
r4 Sir 35 74 


stellt die ersten drei Glieder der Potenzreihenentwicklung fiir die bekannte Liésung ! der linearen 
Theorie (¢ = 9) 


o* =— Ac i fiea\ oa 

java S\Wera) 3 
dar, wobei c eine noch offene Konstante bedeutet. Die Funktion g, erfiillt daher die Rand- 
bedingung ® = konst. am Rand. Von der Funktion gy, verlangen wir, daB sie am Rand des 
Meridianschnittes fir r = + kz verschwindet (k = tg a gesetzt), fiir r = 0 keine Spannungs- 
komponenten liefert und die aus ihr herzuleitenden Schubspannungen umgekehrt proportional 
zur neunten Potenz des Abstandes vom Koordinatenursprung sind. Die einfachste algebraische 
Funktion, die allen Forderungen geniigt, ist 


4 
1 =a (P—B), (6.20) 
Der Ansatz (6.18) befriedigt somit fiir beliebige Werte des Parameters a die Randbedingungen 


des Problems. 
Setzt man die Naherungsfunktion (6.18) in (6.16) ein, so erhalt man auf der rechten Seite 


eine Fehlerfunktion 
315 


8 
e(r, 234) = ey 2a(TB rt + 6 1222 — 55 ker? 2?) 


1 4 
cligae c3(945 78 + 70 r4 2? — 40 r? 24 + 16 28) (35 r*# — 20 r? 2 1. 8 24)2 (r? + 2), 


2 


r 
yi4 


we 


+ 3c ac? (227850 r¥4 — 875 (18 + 203 ) rl? 2 — 15 (2372 — 2065 2) 11° 24 
4. 1500(36 -- 11 2) r8 28 — 128(81 + 295 k2) r8 28 + 384(4 + 25 h2) r4 21° 
4192 (4 — 11 2) v2 al? — 256 be td] 
z 

4 6 ¢.a®c [44100 r!? + 45 (81 — 1540 h2) r1° 22 — (2022 — 30 k® — 27125 kf) r® zt 
+ 4(894 + 390 k® — 545 he) 16 28 + 24(18 — 214k? + 5 bf) 428 

2 b de 
— 32 (12 — 55 BY) ker? 9 4 64 be 232] 
8 ¢ a8 [8424 119 + 54(57 — 370 KY) 18 22 + 9(24 — 679 k® + 1750 kA) 19 24 
__ 3(108 — 1314 he £1375 b4) Rr 28 + 4(36 — 205 Be) heb v2 28 — 16 8 210] 

4 


Man bildet dann das iiber den ganzen Bereich erstreckte Flachenintegral 
s=l,r=ka 
P=2 f [ @@ > 23) drdz 
s=l, 0 
= Kj, ce? + Ky, ae+ Ky a? + ¢ (Ky, ct + Kj, @¢3 + Ky, a? c? + Ks, a? ¢ + Ky a?) 
+ 62(Ky, c& + K,, ac + Ky, a? c4 + Ky, a5 c3 + Kyy a4 c? + Ky, a ¢ + Key a°) 
mit den Abkiirzungen 


5 
Ke = (0,46154 + 0,96923 k2) k19(I-8 — I>8), 


4 
Ky == (4 + 11,758 h® + 24,111 bA) bOI — I), 


u Vel. S. Timoshenko, aa momes0. 
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1A = (682,67 — 2627,23 hk? + 4864 k4 +. 57935 k® — 102879 ks + 49026 k!° +. 907 858 12 
— 1108550 k4 + 897126 k6 +. 2058 000 k8) k9 (I-14 — J 14) , 
3 


Koy == (248,24 + 74,09 k? — 339,69 k4 + 5196 k8 + 7847 8 — 11767 kt + 36140 ki 
“+ 42 287 kl4) 411 (J-20 _ j—20) 
K,=5 (32,597 + 117,36 k? + 43,30 kt + 155,79 k8 + 1519,98 k® + 978,62 k9) ki8 (1-26 _ J—26) 
Kyo = 2(4,740 + 27,41 k? + 51,36 h4 + 38,52 hs) 19 (1-32 _ j-32) 
Kee — (42 366,71 — 313 931,19 k® + 1234 597 k4 — 505 817 k® — 9636316 k® + 41894457 j10 
~~ 48891 039 ki — 78 320 540 kt4 + 477 048 747 16 — 717522 332 kI8 + 350 142.547 420 
+ 984249 247 h22 — 1174753298 k24 + 445 482 838 k26 
+ 1205819263 k2) 11 (1-20 _ J—20) , 


(61 844,19 — 224151,46 k2 + 394517 k4 + 2088 566 k® — 7033375 k8 + 10046445 10 
+ 28978 336 k!2— 84.990 584 k14 + 124763 998 k6 + 65 607923 k}8 — 158227195 j,20 
+ 186 611076 k?2 + 217980398 kt) f13 (26 _ f 28), 


K.=— (23 402,90 — 18 134,95 k® — 49099,89 k4 + 768419 k® — 155 666 k8 — 1541618 f10 
+ 8917412 1? — 1121964 k14 — 3570448 k16 + 21010281 K18 
* + 15 083 696 ®) 45 (32 — J—s2) 


Kyo = 75 (2041,71 + 3735,88 k® — 4815,36 k4 + 34253 8 + 100516 k8 — 44678 410 
t 151320 k¥? + 486.479 h14 + 259262 hls) 27 (I-38 — J—38) , 


K,, = (64,165 + 289,46 k® + 278,20 k4 + 322,73 k® + 2798 k8 + 4256 10 
+ 1593 h12) R19 (1-44 —_ J—44) , 


ea = (9,880 + 68,76 k? + 170,98 k4 + 170,82 h8 + 86,59 KS) h2t (I-50 —_ J—50) . 


_ 


i) 


vy) 


Ky, = 


mel 


9 


Nun verlangt man, daB das mittlere Fehlerquadrat I méglichst klein wird. Die notwendige 
Bedingung fiir einen Kleinstwert dI/da = 0 fithrt auf die Bestimmungsgleichung fiir den 
Freiwert a 


6 6? Key a + 5 ee Ky, a4 + 4 (Ky) + € c? Ky) a3 l 
+ 3 ec(K*! + ¢c? K,5) a? + 2(Kyy + €¢? Kyo + 2 c4 Ky)a > (6.21) 
e(Ky, + ec Ky, + ct K,) =0. J 


Die Konstanie c bestimmt sich aus der Forderung, daB das Moment der im Wellenquer- 
schnitt wirkenden Schubkrafte das vorgeschriebene Torsionsmoment W ergeben muB: 


W = Barf tg, dr = Dek (1 — FB TE), 


woraus 
il M 
5 35 \ Qn 
Cola = = arp a eee 
i (1 3 as 16 i 


c= 


(6.22) 


folgt. 
Die aus dem Ansatz (6.18) heranleitenden Schubspannungen Tpr und t,, berechnen sich 
schlieBlich nach (6.10) zu 


r2 r 35 r4 r? ne 
Ta “(4 —10542 7) 4-207 (6? —5 k 2), 


r, 35 r4 


(6.23) 
r 
ee ee (4 10 +5 


ne 2a, (3 r? — 2 k2 22), 


gt 


Zahlenbeispiel: Kreiskegel mit 1, = 20cm, 1, = 50cm; halber Offnungswinkel 
tg a=—k= 0,2. Torsionsmoment W = 103 880 cm kg. Werkstoff: Kupfer mit G = 0,47- 106 
kg/cm’, g, = 7,26: 108; also nach (6.17) ¢ = 21,9: 10-8 em4/kg?. Nach (6.22) erhalten wir 
e = 11,03- 10° (cmkg). Die Bestimmungsgleichung fiir den Parameter a lautet nach (6.21) 
177+ 10-® a + 262- 10-64 a4 + 1617- 10-48 a3 
+ 4235 - 10~34 a? + 5525- 10-16 a + 1860 — 0. 
Die Auflésung ergibt a = — 0,4857- 1016 (cm? kg). 
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Wir ermitteln die Schubspannung t,, im Langsschnitt am Rand des Endquerschnittes 
(r =4cm,z = 20cm). Aus (6.23) erhalten wir Ty, = 848 kg/cm2. Die lineare Theorie (e = 0) 
liefert hier den Wert 1%, = 1000 kg/cm”. Der Wert der Schubspannung t,, am Rand des 
Endquerschnitts verringert sich demnach gegeniiber dem Wert der linearen Theorie bei 
gleichem Torsionsmoment um T5s2a: 


7. SchluBbemerkung. Soweit man aus den Zahlenbeispielen allgemeine Schliisse ziehen 
darf, kann man zweierlei folgern: einmal, da® der EinfluB der Nichtlinearitat des Elastizitats- 
gesetzes auf die Spannungen dann besonders groB ist, wenn das Problem nur von der Schub- 
funktion abhangt (Ziff. 6); und auBerdem, daf schon eine verhaltnismaBig kleine Abweichung 
des Elastizitatsgesetzes von der Linearitat die Spannungsspitzen erheblich abzubauen vermag 
(Ziff. 4b), sodaB also der in der technischen Praxis oft zu beobachtende Abbau der Spannungs- 
spitzen recht wohl auf geringfiigige Abweichungen vom Hookeschen Gesetz zuriickzufithren 
wire. Dies soll durch weitere Rechnungen noch nachgepriift werden. 


(Eingegangen am 17. Juli 1953.) 
Anschrift des Verfassers: Dr.-Ing. F. Jindra, Stuttgart-O, Gerokstr. 70. 
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Bemerkungen 


zu der Arbeit des Herrn Ki Yamada ..Vibration of turbine disc 
in its plane“. Proceedings of the Second Japan National Con- 
gress for Applied Mechanics 1952, S, 343. 


Von R. Grammel. 


In der angefiihrten Arbeit nimmt Herr Yamada ein 1935 von mir angeschnittenes Problem 
erfreulicherweise von neuem auf, aber leider enthalten seine Ausgangsgleichungen einen wesent- 
lichen Fehler. Herr Yamada schreibt die Bewegungsgleichungen fiir die radiale und azimutale 
Verschiebung g und ry eines Scheibenelements dm = (y/g)y rdrdp, an welchem die 
Spannungskrafte 

dP =d(yro,)dp— yo,drdg, 
dQ =d(yrt) dp + yr dr dp 
angreifen (wobei y/g, y, r, o,, o;, t und dp die iibliche Bedeutung haben), bei der mittleren 


Drehgeschwindigkeit w der Scheibe in der folgenden Form (mit Verbesserung einiger kleiner 
offensichtlicher Druckfehler) : 


gdm =—dP+ ra? dm+ 2rwywdm, (1) 
ry dm —dQ—2m@dm+ rwydm. (2) 


Die richtigen und (im Rahmen der Elastizitatstheorie der kleinen Verformungen) exakten 
Gleichungen aber lauten, wie ich seinerzeit angegeben habe}, 


odm —dP+ rw pdm =dP+rw2dm +2rowypdm- [ry dm], (3) 
rypdm =dQ— 2a + ) @dm = dQ— 2m 6dm— [2 yedm). (4) 


Die Richtigkeit meiner Gleichungen (3) und (4) sieht man rasch ein, wenn man sich entweder 
auf ein Bezugssystem bezieht, das die (veranderliche) Drehgeschwindigkeit w + w besitzt 
und in welchem als d’Alembertsche Krafte auftreten: die Fliehkraft r(@ + )?dm radial, sowie 
die Corioliskraft — 2(w + wp) @ dm azimutal und endlich die azimutale Tragheitskraft — r y dm; 
oder noch schneller, wenn man sich daran erinnert, daB (3) und (4) die bekannten klassischen 
Gleichungen der Zentralbewegung oder die Bewegungsgleichungen in Polarkoordinaten sind?. 
[Streng genommen mite man iberall in den Gleichungen (1) bis (4) genauer TF =r - @ 
statt r schreiben. Man schatzt aber ab, daB dies (wie schon in der statischen Spannungshe- 
rechnung, die man ebenfalls mit r statt 7 durchzufiihren pflegt) die Ergebnisse, also hier die 
Frequenzen, bei den iiblichen Baustoffen und Drehgeschwindigkeiten nur um einen Faktor 
(1 + «) verbessern wiirde, wo « héchstens von der GréBenordnung 10-4 bis 10-5 ware. Dies 
kann man ganz leicht schon an einem 1otierenden Stab abschatzen. | 

Herr Yamada beschrankt sich (allerdings ohne dies zu sagen) auf ,,kleine Schwingungen“, 
so daB man die in (3) und (4) eingeklammerten letzten Glieder weglassen darf, womit wenigstens 
die Yamadasche Gleichung (1) fiir ,,kleine‘’ Schwingungen richtig wird. Nicht verstandlich 
sind aber seine beiden Behauptungen, daf ich selbst in der Gleichung (3) die ,,Corioliskraft*‘ 
[also das letzte Glied in (1)] und in der Gleichung (4) ein Fliechkraftglied [also das letzte Glied 
in (2)],,vernachlassigt“‘ hatte. Beide Behauptungen sind unrichtig. 

DaB die erste Behauptung glattweg falsch ist, ersieht man an den einfach unter- 
strichenen Coriolisgliedern in (1) und (3). Diese unsinnige Behauptung ist umso erstaunlicher, 
als Herr Yamada lange Teile meiner Rechnung von 1935 wértlich iibernimmt, in denen gerade 
dieses Coriolisglied immer wieder explizit erscheint, das zusammen mit dem Coriolisglied 
—2wodm der Gleichung (4) das eigentliche Ker nstiick meiner ganzen Rechnung und ihrer 


1 R. Grammel, Ing.-Arch. 6 (1935) S. 258, Gl. (1). 
2 Vel. z. B. E. Meissner u. H. Ziegler, Mechanik Bd. 2, S. 164, Basel 1947. 
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Ergebnisse bildete. [Freilich tritt im System mit der Drehgeschwindigkeit w +- p keine Co- 
rioliskraft auf, aber ein zusatzliches Fliehkraftglied mit 2 rw y dm, und es ist bekannt, daB 
eben dieses Glied in einem System mit der Drehgeschwindigkeit w (wie es offenbar Herrn 
Yamada vorschwebt) das Coriolisglied ist. ] 

Die zweite Behauptung von Herrn Yamada ist deswegen falsch, weil es in der Glei- 
chung fiir » gar kein Flichkraftglied gibt [weder in einem System w noch 
in einem System w+ yp]. Das doppelt unterstrichene Glied in der Yamadaschen Glei- 
chung (2), das in meiner Gleichung (4) in der Tat nicht vorkommt, muB bei ihm weggelassen 
werden, DaB dieses Glied unmiglich ist, erkennt man beim Ubergang zum starren Korper; 
dabei wird 9 =0 und 1 = 0, und die Yamadasche Gleichung (2) nimmt die Gestalt ==)" p 
an, mit Lésungen von der Form yp = wy Go) w t, welche besagen wirden, daB eine starre Scheibe 
die mit einer Drehgeschwindigkeit @ + 0 umlauft, beim kleinsten zusatzlichen Drehsto8B wp 
,infolge der Fliehkraft* immer rascher und rascher umlaufen mite — ein verbliffend ein- 
faches Perpetuum mobile! 

Die allgemeine Lisung des Systems (3), (4) bereitet wegen der Verkoppelung der beiden 
Gleichungen durch die Drehgeschwindigkeit der Scheibe erhebliche Schwierigkeiten, sodaB 
bisher explizit nur asymptotische Lésungen gelungen sind1, welche Herr Yamada, durch sein 
falsches Flichkraftglied entstellt, wiederholt. 

Gliicklicherweise ist jene Verkopplung bei Dampfturbinenscheiben nur schwach 2, sodaB 
man wohl hinreichend genaue Lésungen erhalt, wenn man die radialen und die azimutalen 
Schwingungen je fiir sich an der nichtrotieren den Scheibe berechnet. Herr Yamada 
versucht diese Rechnung fiir ein spezielles (hyperbolisches) Scheibenprofil, obwohl er aus 
dem von ihm zitierten Band des ,,Ingenieur-Archivs wissen muBte, daB ich sie schon langst 
nicht nur fiir dieses spezielle, sondern fir beliebige Profile und beliebige Randbedingungen 
vollstandig erledigt hake®. 


(Eingegangen am 4, Marz 1954.) 
Anschrift des Verfassers: Professor Dr. R. Grammel, Stuttgart-N, Robert-Bosch- Str. 101. 


1 R. Grammel, a. a. O. S. 259. 
2 C, B. Biezeno u. R. Grammel, Technische D ik B i 

ie aes sche Dynamik Bd. 2, S. 35, 2. Aufl. Berlin 1953 (oder S. 655 
3 R. Grammel, Ing.-Arch. 6 (1935) 5. 442; Z. angew. Math. Mech. 5 (1925) 5, 193. 
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Der Kreisel. Seine Theorie und seine Anwendungen. Von Dr. R. Grammel, o. Pro- 


fessor an der Technischen Hochschule Stuttgart. Zweite, neubearbeitete Auflage. In 
zwei Banden.’ 


Erster Band: Die Theorie des Kreisels. mit 137 abbildungen. XI, 281 Sei- 
ten. Gr.-8°. 1950. DM 30.—; Ganzleinen DM 33.— 


. on ght Itstiber sicht: Einleitung. —-Grundlagen: Grundlagen der Vektorrechnung. — Grundlagen der Me- 
ge . — Der Tragheitstensor. —Der symmetrische Kreisel: Der kriftefreie symmetrische Kreisel. — 
ie gefiihrte Bewegung des symmetrischen Kreisels, — Der symmetrische Kreisel unter Zwang und StoB, — Der schwere 
symmetrische Kreisel..— Der EinfluG§ der Reibung. —Der unsymmetrische Kreisel: Der kriftefreie un- 
symmetrische Kreisel. — Der schwere unsymmetrische Kreisel. —Besondere Probleme: Kreiselim erweiterten 
Sinne. Gyroskopische Systeme. — Darstellung der Kreiselbeweg ungen durch Thetafunk- 
tionen, Namen- und Sachverzeichnis. ¢ , : 


) 


Zweiter Band: Die Anwendungen des Kreisels. mit 133 Abbildungen. 
VI, 268 Seiten. Gr.-8°. 1950. DM 30.—; Ganzleinen DM 33.— 


Inhaltsiitbersicht:; Kreiselwirkungen bei Radsatzen: Kollermtthlen — Kritische Drehzahlen 
von Rotoren — Fahrzeuge — Flugzeuge. —Kreiselgerate: Gerite mit KompaSkreiseln — Der KreiselkompaB — 
Kiinstliche Horizonte mit Pendelkreisen — Wendekreisel und Lagekreisel — Sonstige Kreiselgerite. —Unmittel- 
bare Stabilisatoren: Richtkreisel — Stiitzkreisel — Dimpfkreisel. — Namen- und Sachverzeichnis. 


»Das Werk stellt in seiner neuen Auflage nach Inhalt und Form eine sehr weitgehende Neuschépfung dar, deren Eigenart in 
einer kurzen Besprechung nur sehr unyollkommen gekennzeichnet werden kann. Besonders hervorzuheben ist der neu 
gestaltete Aufbau von den Grundlagen der Dynamik aus und die erzielte Klarheit und Anschaulichkeit des gesamten aus- 
gedehnten Stoffgebictes, bei der in vielen Einzelheiten neue Gesichtspunkte und Zusammenhange zum Ausdruck kommen. 
Ist es schon ein Beweis hoher Kunst, den groBen Stoffin dieser Vollstandigkeit in zwei schlanken Banden unterzubringen, 
so mu8 auch auf die Art der Darstellung, die Diskussion der erhaltenen Ergebnisse und deren Veranschaulichung durch 
Schaubilder und Kurventafeln besonders hingewiesen werden, ... (Archiv der Mathematik) 


$ 


Stabilititsprobleme der Elastostatik. you pst. rai. 


Alf Pfliiger, Professor an der Technischen Hochschule Hannover. Mit 389 Abbildungen. 
VIII, 339 Seiten. Gr.-8°. 1950. Ganzleinen DM 34.50 


X\ 
Inhaltsibersicht: Crundsitzliches tiber Stabilititsprobleme. — Exakte Lésungen. — Kriterien fiir die Gleich- 
gewichtsarten.— Zwei-und dreidimensionale Probleme. — Naherungsverfahren fiir Verzweigungsprobleme. — Giiltigkeits- 
grenzen des Naherungsverfahrens. — Naherungslésungen fiir Kigenwertprobleme. — Praktische Bedeutung der Stabilitats- 
ieee Sea Formelzusammenstellung fiir kritische Werte von Verzweigungsproblemen. — Literaturverzeichnis. 
— Sachverzeichnis. 


y+ «- Das, Buch zeichnet sich durch einen auBerordentlich anschaulichen und sorgfaltigen Aufbau aus. Fiir den Praktiker 
wird besonders der letzte Abschnitt des Buches bedeutungsvoll sein, in dem von Widerspriichen zwischen Theorie und 
Versuch, von Abweichungen vom Hookeschen Gesetz und auch vom Wert der Theorie fiir die Anwendung gesprochen wird. .‘ 

‘ (Forschung auf dem Gebiete des Ingenieurwesens.) 
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Dynamik des Bogentragers und Kreisringes. ., 


Dr. Karl Federhofer, Professor an der Technischen Hochschule Graz. Mit 35 Textabbil- 
dungen und 26 Zahlentafeln. XII, 179 Seiten. Gr.-8°. 1950(W) Steifgeheftet DM 23,.— 


Inhaltsibersicht: Kinematik der Biegung und Drillung des von Haus aus gekriimmten und tordierten Stabes. — 


Die Differentialgleichungen der Eigenschwingungen des réumlich schwingenden offenen Kreisbogens und ihre Rand-. 


bedingungen, — Sonderfall des Kreisbogentragers mit einfach-symmetrischem Querschnitt, dessen eine Hauptachse in der 
Kreisebene liegt. — Einflu8 der Querschnittsform auf die Schwingungsgleichungen in den Sonderfillen des vorhergehenden 
Abschnittes. — Der Grenzfall des geraden Stabes. — Kreisbogentriger mit entlang der Bogenachse verinderlichem Quer- 
schnitte. — EinfluB eines radialen AuBen- oder Innendruckes auf die Frequenzen eines Kreisbogentrigers und Kreisringes. 
— Der Kreisring mit diinnwandigem geschlossenem Hohlquerschnitt (Hohlreifen). — Frequenzen des Zweigelenk-Bogen- 
trigers mit parabolischer Achse. 


+» » Das yorliegende Werk des bekannten Autors stellt im wesentlichen cine Erweiterung unserer Kenntnisse in der Schwin- 
gungsforschung dar, Wird doch im Gegensatz zu den bekannten Standardwerken tiber technische Schwingungsprobleme 
der Fall gekrimmter Stabe einer éingehenden Behandlung unterzogen. Der Verfasser beschrankt sich aber nicht nur auf die 
allgemeine Behandlung der Probleme, sondern er ist bemiht, bis zur zahlenmaBigen Auswertung jedes Problems vorzu- 
dringen. So findet man auch die Naherungsverfahren in vollkommener Art durchgefiihrt, und es ist dem Leser méglich, 
mittels der zahlreichen Tabellen und Darstellungen jedes einschlagige Problem auch rasch tiberschlagsweise zu beurteilen.“ 

(Acta Physica Austrica.) 
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SPRINGER-VERLAG IN WIENI 


Die Methode der Festpunkte. Von Ernst Suter, Vorelabactree Ver- 


fahren zur Berechnung statisch unbestimmter Konstruktionen mit Beispielen aus der 
Praxis; insbesondere von Stahlbetontragwerken. Dritte, neubearbeitete Auflage von 
Dipl.-Ing. Ernst Traub. Mit 232 Abbildungen und 7 Tafeln im Text. XII, 216 Seiten. 
Gr.-8°, 1951. Ganzleinen DM 21.— 


Inhaltsuibersicht: Einleitung. Allgemeines. Die Mohrschen Satze. Graphische Bestimmung der Festpunkte fir 
durchlaufende Trager. —Erster Abschnitt: Tragwerke mit unverschieblichen Knotenpunkten und mit von Stab 
au Stab verinderlichem, aber auf Stablange konstantem Tragheitsmoment: Das k-Verfahren zur Bestimmung der Fest- 
punkte, VerteilungsmaBe und Ubergangszahlen, Das k-Verfahren zur Bestimmung der Kreuzlinienabschnitte und Momente 
infolge beliebiger Belastung des Tragwerkes. Rechnungsabschnitt I. Ubersichtliche Zusammenstellung des 
Rechnungsvorganges. Bestimmung der Querkrafte, Normalkrafte und Auflagerkrafte. Bestimmung der Festhaltungs- 
krifte. — Zweiter Abschnitt: Rechnungsabschnitt I. Tragwerke mit verschiebbaren Knuoten- 
punkten und mit von Stab zu Stab veranderlichem, aber auf Stablinge konstantem Dyagheitsmoment: Biegungsmomente 
infolge gegenseitiger Verschiebung der beiden Endpunkte eines Stabes und die zugehorige Erzeugungskraft, Tragwerke mit 
senkrechten Stielen und waagerechten Balken, Tragwerke mit schiefen Stielen oder schrigen Balken. —Dritter Ab» 
schnitt: Tragwerke mit auf Stablinge veranderlichem Tragheitsmoment: Das k-Verfahren zur Bestimmung der Fest- 
punkte und Kreuzlinienabschnitte. Das Drehwinkel-Verfahrea zur Bestimmung der Festpunkte und Kreuzlinienab- 
schnitte. —Vierter Abschnitt: Grenzwerte der Momente und Querkrafte: Grenzwerte der Momente‘und Quer; 
krafte fiir standige Last und Nutzlast. EinfluBlinien der Momente und Querkrafte fiir bewegliche Lasten. Momente 
infolge Temperaturanderung und Stiitzensenkung. — Finfter Abschnitt: Hilfstafeln zur Bestimmung der 
Festpunkte und Kreuzlinienabschnitte, Momententafeln fiir durchlaufende Trager. — Sechster Abschnitt: 
Beispiele aus der Praxis. . 
,,.., Dem Studierenden wird in klar aufgebauter Darstellung, ohne Wiederholungen, alles Notwendige gesagt, wobei einige 


ganz durchgerechnete ,Beispiele aus der Praxis’ den Eingang in die praktische Anwendung der Methode erleichtern.... Be 
Das Buch kann jedem Rahmeustatiker und solchen, die es werdén wollen, wirmstens empfohlen werden.‘** é. 
. : ~ (Schweiz. Technische Zeitschrift.) % 

Die Metheden der Rah alc 
ie Metheden der Rahmenstatik. augax, Zusammenfassung | 
und Kritik. Von Dr.-Ing. habil. Otto Luetkens. Mit 38 Abbildungen und 9 Zahlentafeln. A 
VII, 281 Seiten. Gr.-8°. 1949. DM 33.—; Ganzleinen DM 36.— be 4 

4 fon 

Inhaltsiibersicht: Linleitung. —I. Grundlagen der Berechnung: A. Einfihrung einheitlicher Bezeichnungen. 2 
B. Beziehung zwischen Kraften und Formanderung. C, Aufstellung der Elastizitatsgleichungen. D. Auflésung-der Elasti- a 
zititsgleichungen. E. Umordnung. — II. Beschreibung der statischen Verfahren: A. Daten des Systems der Vergleichs- 3 
rechnung. B. Verfahren am Hauptsystem A. -C. Verfahren am Hauptsystem B, D. Verfahren am Hauptsystem C. — ri 


Ill. Eignung*der Verfahren: A. Betrachtung der Grundlagen. B. SchluSfolgerung auf die Berechnung der einzelnen Arten 
von Systemen, — IV. Literaturverzeichnis. — V. Anhang, Tafeln Avl bis-A 9. 


Die vorliegende Arbeit bezweckt nicht nur cine Zusammenfassung der statischen Verfahren, sondern auch einen Ausgleich 
zwischen den beiden Anschauungen, wonach die Rahmenstatik einerseitsin der Lésung algebraischer Gleichungen, anderer- 
seits in der Auslegung statisch-geometrischer Zusammenhiange besteht. 

Der Umfang der Arbeit beschrinkt sich bewuBt auf die einfachen Rahmensysteme mit geraden Staiben von stabweise 
konstanten Triagheitsmomenten und damit auf das praktische Bedirfnis fiir die Stabwerke im Hochbaun. Der Hauptwert 
wurde nicht auf die Entwicklung y neuer‘ Methoden, sondern auf die gedankliche Zusammenfigung vorhandener Verfahren 
gelegt. Am SchluB dieses Buches befindet sich ein umfangreiches Schrifttumsverzeichnis. ay 


Statik der rahmenartigen Tragwerke. you pros. Drang. 


J. Pirlet, Kéln, ehem. Honorarprofessor der Technischen Hochschule Aachen, Mit 80 Ab- 
bildungen im Text und 5 Tafeln in 1 Tasche. VII, 168 Seiten. 8°. 1951. Canzl. DM 24,— 


Inhaltstibersicht: Erster Teil: Grundlagen der Berechnung rahmenartiger Trag- 
werke. Der elastisch eingespannte Stab. Das Stabgefiige. Der Knoten, Der Knotenstabzug. Die Einspanngrade e, 
Berechnung der Einspannmomente M; und Mj, d. h. der Unbekannten in einem Stabzugfelde 1;. Gleichungen fiir M; und 


My;, bei einzelnen Belastungsarten. — Anhang: Zusammenstellungen der Ergebnisse. — Zweiter Teil: Der 


elastisch eingespannte kontinuierliche Balken in rahmenartigen Tragwerkem — 
Das Berechnungsyverfahren. Einleitung. Die Einspanngradee. Die Einspannmomente M; und Mj, des 


beiderseits elastisch eingespannten Balkens, Die Einspannmomente M in kontinuierlichen Trigern. —Dritter Te il: 
Zah ie nbeispiel. Von der Belastung unabhingige Werte. Von der Belastung abhingige Werte. -Sachver- 
zeichnis. : 


» Der Verfasser bringt in seinem itibersichtlich gehaltenen Buch ein Berechnungsverfahren fiir rahmenartige Tragwerke, das 
geeignet ist, dem praktischen Statiker die Arbeit bei der Bemessung dieser vielfach statisch unbestimmten Systeme wesent- 
lich zu erleichtern. Hier wird ein Verfahren entwickelt, das zunichst im geschlossener Form genaue Ergebnisse unter der 
Voraussetzung der Unverschieblichkeit der Knotenpunkte und konstanter Querschnitte liefert. Es geht demnach im 
wesentlichen um den elastisch eingespannten Balken bzw. um den Durchlauftrager mit elastischer Einspannung dér Knoten 
Die gesuchten Momente werden als Funktion der ,,Einspanngrade“ der Stahe gefunden. In dem Umstand, da das Ziel dee 
Lésung die EinfluBlinien sind, liegt der besondere Wert der Arbeit fiir die Praxis. : amy 

AnschlieBend an das strenge Verfahren wird eine Naherungsberechnung entwickelt und an Hand von Beispielen nach- 
pe werner o deren Ergebnisse fiir die meisten praktischen Fille genau genug mit der exakten Lésung iibereinstimmen. 
ke aeee Pechan ate ame eine groBere Anzahl von Formeln fir Sonderbelastungsfalle und Zahlenbeispiele eqletohtcrs 
Die Art der behandelten Tragwerke bringt es mit sich, da8 das Buch zunichst in erster Linie den Sta i - 
spricht. Da das Verfahren aber geeignet ist, als Grundlage fiir die Berechnung der eigentlichen Harenapiar se Hi 
und der Verfasser eine spitere Erweiterung seines Werkes in dieser Richtung in Aussicht stellt, ist die Arbeit eine willkommene 
Bereicherung des Riistzeuges fiir die gesamte Statik." (Zeitschrift des dsterreichischen Ingenieur- u. Architekten-Vereins) 
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EinfluBfelder elastischer Platten, ,. pipi-tng. prof, Dr. teenn, 


Adolf Pucher, Graz. 52 Tafeln mit VIIT, 13 Seiten Text und 10 Textabbildungen. Quer-4°. 
1951. (W) Ganzleinen DM 27.70 


y+. Das vorliegende Werk diirfte cine begriiBenswerte Bereicherun i i ir di 

: : c g des Hilfsmittel fiir die schnell ; 
isotropen Platten sein und wird u. a. bei der nach der neuen DIN 1075 geforderten Berechnung a Pabcboheseuarateh 
der Plattentheorie wertvolle Dienste leisten. Das gut ausgestattete Buch kann daher bestens empfohlen wenden 2 


(Der Bauingenieur) 
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